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I.   INTRODUCTION 

As  in  many  other  regions  of  the  country,  the  supply  of 
natural  gas  in  Montana  may  not  be  sufficient  to  meet  future 
needs.   One  method  to  avoid  this  problem  in  Montana  is  to 
produce  SNG  (synthetic  natural  gas)  using  the  state's  exten- 
sive coal  reserves.   The  purpose  of  this  study  is  to  deter- 
mine the  demand  for  SNG  in  Montana  through  2000  A.D. 

Determining  the  demand  for  SNG  essentially  is  a  balancing 
problem.   The  SNG  demand  is  the  additional  production  required 
to  balance  natural  gas  requirements  with  the  supplies  avail- 
able.  Thus,  the  demand  for  and  supply  of  natural  gas  must 
be  determined  before  the  demand  for  SNG  can  be  found.   Natural 
gas  cannot  be  considered  alone  as  it  competes  with  other  fuels. 
Major  conversion  of  industry  gas  uses  to  coal  and  electricity 
certainly  are  possible  and  can  affect  the  demand  for  natural 
gas  significantly.   The  cost  of  producing  SNG  is  likely  to 
be  different  from  the  natural  gas  obtained  from  conventional 
sources.   Thus,  SNG  may  enter  the  market  place  and  compete 
at  different  price  levels  than  other  natural  gas. 
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The  approach  used  here  was  to  develop  a  model  of  Mon- 
tana's economy  on  which  energy  demands  can  be  based.   The 
complete  documentation  of  the  model  is  contained  in  appendices 
A-D  of  this  report.   The  industrial  and  commercial  demands 
for  energy  are  based  upon  the  levels  of  activity  in  the  corre- 
sponding sectors  of  the  economy  while  residential  demands  are 
based  upon  the  number  of  households  and  their  income.   In  addi- 
tion, the  effects  of  inter-fuel  competition  and  conservation 
are  considered. 
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The  primary  modeling  technique  used  is  dynamic  systems 
simulation.   This  methodology  deals  with  the  forces  causing 
change  in  the  economy  and  should  be  well  suited  for  the  tran- 
sients which  are  likely  in  Montana's  future.   It  also  allows 
the  model  to  be  based  more  closely  on  the  actual  processes 
and  actions  which  occur  in  Montana's  economy.   The  dynamic 
simulation  is  combined  with  input/output  analysis  in  deter- 
mining levels  of  economic  activity.   The  use  of  the  input/ 
output  analysis  assists  in  obtaining  more  detail  in  this 
critical  area. 


Once  developed,  the  model  was  used  to  analyze  the  effects 
of  inter-fuel  substitution  and  various  energy  conservation 
scenarios.   Some  industries  in  Montana  recently  switched  from 
gas  to  other  fuels,  and  more  substitutions  are  projected  for 
the  future.   Residential  heating  is  a  major  gas  use  in  Montana; 
conservation  in  that  area  was  modeled  for  several  projected 
gas  price  levels.   Fuel  substitution  and  energy  conservation 
will  be  significant  factors  in  the  gas  demand  picture  for 
Montana. 
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II.   DYNAMIC  SYSTEM  SIMULATION  MODEL  OF  MONTANA'S  ECONOMY 

The  demand  for  natural  gas  (including  SNG)  in  Montana 
is  determined  by  the  level  of  activity  in  the  Montana  economy. 
The  SNG  study  group  felt  that  it  would  be  worthwhile  to  de- 
velop a  computer  model  of  the  economy  to  aid  in  demand  analy- 

i 
sis.   A  dynamic  systems  simulation  approach  was  selected  as 

most  appropriate  for  the  model.   This  approach  takes  into 

account  the  various  interrelationships  and  feedback  effects 

in  the  economy. 

The  model  is  composed  of  several  sub-models  each  dealing 
with  a  certain  area  of  the  economy.   Overall  structure  for 
this  model  is  illustrated  in  Figure  1.   The  model  is  similar 
in  concept  to  the  Oregon  State  Simulation  Model,  and  the 
structure  of  some  of  the  sub-models  basically  is  the  same  as 
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thexr  Oregon  counterparts.   The  Montana  Futures  Project  input/ 
output  model  also  was  used  in  the  formulation  of  the  model. 
Each  sub-model  is  discussed  separately  in  the  following  sec- 
tions . 


A.    Demographic  Sub-Model 

The  demographic  sub-model  predicts  both  the  number  of 
people  that  will  be  in  Montana  and  the  number  of  households 
in  the  state.   Projections  are  based  upon  current  population 
and  predicted  net  migration.   Migration  is  affected  by  the 
difference  between  Montana  unemployment  levels  and  national 
unemployment  levels. 


to 


o 


cu 

E 

ZS 


"O  -I- 

III  to 

on  3 

=3  O 

■o 

c  5- 

(O  o 


Demographic  Sub-Model 
^ — 


Land  Use 
Sub-Model 


Agricultural    & 
Forest  Land 
in  Production 


o. 

Q 
CO 

s_ 
o 

jO 


o 

Q. 

E 


.aland  Mined. 


to 
<u 
+-> 
<a 
or 

cu 

CD 


o 

■r— 

o 


0  National  Unemployment 
•o  National  Wage  Rates 


1 


£ 


Prices  for 
° Primary  Products 

Export  &  Gov't  Demand 


Economic  Sub-Model 


Land  "Urbanized" 


Energy  Sub-Model 


T 


T 


•6 


SNG 
Production 


Imports 
Available 


Energy 
Prices 


+-> 

o 

<*- 

o 


00 
OJ 
CD 

a 
o 

CO 


Figure  1. --Overview  of  Model  Structure 
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Land  Use  Sub-Model 


The  land-use  component  maintains  a  current  inventory  of 
land  in  each  of  six  categories.   Residential  land  demand  is 
a  function  of  the  rate  of  new  household  formation  and  the 
residential  conversion  density  (which  originates  in  the  demo- 
graphic component).  Demand  for  commercial/industrial  land  is 
determined  by  the  rate  of  employment  and  the  employees/acre. 
This  is  generated  by  activity  levels  in  the  industrial  sectors 
of  the  economic  component.   The  figure  for  total  land  converted 
from  agricultural  to  urban  land  is  obtained  from  the  sum  of 
the  demand  for  residential  and  commercial/industrial  land. 


Mining  land  is  converted  from  agricultural  land  by  a 
scenario  variable  derived  from  standard  coal  production  per 
year.   Waste  land  is  derived  by  taking  the  product  of  the 
land  mined  and  the  fraction  of  revegetation  (this  fraction 
being  a  technical  policy  input) . 

Forest  land  accounts  for  all  the  land  upon  which  com- 
mercially harvestable  timber  is  growing.   Agricultural  land 
can  be  converted  to  forest  land;  however,  this  is  unlikely 
in  Montana.   Forest  land  may  be  converted  to  wilderness  land 
and  vice  versa  by  policy  input. 

C .    Economic  Sub-Model 

Five  separate  models  are  used  for  the  primary  economic 
activities:   agriculture,  forestry,  coal  mining,  oil  and  gas 
production,  and  non-energy  mining.   Since  the  output  of  the 
petroleum  refining  and  primary  metal  processes  depends  mainly 
on  the  amount  of  Montana  oil  produced  plus  crude  oil  imports 
and  metal  mined,  they  also  are  considered  separately.   This 
part  of  the  model  can  be  considered  supply  limited.   The  out- 
put of  these  sectors  plus  consumer  expenditures,  export  de- 
mands, government  spending,  and  capital  expenditures,  drive 
an  input/output  model  for  the  rest  of  the  economy.   These 
activities  can  be  considered  demand  limited.   The  outputs 
then  are  used  to  determine  total  employment  and  income. 

Each  of  the  models  for  the  primary  activities  and  for 
petroleum  refining  and  primary  metal  processing  will  be  dis- 
cussed separately.   The  capital  investment  and  employment 


models  also  are  discussed  separately. 

D.    Agricultural  Sub-Model 

The  agricultural  sub-component  of  the  economic  model  is 
divided  into  two  parts,  agricultural  production  and  capital 
formation. 

Agricultural  Production:   The  amount  of  agricultural 

land  and  the  level  of  capital  employed  are  combined  (utili- 

2 
zing  a  modified  Cobb-Douglas  production  function)  to  generate 

the  production  level.   The  gross  agricultural  production  is 

converted  to  net  farm  receipts  by  subtracting  the  aggregate 

expenses  incurred  in  production.   The  production  expenses 

fall  into  four  categories:   direct  operating  expenses,  labor 

expense,  interest  expense,  and  property  tax  expense. 

Capital  Formation:   The  agricultural  capital  formation 
model  starts  with  the  value  of  net  farm  receipts  developed 
in  agricultural  production  and  determines  the  resulting  change 
in  capital.   If  net  farm  receipts  are  greater  than  zero,  a 
portion  is  used  to  purchase  new  capital.   Depreciation  of 
capital  stock  is  assumed  to  occur  at  a  rate  that  is  propor- 
tional to  the  installed  level. 


In  the  event  that  net  farm  receipts  are  less  than  zero, 
liquidation  can  occur.   If  the  magnitude  of  the  loss  exceeds 
a  specific  fraction  of  gross  production,  it  is  assumed  that 
the  agricultural  community  will  be  forced  to  further  reduce 
the  level  of  owned  capital  stock  by  liquidating  it  at  its 


salvage  value;  otherwise,  capital  is  simply  allowed  to  depre- 
ciate without  replacement. 

E.  Logging  Sub-Model 

The  principal  assumption  in  the  logging  sub-model  is 
that  the  logging  sector  of  Montana's  economy  in  the  long  run 
is  driven  by  the  projected  inventory  of  standing  mature  saw- 
timber  available  to  be  cut  in  each  year.   The  total  capital 
to  harvest  the  allowed  timber  cut  is  adjusted  by  technology 
and  the  future  capital  required.   Additional  capital  is  a 
relationship  between  the  total  capital  required  and  the 
present  level  of  capital. 

The  annual  cutting  rate  is  determined  by  calculating 
the  difference  between  the  standing  timber  available  to  cut 
and  the  maximum  possible  cut  given  the  current  level  of 
installed  capital. 

F .  Coal  Mining/Non-Energy  Mining  Sub-Model 

The  models  for  coal  mining  and  non-energy  mining  are 
identical  in  structure.   Only  the  parameters  are  different. 
For  the  sake  of  simplicity,  they  will  be  discussed  together 
here . 

The  productivity  and  cost  for  mining  are  determined  by 
relating  mining  cost  and  capital  cost  functions  to  the  total 
quantity  that  has  been  mined.   A  cost  distribution  function 
determines  what  fraction  of  total  possible  output  is  economi- 
cal.  This  is  used  to  determine  the  actual  output.   Profit 
rate  and  the  rate  of  return  are  calculated  by  subtracting 
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mining  expenses  from  mining  income. 

G.    Oil  and  Gas  Sub-Model 

The  oil  and  gas  model  consists  of  parallel,  coupled 
models — one  for  oil  discovery  and  production  and  one  for  gas 
discovery  and  production.   The  models  are  coupled  since  oil 
production  normally  has  some  gas  production  associated  with 
it  and  production  from  gas  fields  usually  results  in  the  pro- 
duction of  some  liquids  as  well.   Similarly,  when  exploration 
for  oil  fields  is  carried  out,  there  is  some  probability  of 
finding  a  gas  field  (or  vice  versa) .   The  following  discussion 
deals  only  with  oil  exploration  and  production  since  the 
structure  is  the  same  for  gas. 

When  oil  is  discovered,  it  is  placed  in  the  producing 
level  and  information  is  retained  as  to  the  age  of  each  dis- 
covery.  Four  functions  are  required  to  describe  the  perfor- 
mance of  producing  reservoirs. 

The  production  function  traces  a  typical  reservoir 
through  the  various  stages  of  discovery,  steady-state  pro- 
duction, and  finally,  a  long  period  of  steady  decline  to  de- 
pletion.  The  operating  expense  function  shows  the  expendi- 
tures required  to  operate  a  producing  field  of  a  given  age. 
The  investment  function  gives  the  rate  at  which  investment 
must  be  made  in  order  to  develop  and  maintain  a  field. 

The  variable  cost  of  production  is  similar  to  the 
operating  expense  function  except  that  it  describes  the  cost 
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on  the  basis  of  unit  of  production  for  fields  of  each  age. 
This  function  is  used  to  compare  the  price  to  the  variable 
cost. 

The  total  oil  production  is  the  sum  of  the  oil  pro- 
duced in  oil  fields  plus  the  liquids  produced  in  natural 
gas  fields.   The  total  income  for  oil  production  is  then  the 
oil  produced  in  oil  fields  multiplied  by  the  oil  price--plus 
the  gas  produced  in  the  oil  fields  multiplied  by  the  wellhead 
gas  price.   The  expenses  for  the  production  come  from  the 
total  operating  expenses  plus  the  costs  for  capital.   The  net 
profit  rate  per  barrel  of  oil  is  determined  by  subtracting 
the  expenses  from  the  income  and  dividing  by  the  production 
rate. 

The  rest  of  the  model  deals  with  exploration.   Two 
additional  functions  are  required  here.   The  depletion 
function  gives  the  increase  in  drilling  required  for  addi- 
tional discoveries  as  more  and  more  oil  is  discovered.   The 
investment  function  gives  the  rate  of  investment  in  terms  of 
the  profit  rate. 

H .    Petroleum  Refining  Sub-Model 

Montana's  petroleum  refining  is  very  closely  related 
to  the  output  of  the  petroleum  sector.   It  also  is  dependent 
upon  imports  of  crude  oil.   These  factors  are  more  important 
than  actual  demand  for  refinery  products  in  Montana.   There- 
fore, it  must  be  considered  separately  rather  than  with  the 
economic  output  model. 
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The  vast  majority  of  the  refinery  capacity  in  Montana 
is  located  around  the  Billings  area.   However,  the  pipeline 
system  is  such  that  the  oil  produced  in  the  eastern  part  of 
the  state  is  not  available  to  these  refineries. 

In  1975  81%  of  the  total  oil  refined  in  Montana  was 
imported  crude  oil.   However,  it  is  difficult  to  predict  the 
availability  of  crude  oil  imports  in  the  future,  especially 
Canadian  imports.   Thus,  the  refinery  output  will  depend  upon 
the  available  Montana  production  plus  the  available  imports 
or  the  maximum  output  of  current  capacity. 

I .    Primary  Metals  Sub-Model 

The  output  of  the  primary  metals  sector  of  the  economy 
is  dependent  primarily  upon  the  output  of  the  metal  mining 
sector.   It  also  is  very  energy  intensive.   Primary  metals 
output  is  assumed  to  be  proportional  to  the  non-energy 
mining  output. 

J.    Economic  Output  Model 

The  economic  output  model  is  simply  an  input/output 
model  for  Montana.   Outputs  for  agriculture,  forestry,  coal 
mining,  oil  and  gas,  non-energy  mining,  petroleum  refining, 
and  primary  metals  are  all  determined  elsewhere.   The  direct 
input/output  coefficients  are  used  to  determine  the  demands 
for  other  products  due  to  these  outputs.   These  demands  then 
are  treated  as  demands  exogenous  to  the  input/output  model. 
Additional  exogenous  demands  arise  from  consumer  spending, 
exports,  and  government  spending. 
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K.    Capital  Investment  Sub-Model 

The  capital  investment  for  the  primary  sectors  is  deter- 
mined in  their  individual  models.   The  same  capital  invest- 
ment structure  is  used  for  all  non-primary  sectors. 

The  model  uses  exponential  delays  to  average  both  the 
output  and  the  rate  of  change  for  a  sector.   The  averaged 
rate  of  change  is  added  to  the  depreciation  rate  to  get  the 
total  investment  rate  required.   However,  if  this  rate  is 
negative,  there  is  no  investment.   If  it  is  positive,  it  is 
multiplied  by  the  averaged  output  and  the  capital  output 
ratio  to  determine  the  actual  investment. 


L .    Employment  Sub-Model 

Employment  is  divided  into  three  groups:   agricultural, 
non-agricultural  primary,  and  derivative  employment.   The 
primary  reason  this  division  is  made  is  because  wage  rates 
for  the  different  groups  vary  considerably.   Levels  of  em- 
ployment are  determined  from  economic  output  and  employee/ 
output  rates.   Government  employment  should  relate  to  govern- 
ment spending  in  much  the  same  way  that  private  employment 
relates  to  output. 

The  employment  for  the  various  sectors  are  summed  to 
determine  the  total  employment.   The  total  work  force  avail- 
able is  determined  by  multiplying  the  total  population  by 
the  employee/population  projection.   The  total  employment  and 
the  total  work  force  are  compared  to  determine  the  unemploy- 
ment rate. 
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M .    Energy  Sub-Model 

The  energy  sub-model  is  divided  into  a  fairly  detailed 
demand  section  and  a  simple  supply  section  that  obtains  much 
of  its  information  from  the  economic  sub-model. 

The  demand  for  energy  is  analyzed  along  the  traditional 
lines  of  residential,  commercial,  and  industrial  usage.   Each 
part  of  the  demand  is  treated  somewhat  differently  and  is 
discussed  separately. 

Residential  Demand The  demand  for  residential  energy 

is  divided  among  its  various  end  uses.   These  uses  include 

the  following  factors: 

Space  Heating 
Water  Heating 
Kitchen  Ranges 
Refrigeration 
Clothes  Drying 
Air  Conditioning 
Other 

The  actual  residential  demand  for  a  particular  form  of 

energy  by  a  particular  end  use  is  determined  by  the  following 

relationship: 

E .   :  =  A .   .XS.   .  x  TH 
!»]     1/3     i,D 

where 

E^  ■    is  the  consumption  of  energy  type  i  by  end  use  j , 

A^  •  is  the  average  consumption  of  energy  for  this  type 
of  end  use  in  a  typical  household, 

S.  .  is  the  saturation  of  type  j  end  use  using  energy 
,J  type  i,  and 

TH  is  the  total  number  of  households. 
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Industrial  Demand Industrial  energy  demand  is  deter- 
mined separately  for  each  sector.   The  use  of  a  particular 
type  of  energy  in  a  given  sector  is  determined  by  the 
following  relationship: 


Ki.j  -  Ai,j  *  °j 


whe  re 


E.  .is  the  use  of  energy  type  i  in  sector  j, 

A-  •  is  the  energy  type  i  consumption  per  unit  of  out- 
put,  and 

0  •  is  the  output  of  sector  j . 

Commercial  Demand There  appears  to  be  very  little 

data  available  which  gives  details  on  commercial  energy 
demand.   This  is  a  common  problem  for  all  energy  studies. 
Usually  everything  that  is  left  over  is  lumped  into  the 
commercial  category. 

The  output  for  the  commercial  sector  also  was  treated 
as  a  residual.   It  consists  of  all  trade,  service,  and  govern- 
ment sectors;  the  construction  sector;  and  the  communications 
sector.   Total  output  for  these  sectors  was  1,924  million 
dollars  in  1974.   The  Montana  Business  Quarterly  reports  that 
income  in  Montana  increased  approximately  5%  from  1972  to 
1974.   It  was  assumed  that  commercial  output  increased  simi- 
larly.  An  adjusted  output  of  2,021  million  dollars  then  was 
calculated.   This  value  for  output  was  used  to  calculate  the 
energy  coefficients. 
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N.    Energy  Supply 

The  energy  supply  part  of  the  model  is  relatively  simple 
compared  to  the  energy  demand.   Most  of  the  actual  supply 
modeling  is  done  in  the  economic  part  of  the  model  which 
gives  the  total  production  of  gas,  oil,  and  coal.   However, 
most  of  the  coal  production  is  destined  for  export,  and  a 
significant  part  of  the  oil  and  gas  production  is  not  avail- 
able for  Montana  use. 

0.    Energy  Conservation  in  the  Residential  Sector 

The  amount  of  natural  gas  used  in  the  residential  sector 
could  change  considerably  as  prices  increase.   There  may  be 
reductions  in  energy  consumption  due  to  fuel  shifts,  invest- 
ments in  home  heating  and  energy  saving  devices,  or  merely 
by  an  adjustment  in  living  habits.   These  changes  could  af- 
fect the  total  demand  for  natural  gas  significantly. 

Our  study  of  residential  energy  conservation  focused  on 
home  heating  since  it  dominates  natural  gas  consumption  in 
this  sector.   The  substitution  of  other  fuels  for  natural  gas 
was  eliminated  since  other  fuels  are  not  likely  to  be  cost 
competitive  to  the  point  of  causing  major  fuel  shifts  in 
home  heating  even  when  price  increases  are  considered.   How- 
ever, existing  methods  for  reducing  heating  requirements 
through  home  insulation  were  found  to  be  cost  effective. 


Energy  savings  and  cost  data  can  be  used  to  determine 
the  return  for  installing  insulation  for  various  price 
scenarios  in  terms  of  a  payback  period. 

14 


A  simple  conservation  model  simulates  the  installation 
of  insulation  and  the  impact  it  has  on  natural  gas  for 
heating  a  typical  Montana  home. 

In  addition  to  the  reductions  resulting  from  insulation, 
significant  energy  saving  could  result  from  modified  living 
habits  due  to  high  natural  gas  prices. 
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III.   SCENARIO  PRESENTATION 

A.    Baseline  Gas  Demand  in  Montana 

A  baseline  gas  demand  in  Montana  was  established  for 
the  period  between  now  and  the  year  2000  (see  Figure  2) . 
The  baseline  demand  assumes  very  little  home  energy  conser- 
vation.  Industrial  fuel  substitution  is  assumed  to  be  mini- 
mal also.   This  baseline  represents  a  continuation  of  current 
policies  and  practices. 

The  demand  for  natural  gas  is  not  considered  to  be  supply 
limited.   Significant  overall  economic  growth  is  anticipated 
in  Montana,  but  not  all  sectors  of  the  economy  will  experience 
the  growth.   The  petroleum  industry  and  metal  mining  and  pro- 
cessing are  expected  to  provide  a  smaller  proportion  of  eco- 
nomic activity  in  the  year  2000  than  now.   This  decline  is 
due  to  a  shrinking  resource  base  and  steadily  increasing 
costs  of  extraction. 


Forestry  and  kindred  industries  are  expected  to  grow 
modestly  but  steadily  through  the  end  of  this  century. 
Agriculture  will  remain  the  dominant  industry  throughout 
the  remainder  of  this  century  and  shows  significant  eco- 
nomic growth.   The  coal  industry  should  grow  dramatically 
and  will  be  the  second  largest  industry  in  Montana' by  the 
year  2000. 

This  mixed  economic  activity  will  change  the  proportion 
of  fuel  demand  in  Montana.   One  rather  dramatic  effect  of  this 
change  in  fuel  usage  is  a  much  smaller  increase  in  natural  gas 
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demand  than  has  been  experienced  in  the  past.   The  demand 
for  the  baseline  case  would  be  83  billion  cubic  feet  in 
1980  and  approximately  the  same  in  1985.   Demand  would  be 
about  89  billion  cubic  feet  in  the  year  2000,  with  peak 
demand  in  the  mid-90 's  of  around  94  billion  cubic  feet  per 
year. 

B.    Residential  Conservation  Effects 

The  baseline  case  for  gas  demand  assumes  very  little 
residential  energy  conservation.   The  first  case  to  be 
considered  after  establishment  of  the  baseline  is  the  prob- 
able effect  of  significant  home  energy  conservation.   Since 
the  greatest  volume  use  of  residential  gas  is  in  home 
heating,  the  focus  was  concentrated  in  this  area.   Reduced 
energy  use  per  household  is  due  to  increased  insulation, 
improved  weather-stripping,  and  the  use  of  storm  windows 
and  doors.   Additional  savings  are  realized  from  altered 
living  habits.   Home  energy  conservation  primarily  is  a 
function  of  gas  price. 

Figure  2  shows  the  effect  of  home  energy  conservation 
on  gas  demand  in  Montana.   The  gas  price  is  assumed  to 
rise  from  current  levels  to  $4.5  0/MCF  in  the  year  2000.   Gas 
demand  approximates  80  billion  cubic  feet  in  1980  and  declines 
to  77  billion  cubic  feet  by  1985.   Demand  for  natural  gas 
will  be  about  74  billion  cubic  feet  in  the  year  2000,  with  a 
peak  in  the  mid-1990' s  of  approximately  85  billion  cubic  feet 
per  year. 


18 


C.    Industrial  Energy  Substitution 

No  industrial  energy  substitution  is  assumed  in  the 
baseline  gas  demand  for  Montana.   The  next  step  in  the  analy- 
sis is  to  see  how  substantial  substitution  of  other  fuels 
for  gas  affects  gas  demand.   Fuel  switches  were  divided  into 
two  types.   Industrial  fuel  substitutions  considered  firm 
are  those  that  currently  are  underway  or  are  scheduled  for 
completion  by  1980.   These  switches  in  fuel  usage  are  caused 
by  various  factors  such  as  environmental  concern,  fuel  sup- 
ply uncertainties,  and  relative  fuel  costs.   The  industrial 
fuel  substitution  then  was  extrapolated  beyond  19  80.   The 
extrapolated  fuel  substitution  is  a  function  of  relative 
fuel  costs  and  the  payback  period  for  conversion. 

The  effects  of  both  types  of  industrial  fuel  substitu- 
tion are  shown  in  Figure  3.   One  curve  shows  the  effects  of 
firm  substitutions  only,  while  the  lower  curve  represents 
demand  when  firm  substitutions  and  extrapolated  substitutions 
are  combined.   The  combined  substitutions  project  a  gas  de- 
mand of  69  billion  cubic  feet  in  1980.   Demand  falls  off  to 
61  billion  cubic  feet  in  1985  and  is  about  62  billion  cubic 
feet  by  the  end  of  this  century.   Peak  usage  in  the  mid-1990' s 
is  approximately  67  billion  cubic  feet.   The  price  of  gas  is 
assumed  to  rise  from  today's  prices  to  $4.00/MCF  by  the  year 
2000. 
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D.    Gas  Price  Effects 

After  showing  the  general  effects  residential  gas  con- 
vation  and  industrial  fuel  substitution  can  have  on  gas 
demand,  the  SNG  study  group  used  several  gas  price  scenarios 
to  determine  a  range  of  possible  gas  demand  in  Montana.   Gas 
demand  curves  shown  so  far  have  assumed  a  rise  in  gas  prices 
from  the  current  level  to  $4.50/MCF  in  the  year  2000. 

Three  price  scenarios  were  selected  as  representative 
of  low,  medium,  and  high  gas  prices.   Prices  are  for  residen- 
tial gas  use;  industrial  gas  prices  were  assumed  to  be  lower 
by  $.50/MCF.   Figure  4  illustrates  the  gas  demand  expected 
at  the  specified  price  levels.   Also  shown  on  the  graph  are 
the  baseline  demand  and  a  curve  representing  the  effects  of 
firm  industrial  fuel  substitution  (no  residential  conservation 
or  extrapolated  fuel  substitution  after  1980) .   Prices  assumed 
for  the  low,  medium,  and  high  gas  prices  are  shown  in  Figure 
5.   Table  1  summarizes  the  results  shown  in  Figure  4. 


E.    Conclusions 

Natural  gas  demand  in  Montana  (including  demand  for  SNG) 
will  decline  from  the  present  level  of  about  80  billion  cubic 
feet  per  year  to  50-70  billion  cubic  feet  per  year  by  the 
year  2000.   The  decline  in  gas  demand  after  1980  is  caused 
primarily  by  projected  increases  in  the  gas  price  to  consumers 
Even  for  modest  increases  in  price  (to  $3/MCF  in  the  year 
2000) ,  the  deamnd  is  reduced  by  about  30%  from  current  demand. 
A  price  of  $3/MCF  by  the  year  2000  is  an  average  price  growth 
rate  of  only  about  2.5%  per  year.   An  average  price  growth 
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Figure  5. --Natural   Gas   Price  Scenarios 
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Table  1. --Natural  Gas  Demand  Level 

(billions  of  standard  cubic  feet) 


PRICE 


Year 

Low 
($2.00  in 

2000) 

($3 

Medium 
.00  in 

2000) 

High 
($6.50  in  2000) 

1980 

67 

67 

66 

1985 

67 

56 

54 

2000 

69 

55 

51 
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rate  of  5.75%  (price  goes  to  $6.50/MCF  in  2000)  will  drive 
gas  demand  down  to  about  60%  of  current  levels  by  the  year 
2000. 

It  should  be  noted  that  the  natural/synthetic  gas  demand 
projections  shown  here  are  probably  maximum  or  "upper  bound" 
estimates.  Inter-fuel  substitution  for  heating  in  the  resi- 
dential and  commercial  sectors  was  not  modeled,  nor  was  con- 
servation in  the  commercial  sectors.  As  the  price  of  natural/ 
synthetic  gas  increases  in  the  future,  these  factors  may  be- 
come important  and  hence  the  demand  for  natural  gas  would  be 
lower  than  that  reported  here.  Unfortunately,  the  impact  of 
these  variables  cannot  be  quantified  at  the  present  time. 

Although  we  project  a  decreasing  gas  demand,  our  "most 
likely"  projections  still  indicate  a  Montana  natural/synthetic 
gas  demand  approaching  55-60  billion  cubic  feet  per  year  in 
the  year  2000. 
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IV.   COMPETITIVE  SUBSTITUTE  AND  SUPPLEMENTAL  ENERGY  FORMS 

A  number  of  substitute  energy  forms  which  are  or  could 
be  competitive  with  natural/synthetic  gas  were  identified  in 
the  course  of  this  study.  Only  in  a  few  industrial  processes 
is  gas  actually  needed.  Natural/synthetic  gas  is  primarily 
used  as  a  source  of  heat  or  as  a  boiler  fuel.  Other  energy 
forms  can  certainly  be  substituted  for  natural  gas  although 
they  may  not  be  as  convenient  or  clean  burning. 

Both  supplemental  natural/synthetic  gas  supplies  which 
could  be  utilized  in  the  current  natural  gas  delivery  system 
and  fuels  requiring  new  energy  delivery  systems  were  identified 
and  evaluated.   Our  inter-fuel  substitution  research  focused 
primarily  on  the  industrial  sector  in  Montana  since  significant 
fuel  substitution  is  currently  ongoing  in  that  sector  and  more 
is  expected.   Industries  also  require  quantities  of  energy 
sufficient  to  make  the  utilization  of  new  fuel  delivery 
systems  particularly  attractive. 

It  was  realized  early  that  coal  could  be  synthesized  into 
substitute  natural  gas  which  could  be  delivered  through  existing 
pipelines  and  combusted  in  existing  furnaces  and  burner  heads 
without  modification.   Detailed  technical  and  economic  evalu- 
ations of  the  various  processes  for  producing  synthetic 
natural  gas  are  being  evaluated  elsewhere. 

Several  additional  supplemental  sources  of  natural  gas 
were  also  identified.   As  the  price  of  natural  gas  increases, 
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Montana  natural  gas  producers  expect  to  be  able  to  discover 
and  produce  more  natural  gas  than  in  the  recent  historic  past 
when  natural  gas  prices  were  constrained  to  artificially  low 
levels.   Similarly,  natural  gas  exploration  is  expected  to 
increase  in  Montana's  neighboring  states.   Montana  could  con- 
ceivably receive  additional  supplies  of  natural  gas  from  these 
sources;  however,  neither  the  quantity  nor  the  price  of  these 
additional  supplies  can  be  determined  accurately  at  the  present 
time. 

A  major  reason  for  Montana's  potential  natural  gas  supply 
shortfall  is  a  projected  curtailment  of  natural  gas  supplies 
which  are  currently  imported  from  Canada.   Montana  Power  has 
the  right  to  import  34.2  billion  cubic  feet  (BCF)  of  natural 
gas  from  Canada  during  the  year  which  will  end  May  13,  1977. 
They  then  expect  to  import  29.2  BCF  annually  until  1985.   Under 
current  circumstances  Canadian  gas  will  be  completely  eliminated 
by  about  1989.    It  is  impossible  to  assess  what  the  Canadian 
government  may  or  may  not  do  in  the  future.   Large  reserves 
of  natural  gas  may  exist  in  the  Canadian  Arctic  which  could 
be  developed  in  the  future.   Based  on  current  policies,  however, 
it  seems  unreasonable  to  expect  that  additional  supplementary 
gas  supplies  from  Canada  will  become  available  to  Montana. 


Large  reserves  of  natural  gas  are  also  known  or  believed 
to  exist  in  the  state  of  Alaska.   These  reserves  are  expected 
to  be  developed  in  the  near  future  with  the  first  gas  supplies 
from  this  source  possibly  reaching  the  continental  United 
States  by  1982.   Several  routes  have  been  suggested  to  bring 
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these  supplies  to  the  Continental  United  States  and  the  matter 
is  currently  before  the  Federal  Power  Commission.   Routes  under 
consideration  include  pipelines  through  Canada  along  the  Alcan 
Highway  or  down  the  MacKenzie  Valley  or  a  pipeline  entirely 
through  Alaska  with  Liquified  Natural  Gas  (LNG)  transport  from 
Alaska  to  the  Continental  United  States.   Montana  Power  cur- 

5 

rently  has  an  option  to  purchase  15  BCF  of  this  gas  annually. 
However,  the  price  of  this  gas  may  be  $3/MCF  or  higher.   The 
United  States  could  also  import  LNG  from  foreign  producers; 
but,  this  gas  is  also  likely  to  cost  over  $3/MCF  delivered 
and  result  in  problems  of  supply  security. 

Several  substitute  energy  forms  which  would  require  new 
or  modified  energy  delivery  systems  were  also  identified. 
After  preliminary  review  it  was  determined  that  low  Btu  coal 
gasification,  direct  coal  combustion,  and  wood  or  wood  waste 
combustion  are  possible  economic  substitutes  for  natural  gas 
in  the  Montana  industrial  sector. 


The  wood  products  industry  in  Montana  currently  utilizes 
a  large  quantity  of  wood  materials  for  energy  generation.   It 
is  certainly  possible  that  more  could  be  used  in  the  future. 
Montana  currently  produces  over  8  million  tons  of  wood  log- 
ging residue  per  year  which  is  either  burned  at  the  logging 

6 

site  or  simply  left  piled  in  place.    Unfortunately,  the  long 
term  delivered  cost  of  this  fuel  is  uncertain;  and,  the  de- 
tailed assessment  required  to  fully  and  accurately  evaluate 
the  potential  for  industrial  utilization  of  this  energy  source 

7 

was  beyond  the  scope  of  this  study. 
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The  direct  combustion  of  coal  as  a  substitute  for  natural 
gas  is  an  alternative  that  several  industries  in  Montana  have 
recently  turned  to.   In  April  1976  the  Ideal  Cement  Plant  at 
Trident,  Montana  converted  from  natural  gas  to  coal  at  a 
capital  cost  of  $2.3  million,  or  approximately  $1.20  per 

8 

million  Btu.    The  Kaiser  Cement  plant  near  Helena,  Montana, 
has  a  similar  conversion  to  coal  planned  for  the  fall  of  1977 

9 

at  a  capital  cost  of  $1.5  million  (%$ . 90/million  Btu).   Finally, 
the  Great  Western  Sugar  Company  in  Billings,  Montana  has 
announced  plans  to  convert  their  sugar  beet  processing  plant 
to  coal.   The  capital  cost  of  this  conversion  is  reported  to 

10 

be  $3.21  million  (^-$2  .  70/million  Btu)  .     Air  pollution  con- 
trol costs  for  this  facility  are  over  $1.2  million.   As  can 
be  seen  above  the  capital  costs  of  conversion  to  coal  from 
natural  gas  are  highly  variable  depending  upon  the  particular 
industrial  facility  considered.   Normalized  conversion  costs 
also  vary  inversely  with  the  size  of  the  facility  converted. 
Fuel  costs  for  coal  delivered  in  Montana  vary  with  location. 
Reasonable  fuel  cost  (1975  dollars)  estimates  for  Colstrip, 
Montana  coal  delivered  to  various  Montana  cities  are:   Butte 

$.70/million  Btu;  Missoula  $.80/million  Btu;  Billings  $.50/ 

i  i 
million  Btu.     These  coal  fuel  costs  are  certainly  favorable 

when  compared  to  natural/synthetic  gas.   Coal  fuel  costs  do 

not  seem  likely  to  escalate  rapidly  in  the  future  either. 


Coal  fired  boilers  are  currently  available  utilizing 
either  spreader  grate  or  pulverized  coal  firing.   Spreader 

c 
grate  boilers  can  be  employed  up  to  abut  0.5  x  10  lb  of  steam 
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per  hour;  however,  the  sized  coal  required  for  spreader  grate 

12 
operation  may  not  be  readily  available.     Both  boiler  design 

and  cost  are  particularly  sensitive  to  the  specific  type  of 
coal  to  be  used.   The  Oak  Ridge  National  Laboratory  estimates 
that  the  annual  cost  of  coal  fired  steam  generation  using 
western  coal  (8500  Btu/lb,  0.5%S)  to  be  $1 . 53/million  Btu. 
They  also  estimate  that  flue  gas  cleaning  processes  are  avail- 
able for  ^$ . 50/million  Btu  additional.   The  Oak  Ridge  esti- 
mate is  split  almost  equally  between  capital  and  fuel  costs. 

The  favorable  economics  of  direct  coal  combustion  is 
somewhat  contrasted  by  its  unfavorable  environmental  impact. 
Even  with  the  "best  available"  air  pollution  control  equipment, 
a  coal  fired  boiler  will  emitt  more  pollutants  than  a  conven- 
tional gas  fired  boiler.   As  previously  noted  a  significant 
portion  of  the  cost  of  coal  fired  boilers  will  be  allocated 
to  air  pollution  control.   Although  direct  coal  fired  combus- 
tion (with  proper  pollution  control)  appears  suitable  for 
areas  with  adequate  pollution  dispersion  potential,  there  are 
certain  areas  within  the  state  (particularly  in  western  Montana) 
where  conventional  coal  combustion  may  simply  be  environmentally 
unacceptable. 


Fluidized  bed  combustion,  a  new  alternative  process  for 
direct  coal  conversion  is  currently  under  development  and 
holds  high  promise  for  the  future.   Fluidized  bed  combustion 
is  undertaken  in  an  inert  bed  of  coal  ash  and  crushed  lime- 
stone which  rests  on  a  bed  of  nozzles.   Combustion  air  enters 
the  bed  through  the  nozzles,  expanding  the  bed  and  causing  it 
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to  flow  much  like  a  liquid.   Combustion  is  accomplished  by 
introducing  crushed  coal  and  raising  the  bed  to  the  ignition 
temperature  of  the  coal.   Volumetric  heat  release  rates  for 
fluidized  bed  combustion  are  on  the  order  of  10  times  those 
of  conventional  pulverized  coal  furnaces.   Heat  transfer  sur- 
faces are  placed  within  the  bed  itself  to  absorb  at  least 
half  of  the  heat  released;  the  remainder  of  the  heat  released 
is  collected  by  conversion  surfaces.   The  combustion  temperature 
of  the  bed  is  controlled  to  1600-1800°F ,  considerably  lower 
than  conventional  coal  combustion.   Since  this  combustion  tem- 
perature is  below  the  ash  fusion  temperature  of  coal,  fluidized 
bed  boilers  can  be  designed  independent  of  the  type  of  coal 
to  be  used. 

One  of  the  major  attractive  features  of  fluidized  bed 
combustion  is  its  potential  for  air  pollution  emission  control. 
Emission  control  is  inherent  in  the  fluidized  bed  combustion 
process.   If  limestone  is  injected  into  the  fluidized  bed, 
the  bed  turbulence  will  promote  the  reaction  of  sulfur  dioxide 
with  limestone  and  effectively  remove  sulfur  from  the  system 
as  dry  calcium  sulfate.   Studies  on  the  potential  for  lime- 
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stone  regeneration  are  currently  underway.     The  relatively 
low  combustion  temperature  is  the  fluidized  bed  combustion 
temperature  in  the  fluidized  bed  combustion  process  also  limits 
the  formation  of  oxides  of  nitrogen  as  does  the  minimum  of 
excess  air  required  for  combustion. 


There  are  no  fluidized  bed  boilers  commercially  available 
in  the  United  States  today.   A  30  Megawatt  pilot  fluidized  bed 
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boiler  was  recently  built  in  Rivesville,  West  Virginia  and  is 
currently  beginning  to  undergo  testing.     The  Foster-Wheeler 
Company,  manufacturer  of  some  of  the  components  of  the  Rivesville 
plant,  indicate  that  they  will  warranty  units  of  comparable  size 
after  6  months  of  successful  operation.   The  city  of  Linkoping, 
Sweden  is  currently  evaluating  bids  from  several  vendors  to 
build  a  15  MW  boiler  for  its  district  heating  system;  hence,  at- 
least  intermediate  sized  fluidized  bed  boilers  must  be  available 
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in  Europe.     Cost  data  for  fluidized  bed  boilers  are  not  yet 
fully  developed  in  the  United  States.   In  1970,  Pope,  Evans, 
and  Rogers  estimated  that  a  fluidized  bed  boiler  would  cost 
approximately  1.4  times  that  of  an  equivalent  gas  fired  boiler 
while  a  conventional  spreader  grate  coal  fired  boiler  would 
cost  2.35  times  more  than  a  gas  fired  boiler.     Oak  Ridge 
National  Laboratory  predicts  slightly  lower  steam  costs  ($1.65/ 
million  Btu  versus  $1 . 84/million  Btu)  for  a  fluidized  bed  boiler 
compared  to  a  conventional  boiler  with  stack  gas  desulf urization 
at  an  equivalent  location.   The  lower  costs  for  fluidized  bed 
combustion  are  projected  because  of  their  more  compact,  stan- 
dardized design,  their  reduced  surface  area,  and  the  fact  that 
they  can  be  factory  rather  than  field  assembled.   Applications 
of  fluidized  bed  combustion  for  residential  and  commercial 
use  are  also  being  researched.   The  Montana  Energy  and  MHD 
Research  and  Development  Institute  is  particularly  interested 
in  a  novel  fluidized  bed  furnace  that  is  applicable  to  home 
heating  uses. 
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Fluidized  bed  coal  combustion  appears  to  be  a  particularly 
useful  alternative  to  natural  gas  in  Montana.  While  the  process 
is  not  widely  available  now  it  seems  to  be  ready  for  commercial 
development  and  must  be  considered  when  evaluating  future  alter- 
natives. A  detailed  investigation  of  its  suitability  and  poten- 
tial a:  an  alternative  energy  form  appears  warranted  on  the 
basis  of  this  investigation. 

Finally,  the  potential  applicability  of  low  Btu  coal 
gasification  in  the  Montana  industrial  sector  was  evaluated. 
Low  Btu  coal  gasification  is  essentially  the  first  step  in  the 
process  of  producing  synthetic  pipeline  quality  natural  gas. 
Coal  is  reacted  with  steam  at  high  temperatures  and  pressures 
to  produce  carbon  monoxide  and  hydrogen.   If  the  reaction  oc- 
currs  in  the  presence  of  air  the  resultant  gas  is  diluted  with 
nitorgen  and  has  a  heating  value  of  approximately  125-150  Btu/ 
cubic  feet  (SCF) .   If  the  initial  reaction  occurrs  in  an  oxygen 
atmosphere  the  heating  value  of  the  resulting  gas  is  increased 
to  nearly  300  Btu/SCF.   Conversion  from  natural  gas  to  low  Btu 
gas  may  require  boiler  derating  and/or  burner  head  modification 
depending  on  the  exact  heating  value  of  the  gas  produced. 


The  technology  of  low  Btu  coal  gasification  is  well  de- 
veloped.  Before  natural  gas  was  widely  available,  low  Btu 
gas  (town  gas)  was  used  for  home  lighting  and  heating  and  in 
industrial  processes  throughout  Europe  and  the  eastern  United 
States.   There  is  a  single  brick  manufacturer  in  Pennsylvania 
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that  currently  gasifies  about  85  tons  of  coal  per  day.     World- 
wide there  are  many  low  Btu  gasifiers  operating.   There  are 
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several  manufacturers  that  currently  offer  low  Btu  gasification 
reactors . 

Low  Btu  coal  gasification  can  be  considerably  cleaner 
than  conventional  direct  coal  combustion.   Particulate  air 
pollution  control  is  to  a  large  degree  inherent  in  the  gasi- 
fication process  since  primarily  gaseous  products  leave  the 
reaction  vessel.   The  gaseous  stream  is  usually  further  scrubbed 
to  remove  most  of  the  remaining  entrained  particulate  matter. 
Sulfur  also  can  be  removed  effectively  since  it  will  be  in  the 
gas  stream  as  hydrogen  sulfide.   There  are  a  number  of  commer- 
cial processes  available  to  remove  hydrogen  sulfide;  however, 
because  sulfur  removal  is  not  inherent  in  the  process,  addi- 
tional costs  are  incurred  for  this  sulfur  removal.   Some  manu- 
facturers claim  that  by  gasifying  low-sulfur  coal  they  can 
meet  applicable  federal  air  pollution  regulations  without 
sulfur  removal.   This  tactic  is  probably  inappropriate  for 
Montana  industry  since  many  of  the  industries  which  could 
utilize  low  Btu  gasification  are  in  areas  where  ambient  air 
quality  standards  are  already  approached  or  exceeded. 


Several  cost  estimates  for  low  Btu  gasification  have 
recently  been  developed.   Low  Btu  gasification  should  inher- 
ently be  cheaper  than  producing  pipeline  quality  synthetic 
natual  gas  since  neither  the  water  gas  shift  reaction  for  rela- 
tive hydrogen  enrichment  or  methanation  are  required.   Low  Btu 
gasification  facilities  are  also  generally  located  at  the  load 
center  so  that  gas  transportation  costs  are  minimized.   Unfor- 
tunately, most  recent  reviews  indicate  that  low  Btu  gas  will 
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cost  in  the  range  of  $3 . 00-$4 . 00/million  Btus.      Costs  in- 
crease as  the  size  of  the  facility  decreases.   However,  a 
number  of  sources  still  do  predict  that  low  Btu  gasification 
will  be  less  expensive  than  pipeline  quality  synthetic  gas  pro- 
duction. 

In  summary  one  can  identify  a  number  of  supplemental  or 
substitute  energy  forms  which  are  applicable  to  Montana.   These 
supplements  or  substitutes  will  all  be  much  more  expensive  than 
natural  gas  historically  has  been.   Some  of  these  substitutes/ 
supplements  will  undoubtedly  be  utilized  in  the  future;  however, 
their  high  costs  signal  that  conservation  will  also  have  a 
significant  impact  on  the  future  demand  for  natural/synthetic 
gas  in  Montana. 


34 


ACKNOWLEDGEMENTS 

Our  special  thanks  go  out  to  the  many  individuals  and 
organizations  who  have  aided  this  project.   Particular  appre- 
ciation is  extended  to  the  people  of  the  Montana  Department 
of  Community  Affairs,  the  Montana  Energy  Advisory  Council, 
the  Montana  Power  Company,  and  the  V7illamette  Simulation 
Unit. 


35 


REFERENCES 


1.  J.  W.  Forrester,  Industrial  Dynamics,  The  MIT  Press, 
Cambridge,  Massachusetts,  1961. 

2.  Oregon  State  Simulation  Model,  Willamette  Simulation  Unit, 
Oregon  State  University. 

3.  Montana  Futures  Project,  Division  of  Research  and  Infor- 
mation Systems,  Department  of  Community  Affairs,  state 
of  Montana. 

4.  Johnson,  M.C.,  "Montana's  Economy:   Where  Its  Been  and 
Where  Its  Going,"   Montana  Business  Quarterly,  Volume 
14,  No.  1,  Winter  19~76~7 

5.  Frizzell,  Thomas  W. ,  Montana  Environmental  Quality  Council, 
personal  communication.   August  1976. 

6.  Coldiron,  W.  H.,  Letter  to  Montana's  Lieutenant  Governor, 
Bill  Christiansen,  Helena,  Montana.   May  1976. 

7.  Billings,  Roger,  Montana  Power  Company,  Personal  Com- 
munication, October  1976. 

8.  Data  derived  from  a  MERDI  research  proposal. 

9.  Quinney,  Dean  N.  "Economics  of  Utilizing  Residuals  from 
Logging  -  Problems  and  Opportunities,"  in  Forest  Product 
Residuals,  W.  K.  Lautner,  Editor.   American  Institute  of 
Chemical  Engineers,  Symposium  Series  146,  Volume  71, 
1975. 

10.  Montana  Standard,  April  9,  1976,  p.  9. 

11.  Frizzell,  Thomas  W. ,  Montana  Environmental  Quality  Council, 
personal  communication. 

12.  Montana  State  Department  of  Health  and  Environmental 
Sciences,  "Preliminary  Environmental  Review  for  Boiler 
Conversions  at  Great  Western  Sugar  Company,  Billings, 
Montana,"  June  1976. 

13.  Adapted  from  Berg,  Lloyd,  "Conversion  Cost  Estimates  for 
the  Montana  University  System",  Rough  Draft,  Bozeman, 
Montana,  September  1975. 

14.  Anderson,  T.  D.,  et.  al .  "An  Assessment  of  Industrial 
Energy  Options  based  on  Coal  and  Nuclear  Systems."   Oak 
Ridge,  Tennessee,  Oak  Ridge  National  Laboratory,  July 
1975. 


36 


15.  Ibid. 

16.  Hammond,  Allen  L.,  "Coal  Research  (IV):   Direct  Combustion 
Lags  its  Potential",  Science,  Volume  194,  pp.  172-194, 
October  8,  1976. 

17.  Lovins,  A.  B.,  "Energy  Strategy:   The  Road  Not  Taken?" 
Foreign  Affairs,  pp.  65-95,  October  1976. 

18.  Pope,  Evans,  and  Rogers,  "The  Development  of  Coal  Fired 
Fluidized  Bed  Boilers,"  Final  Report,  Volume  1,  OCR  R  &    D, 
Report  No.  36,  February  1970. 

19.  Tillman,  David  A.,  "Status  of  Coal  Gasification,"  Environ- 
mental Science  and  Technology,  Volume  10,  No.  1,  pp.  34-38, 
January  19  76. 

20.  Stokes,  Charles  A.,  "Gas  from  Coal  has  a  Future  as  a 
Boiler  Fuel,"  paper  presented  at  the  American  Chemical 
Society  Annual  Meeting,  San  Francisco,  California,  August 
1976. 

21.  Ferretti,  Emmett  J.,  "Technology  and  Economics  of  Indus- 
trial Fuel  Gas  from  Coal,"  paper  presented  at  the  American 
Chemical  Society  Annual  Meeting,  New  York,  New  York, 
April  1976. 


37 


APPENDIX  A 

Narrative  Description  of  the  Montana 
Dynamic  System  Simulation  Model 

This  appendix  is  divided  into  sections  as  follows: 

Page 

Section  A-I,   Overview  of  the  Model A-l 

Section  A-I  I,   Demographic  Sub-Model A-5 

Section  A-I  II,   Land  Use  Sub-Model A-7 

Section  A-IV,   Economic  Sub-Models A-l  1 

Section  A-V,   Economic   Input/Output  Sub-Model    A-45 

Section  A- VI,   Capital    Investment  Sub-Model    A-48 

Section  A- VI I ,   Employment  Sub-Model A-52 

Section  A-VIII,   Energy  Sub-Model    A- 59 

Section  A-IX,   Residential   Energy  Conservation  Sub-Model A-83 


Section  A- I,  Overview  of  the  Model 

As  in  many  other  regions  of  the  country,  the  supply  of  natural  gas  in 
Montana  may  not  be  sufficient  to  meet  future  needs.  One  method  to  avoid 
this  problem  in  Montana  is  to  produce  SNG  (synthetic  natural  gas)  using  the 
extensive  coal  reserves.  The  purpose  of  this  study  is  to  determine  the  de- 
mand for  SNG  in  Montana. 

Determining  the  demand  for  SNG  is  essentially  a  balancing  problem.  The 
SNG  demand  is  the  additional  production  required  to  balance  natural  gas  re- 
quirements with  the  supplies  available.  Thus,  the  demand  for  and  supply  of 
natural  gas  must  be  determined  before  the  demand  for  SNG  can  be  found. 
Natural  gas  cannot  be  considered  alone  as  it  competes  with  other  fuels. 
Major  conversion  of  industry  gas  uses  to  coal  and  electricity  are  certainly 
possible  and  can  significantly  affect  the  demand  for  natural  gas.  The  cost 
of  producing  SNG  is  likely  to  be  different  from  the  natural  gas  obtained 
from  conventional  sources.  Thus,  the  part  of  the  demand  filled  by  SNG  may 
be  complete  at  different  price  levels  than  other  natural  gas. 

The  approach  used  herewasto  develop  a  model  of  Montana's  economy  on 
which  energy  demands  can  be  based.  The  overall  structure  is  shown  in  Figure 
1.  The  industrial  and  commercial  demands  for  energy  are  based  upon  the  levels 
of  activity  in  the  corresponding  sectors  of  the  economy  while  residential 
demands  are  based  upon  the  number  of  households  and  their  income.   In  addition, 
the  effects  of  inter-fuel  competition  and  conservation  are  considered. 

The  model  is  similar  in  concept  to  the  Oregon  State  Simulation  Model  (1) 
and  the  structure  of  some  of  the  sub-models  is  basically  the  same  as  in  their 
Oregon  counterparts.  The  Montana  Futures  Project  input-output  model  (2)  is 
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also  used  in  the  formulation  of  the  model.  The  primary  modeling  technique 
used  is  dynamic  systems  simulation. (3>  This  methodology  deals  with  the 
forces  causing  change  in  the  economy  and  should  be  well  suited  for  the 
transients  which  are  likely  in  Montana's  future.   It  also  allows  the  model 
to  be  based  more  closely  on  the  actual  processes  and  actions  which  occur  in 
Montana's  economy.  The  dynamic  simulation  is  combined  with  input/output 
analysis  in  determining  levels  of  economic  activity.  The  use  of  the  input/ 
output  analysis  assists  in  obtaining  more  detail  in  this  critical  area. 

Each  sub-model  is  discussed  separately  in  the  following  sections.  An 
effort  has  been  made  throughout  to  design  the  models  such  that  they  can 
easily  be  improved  as  additional  information  becomes  available.   In  this 
manner,  major  model  changes  will  not  be  required  when  better  data  are 
obtained. 
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Section  A- II,  Demographic  Sub-Model 

The  demographic  sub-model  predicts  the  number  of  people  that  will  be 
in  Montana  and  the  number  of  households  in  the  state.  The  projections  are 
based  upon  current  population  and  predicted  net  migration.  The  migration 
is  affected  by  the  difference  between  Montana  unemployment  levels  and 
national  unemployment  levels. 


A- 5 
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Section  A- III,  Land  Use  Sub-Model 


The  land  use  component  maintains  a  current  inventory  of  land  in  each  of 
six  categories.  The  following  description  details  the  flow  between  cate- 
gories. 

Residential  and  commercial-industrial  land  is  directly  converted  from 
agricultural  land.  The  demand  for  residential  land  is  quite  different  from 
the  demand  for  commercial-industrial  land  due  to  the  relative  intensities 
of  use.  Demand  for  residential  land  is  driven  by  the  demographic  component. 
This  rate  is  a  function  of  the  rate  of  household  formation  and  the  residen- 
tial conversion  density.   In  mathematical  form  the  change  in  demand  for 
residential  land  is  given  by: 

ADEMR  =  AHH  i   CD 
where  DEMR  is  demand  for  residential  land,  HH  is  the  number  of  households 

(from  demographic  component),  and  CD  is  the  conversion  density  (total  households 

total  residential  acres). 

The  commercial-industrial  (C-I)  land  use  demand   is  generated  in  the 

levels  of  activity  in  the  industrial  sectors  of  the  economic  component.  The 

rate  of  demand  for  the  C-I  land  is  determined  by  the  rate  of  employment  and 

the  employees/acre  in  the  C-I  sectors.   In  mathematical  form: 

ADEMp  =  AEMC  t  EA 

where  DEMq  is  demand  for  C-I  land  (acres),  EM^  is  employment  in  C-I  sectors, 

and  EA  is  total  employees  in  C-I  sector 
total  acres  in  C-I  land. 

The  demand  for  residential  and  C-I  land  is  summed  to  obtain  the  total 
land  converted  from  agricultural  to  urban  land.   In  mathematical  form: 

ADEMy  =  ADEMC  +  ADEMR 
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where  DEMj  is  total  urban  land  demanded. 

The  mining  land  is  directly  converted  from  agricultural  land,  since  the 
stripped  land  in  Montana  is  almost  all  used  for  crops  or  grazing.  The  con- 
version rate  is  a  scenario  variable  derived  from  standard  coal  production 
per  year.  The  land  mined  per  year  is  calculated  by  dividing  the  tons  of 
coal  produced  by  the  tons  of  coal  per  acre. 

Mined  land  converts  to  waste  land  or  reclamation  land.  The  stripped 
land  which  falls  into  the  reclamation  land  category  is  revegetated  over 
a  period  of  years.  The  time  span  is  a  technical  policy  input.  Waste 
land  is  derived  from  the  product  of  land  mined  and  fraction  of  revege- 
tation,  this  fraction  being  a  technical  policy  input. 

Forest  land  use  accounts  for  all  the  land  upon  which  commercially 
harvestable  timber  is  growing.  Agricultural  land  can  be  converted  to 
forest  land;  however,  in  Montana  this  is  unlikely.  Policy  input  for 
reforestation  is  0. 

Wilderness  land  is  an  area  which  can  be  converted  from  forest  land 
to  preservation  status  by  a  policy  input.  This  relationship  can  also  be 
reversed,  wilderness  land  converted  to  forest  land. 
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LAND  USE  DATA 


Land  Use  (1 ,000  Acres) 

1967 

Agricultural  Land 

(5) 

62,90C 

Forest  Land  Total 
Wilderness  (3) 
Logging 

(2) 

18,805 

2,145 

16,660 

Urban  (2) 

818 

Mining* 

4 

Other  (waste)  ** 

10,562 

TOTAL  LAND  AREA 

93,089 

Coal  Mined  =  11  x  106  tons/yr  (1973)  (4) 

Tons/Acre  =  28320  T/A.  *** 

Residential  Land  -613,500  (assumed  75%  of  Urban)  (3) 

Industrial  Land  =204,500  (assumed  25%  of  Urban)  (3) 

Households/Acre  =.354  (4) 

Employees/Acre  =  .2396  (4) 

Total  Households  =  217,300  (4) 

Total  Employees  (non  farm)  =  49,000  (4) 


*   Mining  Land  =  (acres  mined  in  last  10  vears) 

**   Other  land  =  Total  land  -  (Agriculture  Land  +  Forest  Land  +  Urban  Land  + 
Mining  Land) . 

***  Tons/Acre  =  11  x  106  tons/yr  / 1416  tons/acre  ft  (.8  recovery7)  (25  ft 
average  coal  seam8) . 
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Section  IV,  Economic  Sub-Models 

The  economic  sub-model  is  structured  as  shown  in  Figure  1.  Five 
separate  sectors  are  used  for  the  primary  activities:  agriculture,  forestry: 
coal  mining,  oil  and  gas  production,  and  non-energy  mining.  Since  the  out- 
put of  the  petroleum  refining,  and  primary  metal  processes  depends  mainly 
on  the  amount  of  oil  produced  and  metal  mined  in  Montana,  they  are  also 
considered  separately.  This  part  of  the  model  can  be  considered  supply 
limited.  The  output  of  these  sectors  plus  consumer  expenditures,  export 
demands,  government  spending,  and  capital  expenditures,  drive  an  input- 
output  model  for  the  rest  of  the  economy.  These  activities  can  be  con- 
sidered demand  limited.  The  outputs  are  then  used  to  determine  total 
employment  and  income. 

Each  of  the  sectors  for  the  primary  activities  and  for  petroleum  re- 
fining and  primary  metal  processing  are  modeled  separately.  The 
capital  investment  and  employment  models  are  also  discussed. 


A-ll 


ECONOMIC    MODEL 


CRUDE   OIL 

IMPORTS 

AVAILABLE 


GOVERNMENT 

SPENDING  EXPORT 

DEMAND 


Figure  1. --Overall  Structure  of  the  Economic  Sub-Model 
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Agricultural   Sector 

The  agricultural   sub-component  of  the  economic  model    is  divided  into 
two  parts,  agricultural    production  and  capital    formation.1      Figure  1    illu- 
strates  the  interrelationship  between   the   two  parts  and  how  each   interfaces 
with  the  rest  of  the  model. 

Agricultural    Production:     The  amount  of  agricultural    land   and   the 

level   of  capital    employed  are  combined,   utilizing  a  modified  Cobb-Douglas 

2 
production   function,   to  generate  the  production  level. 

ALP  (1-ALP) 

Y   =  A  x  TK         x  AL 

where  Y  =  Agricultural  production;  total  annual  market  value  of  farm  products 
sold  (in  constant  dollars). 

TK  =  Agricultural  capital;  total  value  of  machinery  and  equipment  em- 
ployed as  a  factor  of  production. 

AL  =  Agricultural  land. 

ALP  =  Parameter:  relative  importance  of  capital  as  compared  to  land. 

A  =  Constant  relating  units  of  capital  and  land  to  units  of  pro- 
duction. 

Gross  production  in  dollars  is  a  function  of  Ag  output  in  constant  dollars 

and  relative  price  (P)  of  Ag  products. 

D  =  YP 

The  remainder  of  the  production  model  converts  gross  agricultural  pro- 
duction to  net  farm  recipts  by  subtracting  the  aggregate  expenses  incurred 
in  production.  The  production  expenses,  all  measured  in  1970  dollars/yr, 
fall  into  four  categories:  direct  operating  expenses,  labor  expense,  inter- 
est, expense,  and  property  tax  expense.  Direct  operating  expenses  (DOE)  in- 
clude the  cost  of  livestock  and  poultry,  feed  and  seed,  fertilizer  and  other 
chemicals,  petroleum  fuels  and  electricity,  machine  hire,  irrigation  and  mis- 
cellaneous farm  services.  Direct  operating  expenses  are  assumed  to  be  a 
fraction  of  the  gross  production: 
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DOE  =  (Y)  x  OP 
where  OP  =  operating  expense  fractional  multiplier. 

The  labor  expense  (ALE)  is  the  cost  of  hired  labor.   It  is  calculated  by 
multiplying  the  labor  demand  by  the  average  annual  v/age  rate: 

ALE  =  WR  x  (Y  v  OPL) 
where  OPL  is  the  output/labor  ratio  ($/man  year). 

The  interest  expense  (EI)  is  the  sum  of  interest  on  the  aggregate  farm 
mortgage  debt  and  interest  paid  on  financial  capital: 

EI  =  (R  x  AKF)  +  (PLF  x  AL  x  PL  x  R) 
where 

R  =  Annualized  interest  rate. 
AKF  =  Level  of  financed  capital  stock  ($). 
PLF  =  %  of  agricultural  lend  mortgaged. 
PL  =  Average  value  of  agricultural  land  ($/acre). 


The  property  tax  expense  (ET)  is  the  product  of  the  average  farm  pro- 
perty tax  rate  (FTR)  and  the  average  value  of  agricultural  land: 

ET  =  FTR  x  AL  x  PL 

The  Net  Farm  Receipts  (ANFR)  are  then  determined  by  reducing  the  value 
of  gross  agricultural  output  by  the  sum  of  the  aggregate  expenses  incurred: 

ANFR  =  PY  -  (DOE  +  EI  +  ET  +  ALE) 

Capital  Formation:  The  value  of  machinery  and  equipment  employed  in 
Montana's  agricultural  sector  has  greatly  increased  in  the  past  two  decades. 
Unfortunately,  due  to  exogenously  imposed  variations  that  are  beyond  their 
control,  farmers  are  generally  forced  to  allocate  their  resources  under 
conditions  of  imperfect  knowledge  of  the  future.   Considerable  uncertainity 
exists  concerning  the  success  of  future  crops  and  the  market  prices  of 
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agricultural  commodities.  There  seems  to  be  a  tendency  for  the  agricultural 
sector,  as  a  whole  to  over-capitalize  in  an  attempt  to  maximize  production, 
with  the  attendant  increased  costs  of  production  offsetting  the  increased 
gross  farm  sales.  Hence,  net  farm  receipts  have  failed  to  keep  pace  with 
the  steady  rise  in  productivity. 

The  capital  formation  model  starts  with  the  value  of  net  farm  receipts 
developed  in  agricultural  production  and  determines  the  resulting  change  in 
capital.  If  net  farm  receipts  are  greater  than  zero,  a  portion  is  used  to 
purchase  new  capital.  Depreciation  of  capital  stock  is  assumed  to  occur  at 
a  rate  that  is  proportional  to  the  installed  level.  The  change  in  the 
level  of  owned  capital  stock  (AKO)  is  given  by:^ 

AAKO  =      [(ANFR  x  PTI)  -  (AKO  x  AKDEP  +  DI)]  x  DT 
where  PTI  =  propensity  to  invest;  scenario  input  with  default  assumption 
that  investment  is  a  linear  function  of  net  farm  receipts. 
AKDEP  =  proportional  rate  of  depreciation  scenario  variable. 
DI  =  Disinvestment. 


It  is  also  assumed  that  farmers  in  the  aggregate  tend  to  increase  the 
sector's  level  of  capital  stock  by  borrowing.  Conventional  financing  (cash 
down  payment,  with  the  purchased  capital  acting  as  collateral)  allows  the 
sector  to  "leverage"  their  available  cash.   In  addition,  there  is  reason 
to  believe  that  the  farm  community  has  utilized  the  deferred  annual  capital 
gain  on  agricultural  land  to  help  finance  capital  investment.  The  change 
in  the  level  of  financed  capital  stock  (AKF)  is  given  by: 

AAKF  =      [(ANFR  x  PTI  x  ALEV1 )  -  (AKF  x  AKDEP)]  x  DT 
where  ALEV1  =  leverage  on  cash,  policy  input  that  accounts  for  the  addition 
to  capital  made  possible  by  using  capital  as  collateral. 
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In  the  event  net  farm  receipts  are  less  than  zero,  disinvestment  can 
occur.   If  the  magnitude  of  the  loss  exceeds  a  specific  fraction  of  gross 
production  (TF),  it  is  assumed  that  the  agricultural  community  will  be 
forced  to  further  reduce  the  level  of  owned  capitol  stock  by  liquidating 

it  at  its  salvage  value;  otherwise,  capital  is  simply  allowed  to  depreciate 

2 
without  replacement.   If  NFR(t)  >TF(Y),  then 

AAKO  -  [NFR  +  TF(Y)]  /salv  -  AKO(t-l)  AKDEP 
AAKF  =  AKF(t-l)  x  AKDEP 


However,  if  NFR    <  TF  x  Y 


,  then 


AAKO  =  -AKO(t-l)  x  AKDEP 
AAKF  =  -AKF(t-l)  x  AKDEP 
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Notes: 


Reference 


Direct  operating  expense/unit  output 
435  x  106/644  x   10&  -  0.675 

Hired  labor/unit  output 

11000  man  yrs./$644  x  105  =   17.1   man  yrs./106$ 

Wage  rates   =   $3,363/yr 

Property  tax  ra  :e 

Total    tax/(total    land  x  value) 

35.7  x  106/(62.9  x  106  x  $60)   -   9.46  x   10"3 

Value  of  land 

$  60/acre   (1969) 
$150/acre   (1976  estimate) 

Interest  rate 
0.06 

Fraction   loss  allowed 
0.20 

Salvage  value 
0.30 

Propensity  to  invest 
0.20 

Leverage  on  cash 
1.00 

Total   capital 

$420  x  10^  (non  land) 

$210  x  10°  (owned  capital) 

$210  x  106  (financed  by  leverage  on  cash) 

Fraction  of  Land  Mortgaged 
.06 

Depreciation 

$97.2  x  106/yr 


(4) 

(3) 
(*) 

(3) 

(4) 

(Assumed) 

(Assumed) 

(Assumed) 

(1) 

(1) 

(4) 


(**) 


(3) 


*     Wage  rate  was  obtained  by  dividing   total    expense  for  hired  farm  labor4 
by  total    labor3. 

Fraction  of  farm  land  mortgaged  was  obtained  by  the  following  relation- 
ship: 

F  =  Total    Interest  Payments3  -   Financed  Capital    x   Interest  Rate 
Total    Land^  x  Average  Land  Valued  x   Interest  Rate 
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Notes : 

Production   function 
A  =  6.994 
ALP  =      .2 

Depreciation  Rate 
.077 


Reference 


(***) 


(5) 


ALP  was  arbitrarily  set  at  0.2.  The  valu 
from: 


e  for  A  was  then  calculated 


ALP     (1  -  ALP) 
Y  =  A  x  TK    x  AL 
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POLICY/SCENARIO    INPUTS 


INPUT 


I    Propensity  to    invest 

Salvage  value 
I   Depreciation    raie 

Investment    leverage 

1 


Capi  t  a  I 
Forma  t  ion 


T 


Net 

Farm 

Receipts 


OUTPUT 


Capita 
Stock 


INPUT 

Agricultural    land 
labor    expense 


1 


Agricultural 
Production 


OUTPUT 


Gross    production 


I    Property    tax    rate 
Interest     rate 

I 


Finure  1 
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Forestry  Sector 

The  principal  assumption  in  the  forestry  sector  is  that  Montana's 

economy  is  driven  in  the  long  run  by  the  projected  inventory  of  standing 

mature  sawtimber  available  to  be  cut  in  each  year.'  The  total  capital 

(TCR)  to  harvest  the  allowed  timber  cut  is  adjusted  by  technology  and 

the  future  capital  required.  Additional  capital  (ACN)  is  a  relationship 

between  the  total  capital  required  and  the  present  level  of  capital  (CAP) 

[a  function  of  depreciation  (DEP)  and  the  capital  formation  rate  (CFR)]. 

In  mathematical  form: 

ACAP  =  CAP  +  (CFR  -  DEP) 
ACN  =  TCR  -  CAP 


where 


ACFR  =  (ACN/CT)  +  DEP 

CT  =  investment  closure  time 


The  annual  cutting  rate  is  then  determined  by  calculating  the  difference 
between  the  standing  timber  available  to  cut  and  the  maximum  possible  cut 
given  the  current  level  of  installed  capital. 

The  structure  of  the  model  is  shown  in  Figure  1. 
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Notes 


1. 


2. 


Year 

Capital   Formation 

Employees 

Growth  Rate 

1967 
1963 
1958 
1954 

$6.4  x  106 
$9.8  x  106 
$6.8  x  106 
$6.7  x  106 

8.5/106$ 
8.3/106$ 
6.7/106$ 
6.2/106$ 

1% 

42 
2% 

7.4  million/yr.   Average 

2.5%  Average 

361.3  x  109  BF  Total  Timber 
1336  Board  Feet  x  106  Timber  Harvest  1963  (2) 


3.  depreciation  6.7%  yr.  (assumed) 

4.  capital  investment  =  actual  capital  (depreciation  +  growth  rate) 

AC  =  7.4  x  106  =  $80  x  106  actual  capital,  1963  ** 
.092 

5.  capital  needed  to  cut  1  BF/yr.  =  actual  capital/timber  harvest 

.06  $BF/yr. 


*  "Forest  Industries,"  Montana  Data  Book,  Chapter  9,  Montana  Department  of 
Planning  and  Economic  Development,  State  of  Montana,  1970. 


** 


0.092  is   the  depreciation  rate  plus   the  average  growth   rate. 
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Coal  Mining  and  Non-Energy  Mining  Sectors 

The  sectors  for  coal  mining  and  non-energy  mining  are  similar  in 
structure.  Only  some  parameters  are  different.  For  the  sake  of  simplicity, 
they  will  be  discussed  together  here. 

Figure  1  shows  the  structure  of  the  non-energy  mining  sector.  The 
total  quantity  which  has  been  mined  is  used  to  determine  the  productivity 
or  cost  for  mining.  The  mining  cost  function  relates  the  average  variable 
production  cost  per  ton  of  output  to  the  total  mined.  The  capital  cost 
function  relates  the  investment  required  for  a  given  rate  of  output  to 
the  total  mined.  The  actual  output  is  limited  by  the  total  existing  invest- 
ment. It  may  also  be  limited  by  the  product  price.  The  cost  distribution 
function  describes  the  distribution  of  cost  for  different  levels  of  output. 
If  the  price  is  divided  by  the  average  cost  then  this  function  can  deter- 
mine what  fraction  of  the  total  possible  output  is  economical.  This  is  then 
used  to  determine  the  actual  output. 

The  total  output  is  multiplied  by  the  product  price  to  determine  mining 
income.  The  average  cost  is  modified  by  integrating  the  cost  distribution 
function  over  the  fraction  actually  produced  to  obtain  a  multiplier.  The 
result  is  multiplied  by  the  actual  output  to  determine  the  mining  expense. 
The  capital  expense  is  determined  by  assuming  the  capital  investment  is 
amortized  over  a  20  year  time  period  at  a  specified  interest  rats.  The 
expenses  are  subtracted  from  the  income  to  determine  the  profit  rate  and 
the  rate  of  return.  The  investment  function  determines  the  rate  of  invest- 
ment from  the  rate  of  return. 

For  non-energy  mining,  most  ore  deposits  (copper  in  particular)  are 
most  easily  mined  at  a  limited  rate.  Higher  production  rates  are  attained 
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only  with  increased  cost  per  unit  of  output.  This  is  simulated  in  the  non- 
energy  mining  sector  by  making  output  capability  a  non-linear  function  of 
total  investment. 

The  model  for  the  coal  mining  sector  of  Montana's  economy  is  shown  in 
Figure  2.  Since  Montana  has  large  coal  reserves  (even  in  comparison  to 
total  United  States  coal  demand)  over  production  and  market  saturation 
could  conceivably  occur.  This  would  tend  to  depress  the  selling  price  of 
the  coal .  The  effect  of  market  saturation  is  simulated  by  determining  a 
price  for  the  coal  being  produced.  An  exogenous  market  level  indicates  the 
size  of  the  market  that  is  available  before  any  saturation  takes  place. 
If  more  coal  is  produced  than  can  be  sold  on  this  market,  the  price  of 
coal  is  decreased.   If  the  coal  production  then  falls  off  to  a  level  that 
is  less  than  the  market  saturation  level,  the  coal  price  will  return  to  a 
normal  value. 

In  a  growing,  profitable  industry,  such  as  coal  mining  is  in  Montana, 
there  will  undoubtedly  be  new  companies  entering  the  industry.  Thus,  new 
investment  will  come  not  only  from  the  profits  of  the  industry,  but  also 
from  other  sources.  The  additional  investment  is  simulated  as  being 
proportional  to  the  expected  rate  of  return  and  the  difference  between 
current  production  and  the  production  level  that  would  saturate  the  market. 
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Oil  and  Gas  Sector 

The  oil  and  gas  sector  consists  of  parallel,  coupled  models;  one  for 
oil  discovery  and  production  and  one  for  gas  discovery  and  production.  The 
models  are  coupled  since  the  production  of  oil  normally  has  some  gas  pro- 
duction with  it  and  the  production  of  fields  which  are  primarily  for  gas 
usually  results  in  the  production  of  some  liquids.  Similarly,  when  explor- 
ation for  oil  fields  is  carried  out,  there  is  some  probability  of  finding 
a  gas  field  or  vice  versa.  The  production  part  of  the  model  is  shown  in 
Figure  la  and  the  exploration  part  in  Figure  lc.  The  following  discussion 
deals  only  with  the  production  and  exploration  for  oil  since  the  structure 
is  the  same  for  gas. 

When  an  oil  or  gas  field  is  discovered,  the  income  obtained  from  the 
field  is  spread  over  a  long  period  of  time,  as  are  the  expenses  of  developing 
and  producing  the  field.  The  value  of  these  cash  flows  must  be  discounted 
if  the  profit  for  discovering  the  field  is  to  be  accurately  determined. 
This  is  shown  in  Figure  lb. 

When  oil  is  discovered  it  is  placed  in  the  producing  level.  Within  this 
level,  information  is  retained  as  to  the  age  of  each  discovery.  Four  functions 
are  required  to  describe  the  performance  of  producing  reservoirs.  Figure  2 
shows  the  production  function.  This  function  traces  the  rate  of  production 
from  a  typical  reservoir  starting  with  the  discovery  of  the  field,  moving 
through  the  development  of  the  field  to  a  period  of  relatively  steady-state 
production,  and  then  on  to  a  long  period  of  steady  decline  as  the  reservoir 
is  depleted.  These  production  rates  are  all  based  on  the  amount  of  initial 
recoverable  oil  in  place.  Therefore,  if  the  total  amount  of  oil  for  fields 
of  each  age  is  multiplied  by  the  production  rate  for  that  age  and  summed  for 
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fields  for  all  ages,  the  total  production  rate  at  any  point  in  time  can  be 
calculated. 

Figure  3  shows  the  operating  expense  function.  This  function  shows  the 
operating  expenditures  required  to  operate  a  producing  field  of  a  given  age. 
This  function  is  also  based  on  the  initial  recoverable  oil  in  place.  The 
total  operating  expense  for  oil  production  can  be  determined  by  multiplying 
the  amount  of  oil  for  fields  of  each  age  by  the  rate  at  which  expenditures 
are  required  for  that  age. 

Figure  4  shows  the  investment  function.  This  function  gives  the  rate  at 
which  investment  must  be  made  in  order  to  develop  and  maintain  a  field.   In  the 
early  years  of  the  field's  production,  the  investment  rate  is  quite  high  since 
the  field  is  under  development.  In  later  years  the  rate  drops  to  a  low  level 
since  little  additional  equipment,  etc.  are  required.  The  operating  expense 
function  and  the  investment  function  combine  to  give  the  total  expense  for  pro- 
duction. As  before,  the  total  rate  of  investment  is  determined  by  multiplying 
the  investment  rate  for  each  age  by  the  amount  of  oil  for  that  age. 

Figure  5  shows  the  variable  cost  of  production.   It  is  similar  to  the 
operating  expense  function  except  that  it  describes  the  cost  on  the  basis  of 
unit  of  production  for  fields  of  each  age.  This  function  is  used  to  compare 
the  price  to  the  variable  cost.  Normally,  as  a  reservoir  is  depleted  the  costs 
per  unit  of  production  increase.  Thus,  the  price  will  determine  how  long  fields 
are  produced  before  being  abandoned.  A  cut-off  age  is  determined  where  the 
price  equals  the  cost.  Fields  older  than  this  are   removed  from  production. 

The  total  oil  production  is  the  sum  of  the  oil  produced  in  oil  fields  plus 
the  liquids  produced  in  natural  gas  fields.  Similarly,  with  the  oil  production, 
a  certain  amount  of  gas  is  produced  determined  by  the  average  gas/oil  ratio. 
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The  total  income  for  oil  production  is  then  the  oil  produced  in  oil  fields 
multiplied  by  the  oil  price  plus  the  gas  produced  in  the  oil  fields  mul- 
tiplied by  the  wellhead  gas  price.  The  expenses  for  the  production  come 
from  the  total  operating  expenses  plus  the  costs  for  capital.  Capital  is 
assumed  to  be  amortized  at  a  specified  interest  rate,  over  a  20  year  time 
span.  The  net  profit  rate  per  barrel  of  oil  is  determined  by  subtracting 
the  expenses  from  the  income  and  dividing  by  the  production  rate. 

The  rest  of  the  model  deals  with  exploration.  Two  additional  functions 
are  required  here.  Figure  6  shows  the  depletion  function.  This  function 
gives  the  increase  in  drilling  required  for  additional  discoveries  as  more 
and  more  of  the  oil  is  discovered.  Figure  7  shows  the  investment  function. 
This  function  gives  the  rate  of  investment  in  terms  of  the  profit  rate.  It 
is  similar  to  the  same  function  for  other  industries. 

From  the  depletion  function,  the  oil  discovered  per  foot  of  drilling 
is  determined.  From  the  investment  in  exploration  drilling,  the  drilling 
cost,  and  the  drilling  investment  life,  the  rate  of  drilling  is  determined. 
These  combine  to  give  the  direct  discovery  rate  of  oil.  The  oil  discovered 
in  gas  exploration  is  added  to  determine  the  total  discovery  rate  for  oil. 
The  income  for  oil  exploration  is  determined  by  multiplying  the  direct 
discovery  rate  for  oil  by  the  oil  production  profit  and  adding  the  gas  dis- 
covered in  oil  exploration  multiplied  by  the  gas  production  profit.  The  ex- 
pense for  the  exploration  is  simply  the  cost  for  the  investment  amortized 
over  the  investment  life.  The  income  and  expense  are  compared  to  determine 
the  profit  rate.  This  is  used  with  the  investment  function  to  determine  the 
investment  rate. 
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The  following  list  summarizes  the  inputs,  parameters,  and  key  outputs 
for  the  oil  and  gas  model.  Figures  2-6  show  functions  required  for  the  oil 
part  of  the  model.  Similar  functions  are  required  for  the  gas  part. 
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List  of  Inputs,  Parameters,  and  Key  Outputs 
for  Oil  and  Gas  Production  and  Exploration 


Inputs: 

Price  of  oil  at  wellhead  ($/bbl) 
Price  of  gas  at  wellhead  ($/ft3) 

Parameters : 

Average  producing  gas/oil  ratio  for  oil  fields  (ft  /bbl) 
Average  producing  oil/gas  ratio  for  gas  fields  (bbl/ft^) 
Capital  cost  for  production  investment  ($/$capital/Yr) 
Gas  fields  discovered  in  oil  exploration  (ft  /bbl) 
Oil  fields  discovered  in  gas  exploration  (bbl/ft^) 
Exploration  drilling  cost  ($/ft) 
Exploration  capital  life  (Yr) 
Exploration  capital  cost  ($/$capital/Yr) 

Key  Outputs: 

Total  oil  production 
Total  gas  production 
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II 


Petroleum  Refining  Sector 

The  petroleum  refining  in  Montana  is  related  to  the  output  of  the 
petroleum  production  sector.   It  is  also  dependent  upon  imports  of  crude 
oil.  These  factors  are  more  important  than  the  actual  demand  for  refinery 
products  in  Montana.  Therefore,  refining  is  considered  separately  rather 
than  with  oil  production  or  oil  demand  in  Montana. 

The  vast  majority  of  the  refinery  capacity  in  Montana  is  located  around 
the  Billings  area.  However,  the  pipeline  system  is  such  that  the  oil  produced 
in  the  eastern  part  of  the  state  is  not  available  to  these  refineries. '^ ) 
In  1975  Montana  produced  32.8  x  10°  barrels  of  crude  oil,  refined  48.1  x  10" 
barrels  of  oil,  but  refined  only  8.9  x  10  barrels  of  oil  which  was  produced 
in  Montana. ^    In  1974  Montana  produced  34.6  x  106  barrels.  Of  this,  10.2  x 
10^  barrels  was  from  central  Montana  and  thus  available  to  Montana  refineries; 


24.4  x  10  barrels  was  produced  in  eastern  Montana  and  was  not  available  to 

•iso 
(2) 


the  refineries.    By  comparison,  9.3  x  10  barrels  of  Montana  production  were 


actually  refined  in  Montana. 

Of  the  imports  refined  in  1975,  19.8  x  10  barrels  came  from  Wyoming 
and  19.4  x  10  barrels  came  from  Canada.  It  is  difficult  to  predict  the 
availability  of  these  crude  oil  imports  in  the  future,  especially  the  Canadian 
imports.  Since  the  companies  who  operate  the  pipelines  and  refineries  are 
interstate  in  nature,  it  is  difficult  to  predict  whether  or  not  they  would 
reverse  past  practices  and  develop  facilities  to  transport  eastern  Montana 
r  ybduction  to  Montana  refineries  if  imports  could  not  be  obtained. 

The  model  for  petroleum  refining  assumes  only  central  Montana  production 
is  available  to  the  refineries.   In  addition,  a  scenario  input  of  imports 
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available  must  be  provided.  It  is  assumed  that  additional  capacity  will  not 
be  added  unless  Montana  production  warrants  it.  This  is  an  unlikely  event. 
Thus,  the  maximum  refinery  output  will  depend  upon  the  Montana  production 
which  is  available  plus  the  imports  available  or  the  maximum  output  of 
current  capacity.  Only  if  available  Montana  production  should  exceed  the 
refinery  capacity,  will  more  capacity  be  installed.  Since  the  oil  and  gas 
model  gives  oil  production  only  in  terms  of  the  total  for  Montana,  a 
parameter  which  states  the  fraction  which  is  available  to  Montana  refineries 
is  required.  Using  1974  data  from  Reference  1,  the  fraction  of  production 
from  the  central  part  of  Montana  is  approximately  29%.  From  Reference  2, 
the  total  refinery  capacity  in  1975  was  161,000  bbl/day. 
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Primary  Metals  Sector 

The  output  of  the  primary  metals  sector  of  the  economy  is  dependent 
primarily  upon  the  output  of  the  metal  mining  sector.  It  is  also  very 
energy  intensive.  For  the  present,  the  primary  metals  output  is  assumed 
to  be  proportional  to  the  non-energy  mining  output.  Since  it  is  a  large 
energy  user,  it  may  be  desirable  to  consider  it  in  more  detail  later  on. 
The  ratio  of  the  outputs   from  Reference  1    is: 

Primary  Metals  Output  439.7 


Non-Energy  Mining  Output  "    177.0 


2.48 
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References  for  the  Economic  Sub-Models 

Agricultural  Sector 

1.  Economic  Component,  Agricultural  Sector,  Oregon  State  Simulation  Model, 
Willamette  Simulation  Unit,  Oregon  State  University. 

2.  Entire  paragraph  was  adapted  from  the  agricultural  sector  of  the  Economic 
Component,  Oregon  State  Simulation  Model,  Willamette  Simulation  Unit, 
Oregon  State  University. 

3.  Montana  Data  Book  (Montana  Department  of  Planning  and  Economic  Develop- 
ment)  Chapter  10  Agriculture. 

4.  U.S.  Department  of  Commerce,  Census  of  Agriculture-Montana,  October  1971. 

5.  H.  H.  Fullerton,  J.  R.  Prescott,  An  Economic  Simulation  Model  for  Regional 
Development  Planning,  Ann  Arbor  Science. 

Forestry  Sector 

1.  "Economic  Component,"  Logging  Sector:  Willamette  Simulation  Unit, 
Oregon  State  University. 

2.  "Forest  Industries,"  Montana  Data  Book,  Chapter  9,  Montana  Department  of 
Planning  and  Economic  Development,  State  of  Montana,  1970. 

Coal  Mining  and  Non-Energy  Mining  Sector 

1.  Montana  University  Coal  Demand  Study  Team,  Projections  of  Northern  Great 
Plains  Coal  Mining  and  Energy  Conversion  Development  1975-2000,  May  1976. 

2.  Montana  Environmental  Quality  Council,  Montana  Energy  Policy  Study, 
February  1975. 

3.  Montana  Energy  Advisory  Council,  Montana  Historical  Energy  Statistics, 
September  1976. 

Oil  and  Gas  Sector 

1.  "Annual  Review  for  the  Year  1974  Relating  to  Oil  and  Gas,"  Volume  18, 
Oil  and  Gas  Conservation  Division,  Department  of  Natural  Resources  and 
Conservation,  State  of  Montana. 

2.  Gustav  Stolz,  Jr.,  Oil  and  Gas  Energy  Resources  of  Montana,  Montana 
Environmental  Quality  Council,  1974. 

3.  Maher  Ibrahim,  Undiscovered  Natural  Gas  Resources  of  Montana,  Montana 
Energy  Advisory  Council,  1976. 
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References  for  the  Economic  Sub-Models  (cont'd) 

Petroleum  Refining  Sector 

1.  "Annual  Review  for  the  Year  1974  Relating  to  Oil  and  Gas,"  Volume  18, 
Oil  and  Gas  Conservation  Division,  Department  of  Natural  Resources  and 
Conservation,  State  of  Montana. 

2.  "Historical  Petroleum  Statistics  for  Montana,"  Montana  Energy  Advisory 
Council,  Office  of  the  Lieutenant  Governor. 


Primary  Metals  Sector 

1.  Input/Output  Model,  Montana  Futures  Project,  Division  of  Research  and 
Information  Systems,  Department  of  Community  Affairs. 


A-44 


Section  A-V,  Economic  Input/Output  Sub-Model 

The  economic  model  is  an  input/output  model  for  Montana.'  The  out- 
puts for  agriculture,  forestry,  coal  mining,  oil  and  gas,  non-energy 
mining,  petroleum  refining,  and  primary  metals  are  all  determined  elsewhere. 
The  direct  input/output  coefficients  are  used  to  determine  the  demands  for 
other  products  due  to  these  outputs.  These  demands  are  then  treated  as 
demands  exogenous  to  the  input/output  model.  Additional  exogenous  demands 
arise  from  consumer  spending,  exports,  and  government  spending.  Each  of 
these  will  be  discussed  separately. 

Consumer  Spending 

Total  consumer  spending  is  assumed  to  be  proportional  to  earned  income. 
In  effect,  the  elasticity  is  assumed  to  be  unity.  The  consumer  spending  is 
assumed  to  be  divided  among  the  sectors  of  the  economy  as  in  the  input/output 
model . 

The  total  non-agricultural  earnings  (in  1972  dollars)  was  2,036  million 
dollars  in  1970  and  2,392  million  dollars  in  1974,  for  an  average  of  2,214 
million  dollars.  It  is  more  difficult  to  estimate  earned  income  in  the 
agricultural  sector.  An  estimate  of  approximately  350-360  million  dollars 
in  1972  was  made.  This  gives  an  initial  estimate  of  2,570  million  dollars 
of  earned  income  in  1972.  The  consumer  spending  in  1972  was  1,629  million 
dollars.'  This  gives  a  spending/earning  ratio  of  0.63. 


Exports 

It  is  not  necessary  to  consider  exports  for  the  sectors  considered  as 
exogenous  to  the  input/output  model.  Of  the  endogenous  sectors,  only  meat 
processing,  railroads,  communications,  utilities,  and  retail  trade  have 
significant  exports.  The  railroad  sector  is  the  only  one  for  which  exports 

A- 4  5 


dominates.  For  the  others  a  constant  ratio  between  consumer  demand  and  ex- 
port demand  is  assumed. 

In  1972  $104  x  106  out  of  a  total  railroad  output  of  $179  x  106  was 
exported.   It  is  difficult  to  predict  how  this  export  demand  may  increase  in 
the  future.  If  there  is  a  large  development  in  the  coal  industry  (which  is 
quite  likely)  it  would  directly  affect  the  export  demand  for  rail  transpor- 
tation. The  current  cost  of  coal  transport  by  rail  is  approximately  lt/ton- 
mile.  If  exported  coal  travels  an  average  of  100  miles  in  Montana,  there  is 
approximately  $l/ton  rail  export  which  results. 

The  exports  demand  for  rail  transportation  is  currently  calculated  in 
the  model  by  the  following  equation. 

Rail  Export  Demand  =  $100  x  106/yr  +  $l/ton  x  Coal  Output 

Government  Spending 

Government  spending  for  both  federal  and  state  and  local  governments 
is  assumed  to  be  proportional  to  earned  income.   Federal  spending  in  1972 
was  1,092  million  dollars.'1'  Using  the  previous  estimate  for  earned  income, 
the  respective  spending/earnings  ratios  are  0.42  and  0.25  (federal,  state, 
and  local ) . 
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Economic  Input/Output  Sub-Model  References 

1.  Input/Output  Model,  Montana  Futures  Project,  Research  and  Information 
Systems  Division,  Department  of  Community  Affairs,  State  of  Montana. 

2.  M.  C.  Johnson,  "Montana's  Economy:  Where  It's  Been  and  Where  It's 
Going,  "Montana  Business  Quarterly,  Vol.  14,  No.  1,  Winter  1976. 
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Section  A-VI,  Capital  Investment  Sub-Model 

The  structure  of  the  capital  investment  model  for  a  single  sector  is 
shown  in  Figure  1.  The  same  structure  is  used  for  all  non-primary  sectors 
The  capital  investment  for  the  primary  sectors  is  determined  in  their  indi- 
vidual models. 

The  capital  investment  model  uses  exponential  delays  to  average  both 
the  output  and  the  rate  of  change  for  a  sector.  The  averaged  rate  of 
change  is  added  to  the  depreciation  rate  to  get  the  total  investment  rate 
required.  However,  if  this  rate  is  negative,  there  is  no  negative  invest- 
ment.  If  it  is  positive,  it  is  multiplied  by  the  averaged  output  and  the 
capital  output  ratios  and  depreciation  rates  are  taken  from  Reference  1. 
They  are  listed  in  Table  1.  An  "averaging  delay  time"  of  5  years  is  used. 
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Table  1.  Capital-Output  Ratios  and  Depreciation  Rates 


Sector 

Capital-Output 
Ratio 

Depreciation 
Rate 

Construction 

0.19 

0.093 

Food  and  Kindred  Products 

0.25 

0.065 

Printing  and  Publishing 

0.39 

0.066 

Chemicals  and  Allied  Products 

0.50 

0.067 

Other  Manufacturing 

0.59 

0.082 

Transportation 

1.76 

0.038 

Communications  and  Utilities 

3.00 

0.033 

Trade 

0.65 

0.052 

Finance,  Insurance,  Real  Estate 

1.05 

0.047 

Services 

0.95 

0.097 

Government 

0.95 

0.097 
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Capital    Investment  Sub-Model    References 

1.      H.   H.    Fullerton  and  J.    R.    Prescott,  An   Economic  Simulation  Model   for 
Regional    Planning,  Ann  Arbor  Science  Publishers   Inc.,  Ann  Arbor, 
Michigan  1975. 
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Section  A- VI  I,  Employment  Sub-Model 

The  structure  of  the  employment  model  is  shown  in  Figure  1.  Employment 
is  divided  into  three  groups  along  the  same  lines  as  in  Reference  1.  These 
groups  and  the  sectors  included  are: 

1 .  Agriculture  - 

2.  Non-Agricultural    Primary- 

Mining 

Manufacturing 

Railroads 

Federal  Government 

3.  Derivative  - 

Utilities 

Construction 

Trade 

Services  and  Finance 

State  and  Local  Government 

Other 

This  division  is  made  to  account  for  the  fact  that  wage  rates  for  the  dif- 
ferent groups  vary  considerably. 

The  employment  levels  are  determined  by  projections  of  the  employee/out- 


put ratios. 


(2) 


These  are  presented  in  Table  1.  A  linear  interpolation  is 


made  for  years  between  those  given.  After  1985  they  are  assumed  to  stay 
constant.  Employee/output  ratios  are  not  given  for  agriculture  or  government. 

Government  employment  should  relate  to  government  spending  in  much 
the  same  way  as  private  employment  relates  to  output.  In  1970  the  federal 
government  provided  11,900  jobs  and  in  1974,  13,000  jobs.'  '  From  reference 
3,  the  federal  government  spent  $1,092  x  10  in  1972.  Using  the  average 
employment  gives  an  employee/spending  ratio  of  11.4  employees/million  dollars. 


In  1970  the  state  and  local  governments  provided  40,700  jobs  and  in 


1974  they  provided  45,200  jobs 


(1) 


In  1972  state  and  local  governments 
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spent  $655  x  10  .    Using  the  same  method  as  above,  the  employee/spending 
ratio  for  state  and  local  governments  obtained  is  65.6. 

As  can  be  seen  in  Table  1  the  employee/output  ratios  are  expected  to 
decline  in  the  private  sector.  Similar  declines  can  probably  be  expected 
for  government  employment.  The  above  coefficients  are  arbitrarily  assumed 
to  decline  by  20%  in  1980,  by  26%  in  1985,  and  stay  constant  thereafter. 
This  is  somewhat  less  decline  than  most  of  the  ratios  for  the  private  sector. 

The  employment  in  the  agriculture  sector  is  treated  similarly.  The 
total  agriculture  output  in  1972  was  $875.3  x  lOv  '     The  agriculture  employ- 
ment was  36,100  in  1970  and  35,100  in  1974. 0)  Using  the  same  method  as 
before,  the  employee/output  ratio  for  agriculture  obtained  is  40.7.  There 
have  been  large  decreases  in  the  agriculture  employee/output  ratio  in  the 
past.  However,  it  is  assumed  that  the  decreases  will  be  more  moderate  in 
the  future.  The  ratio  is  arbitrarily  assumed  to  decrease  by  25%  in  1985, 
by  35%  in  1980,  and  to  hold  constant  thereafter. 

The  employment  for  the  various  sectors  are  summed  to  determine  the  total 
employment.  The  total  work  force  available  is  determined  by  multiplying 
the  total  population  as  supplied  by  the  demographic  model  by  the  employee/ 
population  projection  shown  in  Table  2.  These  ratios  are  obtained  by  altering 
employment  population  ratios  from  Reference  2  to  account  for  the  unemploy- 
ment considered  in  these  ratios.  The  total  employment  and  the  total  work 
force  are  compared  to  determine  the  unemployment  rate. 

The  wage  rates  for  the  non-agricultural  primary  and  the  derivative 
groups  increase  in  real  terms  if  unemployment  drops  below  5%.  When  unem- 
ployment is  higher  than  this  the  wage  rates  are  assumed  to  stay  constant. 
The  wage  rates  for  the  agriculture  sector  are  somewhat  different  since  a 
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large  part  of  the  employment  consists  of  farm  operators.  For  them,  the 
wage  rate  depends  upon  the  net  farm  receipts  which  are  not  reinvested.  For 
hired  hands,  the  wages  are  assumed  to  change  in  the  same  manner  as  for  the 
other  groups.  The  fraction  of  employees  which  are  hired  must  be  supplied 
as  a  parameter.  From  Reference  4  the  average  number  of  hired  farm  workers 
was  11,000  in  1969.  Compared  to  the  total  employment  of  36,100  in  1970 
thisrgives  a  fraction  of  about  30  ',.      This  fraction  is  assumed  to  state  constant. 

In  order  to  start  the  simulation,  it  is  necessary  to  have  initial 
values  for  the  wage  rates.  The  wage  rate  for  derivative  workers  was  $5,072/yr. 
in  1974;'''  for  non-agricultural  primary  workers  it  was  $9,480/yr  in  1974 
(^)  (both  in  terms  of  1967  dollars).  The  wage  rate  for  hired  farm  labor 
was  $3,340  in  1969  as  determined  for  the  agriculture  model. 
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Table  1.  Employee/Output  Ratios' 


Employee  /$10b  Output 
1972     1980     1985 


Sector 


25 

18 

16 

11 

9 

7 

14 

13 

12 

33 

24 

20 

16 

13 

11 

30 

21 

17 

5 

4 

4 

9 

7 

6 

28 

23 

21 

37 

22 

17 

38 

33 

14 

32 

21 

17 

14 

10 

8 

78 

64 

58 

17 

14 

13 

118 

100 

92 

Non-fuel  Mining 

Coal 

Crude  Petroleum 

Construction 

Food  Processing 

Lumber  and  Products 

Petroleum  Refining 

Primary  Metals 

Other  Manufacturing 

Railroad 

Motor  Freight 

Communication 

Util ities 

Trade 

Finance,  Insurance,  Real  Estate 

Service 


*  All  in  1972  dollars. 
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Table  2.  Projection  of  Work  Force/Population  Ratio 
Year         


1975 
1980 
1985 

1990 


Work  Force/ 
Population 


0.48 
0.51 
0.515 
0.52 
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Employment  Sub-Model  References 

1.  M.  C.  Johnson,  "Montana's  Economy":  Where  It's  Been  and  Where  It's 
Going",  Montana  Business  Quarterly,  Vol.  14,  No.  1,  Winter  1976. 

2.  Montana  Futures  Project,  Division  of  Research  and  Information  Systems, 
Department  of  Community  Affairs,  State  of  Montana. 

3.  Input/Output  Model  -  Reference  2. 

4.  "Agriculture,"  Montana  Data  Book,  Chapter  10,  Department  of  Planning 
and  Economic  Development,  State  of  Montana  1970. 
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Section  A- VI II,  Energy  Sub-Model 

The  energy  sub-model  is  divided  into  demand  and  supply  sections.  The 
demand  section  is  fairly  detailed  while  the  supply  section  is  fairly  simple 
and  obtains  a  large  part  of  its  information  directly  from  the  economic  sub- 
model. Each  section  will  be  discussed  separately. 

Energy  Demand 

The  demand  for  energy  is  analyzed  along  the  traditional  lines  of  resi- 
dential, commercial,  and  industrial  usage.   Each  part  of  the  demand  is 
treated  somewhat  differently  and  is  discussed  separately.  Agricultural 
energy  demand  is  considered  separately  also. 

Residential  Demand 

The  demand  for  residential  energy  is  divided  among  its  various  end  uses. 

These  are: 

Space  Heating 
Water  Heating 
Kitchen  Ranges 
Refrigeration 
Clothes  Drying 
Air  Conditioning 
Other 

The  actual  residential  demand  for  a  particular  form  of  energy  for  a 

particular  end  use  is  determined  by  the  following  relationship: 

E-  .=A.  .  x  S-  •  x  TH 

where 

E.  .  is  the  consumption  of  energy  type  i  by  end  use  j. 
i » J 

A.j  j  is  the  average  consumption  of  energy  for  this  type  of  end  use  in 
'  a  typical  household. 

S.j  j  is  the  saturation  of  type  j  end  use  using  energy  type  i,  and 

TH  is  the  total  number  of  households. 
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Current  average  uses  and  saturations  form  a  baseline  rate  of  energy  use. 
Energy  use  can  change  from  this  baseline  in  three  v/ays : 

1.  Change  in  average  energy  consumption  for  an  end  use, 

2.  Change  in  the  relative  saturations  of  an  end  use  between  fuels,  and 

3.  Change  in  the  total  saturation  of  an  end  use. 

There  are  several  factors  which  must  be  considered  in  determining  these 
changes : 

1.  Price  effect  on  energy  consumption/end  use, 

2.  Price  effect  on  substituting  one  type  of  energy  for  another. 

3.  Price  effect  on  total  saturation  of  an  end  use, 

4.  Income  effect  on  energy  consumption/end  use, 

5.  Income  effect  on  total  saturation  of  an  end  use,  and 

6.  Trends  from  a  "market  entry  effect". 

Current  investigation  has  focused  on  determining  the  baseline  coefficients 
The  derivation  of  these  follows. 


Yearly  energy  consumption  data  for  the  electric  end  uses  were  taken  from 
reference  1.  This  was  used  unaltered  except  for  air  conditioning  and  heating. 
Heating  consumption  was  determined  by  multiplying  the  average  consumption  in 
Los  Angeles  by  the  ratio  of  Montana  to  Los  Angeles  heating  degree  days  obtained 
from  reference  2.  Due  to  the  relatively  cool  climate  in  Montana,  energy  con- 
sumption by  air  conditioners  was  assumed  to  be  one-third  of  a  typical  California 
window  unit. 

To  determine  the  annual  energy  consumption  for  non-electric  end  uses, 
the  energy  consumption  by  similar  electric  end  use  was  multiplied  by  the 
ratio  of  electric  efficiency  divided  by  the  efficiency  for  using  the  other 
energy  source.  These  efficiencies  were  obtained  from  reference  3. 
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The  saturation  of  each  type  of  end  use  for  each  fuel  was  assumed  to  be 
the  same  as  those  reported  in  reference  4.  The  total  number  of  households 
in  Montana  was  assumed  to  be  7%  greater  in  1974  than  reported  in  reference 
5  in  1970.  This  is  approximately  equal  to  the  total  population  growth  in 
Montana  during  the  same  period  as  reported  in  reference  6.  The  total  energy 
consumed  in  1974  was  then  estimated  by  multiplying  the  energy  per  end  use  by 
saturation  and  then  by  the  total  number  of  households.  The  results  were  com- 
pared with  the  1974  data  for  total  residential  electricity  and  natural  gas 
consumption  in  references  7  and  8.  Reference  7  was  assumed  to  be  more  accurate. 
An  "other"  category  was  added  for  electricity  to  account  for  the  large  number 
of  small  end  uses  not  considered  initially.  This  was  assumed  to  be  the 
difference  between  the  total  predicted  in  the  initial  estimate  and  the  total 
reported  in  reference  7  less  the  allowance  for  agriculture  in  reference  9. 
There  are  very  few  residential  uses  for  natural  gas  not  considered  in  the  ini- 
tial estimate.  Therefore,  the  energy  consumption  for  each  end  use  was  scaled 

to  make  the  total  agree  with  reference  7.  Coal  and  oil  were  left  the  same 
as  in  the  initial  estimate.  End  uses  of  bottled  gas  were  assumed  to  be  the 
same  as  for  natural  gas. 
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Table  1 .--Residential  Energy/End  Use  and  Saturation 


Bottled  Gas 


Use 

Electricity 
KWH/yr 

Natural  Gas 
103Btu/yr 

Oil 
103Btu/yr 

Coal 
103Btu/yr 

&  Other 
103Btu/yr 

Space  Heating 

27,500 
2.3% 

118,700 
83.8% 

141,300 
7.1% 

161,900 
1.9% 

118,700 
4.9% 

Water  Heating 

4,200 
17.9% 

20,700 

78.7% 

26,500 
0% 

88,300 
0% 

20,700 
2.9% 

Kitchen  Range 

1,180 
78.2 

8,130 
19.3% 

8,130 
2.3% 

Refrigerator 

1,140 

105.1% 

Freezer 

1,200 
65.3% 

Clothes  Dryer 

990 
68 . 8% 

4,110 
3.9% 

4,110 
0.3" 

Air  Conditioner 

460 

13.4% 
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Table  2. --Average  Residential  Energy  per  Household  by  End  Use* 


Bottled  Gas 


Use 

Electricity 
KWH 

Natural  Gas 
103Btu 

Oil 
103Btu 

Coal 
103Btu 

&  OLther 
lO^Btu 

Space  Heating 

63 

99,500 

10,000 

3,080 

5,820 

Water  Heating 

760 

16,300 

600 

Kitchen  Range 

920 

1,570 

190 

Refrigerator 

1,200 

Freezer 

780 

Clothes  Dryer 

680 

160 

12 

Air  Conditioner 

62 

Total  per 
Household 

5,030 

117,500 

10,000 

3,080 

6,620 

Total  Montana** 

1,170  x  106 

2,300  x  105 

2,300  x  106 

720  x  106 

1,540  x  106 

*   (Energy  per  End  Use)  x  saturation 
**  Assumes  232,500  households 
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Table  3. --Comparison  of  Initial  Estimate 
With  Other  Sources  of  Data 

Total  1974  Montana  Use 
109Btu 


Use 

El 

ectrici 

ty 

Natural  Gas 

Oil 

Coal 

Bottled  Gas 
&  Other 

Initial  Estimate 

4,050 

27,300 

2,300 

720 

1,540 

Montana  Energy 
Advisory  Council 
(Reference  37) 

6,400 

21,600 

— 

— 

— 

Los  Alamos 
(Reference  38) 

6,500 

23,600 

4,000 

200 

— 

Electricity  unaccounted  for  put  into  "other"  category: 

Other  Electricity= 
(6,400  -  4,050)  x  109Btu 


232,500  Households  x  3,412  Btu/KWH 

=2,960  KWH/Household  x  0.25  (Adjustment  for  Agriculture) 
=  740  KWH/Household  " 


21  ,600  =  0.79 


Total  Gas  Use  Reference  7 

Total  Gas  Use  Initial  Estimate   27,300 
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Table  4. --Baseline  Estimate  -Energy  per  End  Use 

Residential 


Use 

Electricity 
KWH/yr 

Natural  Gas 
103Btu/yr 

Oil 
103Btu/yr 

Coal 
103Btu/yr 

Bottled  Gas 
&  Other 
103Btu/yr 

Space  Heating 

27,400 

93,800 

141,300 

161,900 

93,800 

Water  Heating 

4,220 

16,300 

26,500 

88,300 

16,300 

Kitchen  Range 

1,180 

6,420 

6,420 

Refrigerator 

1  ,140 

Freezer 

1,200 

Clothes  Dryer 

990 

3,250 

3,250 

Air  Conditioner 

460 

Other 

740* 

Energy  not  accounted  for  per  household  by  above  end  uses  less  the  agriculture 
allowance.  A  100%  saturation  should  be  used. 
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Residential  Demand  References 


1.  C.  C.  Mow,  W.  E.  Mooz,  S.  K.  Anderson,  "A  Methodology  for  Projecting  the 
Electrical  Energy  Demand  of  the  Residential  sector  in  California," 
R-995-NSF/CSRA,  Rand  Corporation,  March  1973. 

2.  "Future  Gas  Consumption  of  the  United  States,"  Vol.  6,  The  University 
of  Denver  Research  Institute,  December  1975. 

3.  Stanford  Research  Institute,  "Patterns  of  Energy  Consumption  in  the 
United  States,"  Office  of  Science  and  Technology,  Executive  Office  of 
the  President,  January  1972. 

4.  "Residential  Appliance  Saturation  Survey,"  The  Montana  Power  Company, 
January  1974. 

5.  "General  Housing  Characteristics,"  1970  Census  of  Housing,  advance  re- 
port, U.S.  Department  of  Commerce/Bureau  of  the  Census,  September  1970. 

6.  M.  C.  Johnson,  "Montana's  Economy:  Where  It's  Been  and  Where  It's 
Going,"  Montana  Business  Quarterly,  Vol.  14,  No.  1,  Winter  1976. 

7.  Energy  use  summary  prepared  by  the  Montana  Energy  Advisory  Council, 
preliminary  version,  1974  data  used  in  calculations. 

8.  A.  J.  Barret,  W.  A.  Beyer,  and  C.  D.  Kolsted,  "Rocky  Mountain  Energy 
1974:  Flow,  Employment,  Prices,"  LA-6122-MS,  Internal  Report,  Los 
Alamos  Scientific  Laboratory,  October  1975. 


9.  See  section  on  agriculture  energy  usage. 
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Industrial  Demand 

Industrial  energy  demand  is  determined  separately  for  each  sector. 
The  use  of  a  particular  type  of  energy  in  a  given  sector  is  determined  by 
the  following  relationship. 

E.;  n-  =  A.  •  *  0. 
where: 

E-j  j  is  the  use  of  energy  type  ^  in  sector  j, 
A-;   i  is  the  energy  type  •  consumption  per  unit  of  output,  and 
0.  is  the  output  of  sector  ,-. 
Where  possible  it  was  desired  to  have  the  energy  use  coefficient  in 
terms  of  physical  measures  of  output  rather  than  dollars  of  output.  However, 
in  most  cases  this  was  impossible  due  to  multiple  product  from  a  single  sec- 
tor or  lack  of  data. 

For  most  sectors,  energy  use  and  output  were  determined  from  separate 
sources  and  compared  to  evaluate  the  coefficients.  In  a  few  cases  energy 
use  per  unit  of  output  could  be  determined  by  other  methods.  Table  5-8 
show  the  data  used  in  calculating  the  coefficients  and  the  coefficients 
calculated.  In  some  cases  the  energy  use  is  for  a  different  year  than  the 
output.  This  is  undesirable,  but  is  necessary  until  more  data  are  available. 
Wherever  possible  the  energy  use  and  output  were  taken  for  years  as  close 
together  as  data  allowed. 

The  coefficients  are  complete  for  gas  and  electricity.  However,  for 
coal  and  petroleum  the  data  are  poor.  Fortunately,  this  will  not  affect 
this  study  significantly  as  most  petroleum  and  SNG  will  see  different  markets, 
With  coal,  the  main  concern  is  with  the  converting  of  existing  natural  gas 
uses  to  coal . 
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Table  5.   Industrial  Use  of  Electricity 


Sector- 

Energy  Use 
(MWH/Yr) 

Output 
(106   $/Yr) 

Energy  Coefficient 
(MWH/106  $) 

Coal   Mining 

8,06c1 

(21, 

105. 02 
,1    x  106  tons)z 

76.8 
(382  MWH/10     ton) 

Metal    Mining 

306, 6791 

143. 52 
(100.7  x  10     tons 
of  copper  dominates) 

2,137 

Petroleum  and  Natural 

Gas 

191, 1911 

253. 42 

786 

Non-Mineral   Mining 

54,09s1 

17. 32 

3,127 

Food  and  Kindred  Prodi 

jets 

69,5773 

257. 74 

270 

Printing  and  Publishing 

6,4143 

32.  34 

199 

Chemical   Products 

489, 8023 

24. 64 

19,911 

Lumber  and  Wood  Product, 
Paper  and  Allied  Products 

386, 2233 

299. 04 

1,292 

Petroleum  Refining 

311, 7673 

232. 54 

1,335 

Primary  Metal 

3,829,0543 

439. 74 

8,708 

Fabricated  Metals 

580  3 

13. 94 

41.7 

Stone,   Clay,   Glass, 
and  Concrete 

81,7043 

42. 44 

1,927 

Other  Manufacturing 

1 ,744] 

255. 02 

6.8 

Agricul ture 

484, 0005 

576. 06 

840 

Utilities 







Railroads 

23,0003 

178. 74 

129 

Motor  Freight 

0 





Superscripts  refer  to  references  at  the  end  of  this  section. 
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Table  6.   Industrial  Use  of  Natural  Gas 


Energy  Use      Output      Energy  Coefficient 
Sector (Ip6ft3/Yr)     (IP6  $/Yr) (ft3/$)  


Total  Mining 

3,3677 

519. 22 

6.48 

Food  and  Kindred  Products 

1,2707 

257. 74 

4.93 

Printing  and  Publishing 

5087 

32 .  34 

15.7 

Chemical  Products 

1,1007 

24. 64 

44.7 

Lumber  and  Wood  Products, 
Paper  and  Allied  Products 

5,6007 

299. 04 

18.73 

Petroleum  Refining 

6,1 957 

233. 54 

26.5 

Primary  Metal 

13,1007 

439. 74 

29.8 

Fabricated  Metals 

4357 

13. 94 

31.4 

Stone,  Clay,  Glass, 
Concrete 

and 

4,4007 

42. 44 

103.8 

Other  Manufacturing 

1,5363 

255. 04 

6.02 

Agriculture 

05* 



0 

Utilities 

1,1097 

212. 54 

5.22 

Railroads 

0 





Motor  Freight 

0 





See  Agricultural  Energy  Demand  section. 
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Table  7.   Industrial  Use  of  Petroleum 


Sector 


Energy  Use      Output     Energy  Coefficient 
(103bb1/yr)    (10b  $/Yr)       (bb1/lp3  $) 


Railroad  Transportation 

1,8408 

178.74 

10.30 

Utilities 

258 

212. 54 

0.118 

Agricul ture 

3,090J 

576. 06 

5.36 

All    Other   Industrial 

2,2888 

1,899.84 

1.20 

Highway  Use 

(Allowing  for  Ag  Highway) 

7,1979 





See  Agricultural  Energy  Demand  Section 


A- 7  0 


■ 


■ 


H 


Table  8.   Industrial  Use  of  Coal 


Sector 


Energy  Use 
(103  tons/Yr) 


Output 
(IP6  $/Yr) 


Energy  Coefficient 


Utilities 

Other  Industrial 


76910 


318 


10 


229.9 4 
2,213.84* 


3.34 
0.144 


*  Does  not  include  agriculture 
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Industrial  Demand  References 

1.  Energy  Use  Summary  prepared  by  the  Montana  Energy  Advisory  Council,  pre- 
liminary version,  (1975  data). 

2.  "Directory  of  Mining  Enterprises  for  1975,"  State  of  Montana  Bureau  of 
Mines  and  Geology,  Bulletin  100,  May  1976. 

3.  Reference  1  (1974  data). 

4.  Montana  Futures  Project  Input  Output  Model,  Research  and  Information 
Systems  Division,  Department  of  Community  Affairs  (1972  data). 

5.  See  section  on  energy  use  in  agriculture. 

6.  "1969  Census  of  Agriculture-County  Data,"  U.S.  Department  of  Commerce/ 
Bureau  of  the  Census. 

7.  "Future  Gas  Consumption  of  the  United  States,"  Volume  6,  The  University 
of  Denver  Research  Institute,  December  1975.   (1974  data) 

8.  Reference  1  (1972  data). 

9.  Reference  1  (1969  data). 

10.  L.  H.  Crump  and  C.  L.  Readling,  "Fuel  and  Energy  Data:  United  States 
by  States  and  Regions,  1972,"  U.S.  Bureau  of  Mines. 
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Commercial  Demand 

There  appears  to  be  very  little  data  available  which  gives  details  on 

commercial  energy  demand.  This  is  a  common  problem  for  all  energy  studies. 

Usually  everything  that  is  left  over  is  lumped  into  the  commercial  category. 

According  to  reference  1,  commercial  energy  use  in  the  United  States  is 

approximated  by  the  following  proportions. 

Heating  and  Air  Conditioning  -  70% 
Lighting  -  13/o 

Refrigeration  -  4% 

Hot  Water  -  3% 

Other  -  10,' 

This  gives  no  information  concerning  the  use  of  various  fuels,  nor 

does  it  indicate  what  fraction  of  it  is  associated  with  various  activities 

(wholesale  trade,  retail  trade,  professional  services,  etc.).   It  does 

indicate  that  heating  and  air  conditioning  dominate.   In  Montana,  heating 

probably  dominates. 

Reference  2  gives  the  following  uses  of  energy  in  Montana  for  the  com- 
mercial sector: 

Gas 19.5  x  1012  Btu 

Oil n.2  x  1012  Btu 

Electricity--  4.7  x  1012  Btu 

Coal 0 

Total 35.4  x  10u  Btu 

For  comparison,  1974  data  from  reference  3  for  commercial  energy  usage 

is  shown  below.  These  are  all  residuals  after  all  other  uses  have  been 

accounted  for. 

Natural  Gas--  14.3  x  10"l2  Btu 

Electricity--  6.2  x  1 0"l 2  Btu 

Coal 2.6  x  10"12  Btu 

Oil 10.3  x  1012  Btu 

LPu 3.1  x  101 2  Btu 

Total 36.5  x  10l2  Btu 

The  discrepancies  are  probably  due  largely  to  differences  in  definitions. 

The  total  for  the  two  sets  of  data  are  within  about  3%  of  each  other.  Since 
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reference  3  was  used  in  the  residential  sector,  it  was  used  in  the  model  to 
maintain  consistency. 

If  the  relative  proportions  for  end  uses  given  previously  were  true 
in  Montana,  the  lighting  and  refrigeration  categories  would  use  6.4  x  10^  Btu; 
more  than  the  total  electricity  use  above.  Reference  4  gives  a  somewhat 
different  end  use  breakdown  for  the  commercial  sector  indicating  the  pre- 
vious numbers  are  by  no  means  certain: 

Heating  and  Air  Conditioning--  60.4 

Water  Heating 7.5 

Refrigeration 7.6 

Cooking 1  .6 

Asphalt  and  Road  Oils 11.2 

Other 11.7 

Using  this  breakdown  as  a  rough  guide  the  energy  use  in  the  commercial 

sector  was  divided  as  shown  in  Table  9.  This  table  is  arbitrary  and 

should  be  treated  as  such.  The  numbers  in  Table  9  were  then  used  to  develop 

energy  coefficients  for  the  commercial  sector.  The  output  for  the  commercial 

sector  was  also  treated  as  a  residual.   It  consists  of  all  trade,  service, 

and  government  sectors;  the  construction  sector;  and  the  communications 

sector.  From  reference  5,  the  total  output  for  these  sectors  is  1924  million 

dollars  in  1972.  From  reference  6  it  is  seen  that  income  in  Montana  increases 

approximately  5%  from  1972  to  1974.   It  was  assumed  that  commercial  output 

increased  similarly.  An  adjusted  output  of  2,021  million  dollars  was  then 

calculated.  Table  9  and  this  value  for  output  were  then  used  to  calculate 

the  energy  coefficients  shown  in  Table  10.  Energy  conversion  factors  were 

taken  from  reference  7. 
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Table  9.      Estimates  of  1974  Energy  Use  in  the  Commercial    Sectors11 


Energy 
Use 


End 
Use 


Gas  Coal  Electricity  Oil  LPG 

(1012  Btu)        (1012   Btu)        (1012  Btu)        (1Q12   Btu)        (1Q12  Btu) 


Space  Heating 

11.8 

1 

8 

0.1 

3.8 

2.5 

Hater  Heating 

2.2 

0 

0.1 

0.1 

0.2 

Refrigeration 

0 

0 

2.6 

0 

0 

Lighting 

0 

0 

2.9 

0 

0 

Asphalt  and 

Road  Oils 

0.1 

0 

0 

3.7 

0.2 

Other 

0.2 

0 

8 

0.5 

2.7 

0.2 

Total 

14.3 

2. 

6 

6.2 

10.3 

3.1 

*     These  estimates  are  crude  and  should  be  used  with  caution. 
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Table  10.  Baseline  Commercial  Energy  Coefficients (Energy  Use  per  Unit  of  Output) 


■ 


Energy 
Coefficient 
End 
Use 

Gas 
ft3/$ 

(Btu/$) 

Coal 
lb/$ 

(Btu/$) 

Electricity 
KWH/$ 

(Btu/$) 

Oil 

gal/$ 
(Btu/$) 

LPG 
gal/$ 

(3tu/$) 

Space  Heating 

5.58 
(5,582) 

0.106 
(851) 

0.014 
(47) 

0.0129 
(1,798) 

0.0129 
(1,183) 

Water  Heating 

1.09 
(1,041) 

0 

0.014 
(47) 

0.0003 
(47) 

0.0010 
(95) 

Refrigeration 

0 

0 

0.360 
(1,230) 

0 

0 

Lighting 

0 

0 

0.402 
(1,372) 

0 

0 

Asphalt  and 
Road  Oils 

0.05 
(47) 

0 

0 

0.0126 
(1,750) 

0.0010 
(95) 

Other 

0.95 
(95) 

0.047 
(378) 

0.069 
(237) 

0.0092 
(1,277) 

0.0010 
(95) 

Total 

6.76 
(6,764) 

0.154 
(1,230) 

0.859 
(2,932) 

0.0351 
(4,872) 

0.0159 
(1,466) 
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Agricultural  Demand 

There  is  little  data  reported  for  agricultural  energy  usage.  It  is 
very  difficult  to  separate  energy  use  for  agriculture  from  that  used  for 
rural  residential.  An  analysis  for  each  fuel  is  made  below. 

Petroleum:  Reference  1  reports  that  in  1969  $3?. 976  x  10^  were  spent 
on  petroleum  fuel  and  oil  for  agricultural  production.   If  we  assume  that 
the  average  price  was  25£/gal  in  1969  this  indicates  3,090  x  10^  bbl  of 
petroleum  was  used  in  agriculture  in  1969. 

Electricity:  No  direct  data  for  agricultural  energy  use  was  available, 
It  is  likely  that  most  of  it  is  reported  in  the  residential  category.  When 
residential  calculations  were  made,  2,960  KWH/yr/household  were  left  unac- 
counted for.  Some  of  this  will  be  used  in  residential  lighting  and  in 
small  appliances,  etc.  However,  a  large  part  of  it  may  represent  agri- 
cultural usage. 

First,  consider  what  would  happen  if  all  of  the  residual  was  assumed 
to  represent  agricultural  usage.  From  reference  2  there  were  217,300 
households  in  1970.  This  gives  643  x  106  KWH/yr  in  the  residual  category. 
Reference  1  reports  the  agricultural  output  in  1969  as  $576  x  10  .  Com- 
paring these  two  gives: 

1.12  KWH/$  Output 
For  comparison  purposes,  data  were  taken  from  reference  3  for  Cali- 
fornia. Using  these  data  the  following  results  were  obtained. 
Agriculture  Energy  Use  in  California 
Year  KWH/$  Output 


1954 
1959 
1964 


1.87 
1.83 
1.54 
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It  is  unclear  whether  or  not  rural  residential  use  is  included  with 
this.  It  is  likely  that  the  values  above  are  somewhat  high.  Also,  agri- 
culture in  California  is  much  different  than  in  Montana.  The  more  inten- 
sive farming  practices  in  California  would  probably  lead  to  higher  energy 
usage.  However,  the  above  results  do  justify  assuming  that  a  significant 
part  of  the  residual  comes  from  agricultural  usage.  For  initial  purposes 
let  us  assume  that  all  agricultural  electricity  is  sold  under  residential 
classification  and  that  it  is  75%  of  the  residual.  This  gives: 
Electric  Energy  Use  =  0.84  KWH/$  Output 

Natural  Gas:  There  is  essentially  no  data  available  for  agricultural 
use  of  natural  gas.  According  to  reference  4  a  significant  part  of  the 
total  comes  from  gas  siphoned  off  of  pipelines  crossing  farms  in  exchange 
for  right-of-way.  This  is  not  metered  or  paid  for.   It  probably  shows  up 
in  our  totals  as  transmission  losses.  Due  to  the  cost  of  installing  dis- 
tribution facilities,  farm  use  is  probably  limited  to  larger  users  and 
shows  up  as  "commercial."  For  the  present  we  will  assume  no  significant 
agricultural  use  of  natural  gas  in  Montana  until  more  data  are  available. 

Coal:  There  is  no  significant  use  of  coal  in  agricultural  production. 

Bottled  Gas:  A  significant  fraction  of  total  LPG  is  probably  used 
in  agriculture.  However,  at  the  present  time  no  good  data  are  available. 
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Energy  Supply 

The  energy  supply  part  of  the  model  is  relatively  simple  compared  to 
the  energy  demand.  The  actual  supply  modeling  is  done  in  the  economic  part 
of  the  model.  It  gives  the  total  production  of  gas,  oil,  and  coal.  How- 
ever, most  of  the  coal  production  is  destined  for  export  and  a  significant 
part  of  the  oil  and  gas  production  is  not  available  for  Montana  use. 


A-  81 


m 


Energy  Supply  References 

1.  Montana  University  Coal  Demand  Study  Team,  Projections  of  Northern 
Great  Plains  Coal  Mining  and  Energy  Conversion  Development  1975-2000, 
May  1976. 

2.  Montana  Environmental  Quality  Council,  Montana  Energy  Policy  Study, 
February  1975. 

3.  Montana  Energy  Advisory  Council,  Montana  Historical' Energy  Statistics, 
September  1976. 


H 


A-82 


& 


Section  A-IX,  Residential  Energy  Conservation  Sub-Model 

Space  heating  dominates  other  uses  of  residential  energy  in  Montana. 
Most  of  this  is  done  with  natural  gas.  For  these  reasons,  the  investi- 
gation of  conservation  in  the  residential  sector  focused  on  space  heating 
requi  rements . 

Space  heating  requirements  for  a  typical  residence  may  be  reduced  by 
making  investments  in  insulation  and  devices  which  reduce  the  heat  loss 
or  by  changing  the  life  style  of  occupants  so  that  less  energy  is  needed. 
The  first  option  tends  to  be  an  economic  decision  and  can  be  analyzed  in  a 
straightforward  manner.  The  liklihood  of  changed  life  styles  and  the 
effect  are  more  difficult  to  assess.  The  potential  energy  conservation 
from  both  investment  and  life  style  change  was  determined  and  a  relatively 
simple  dynamic  model  (which  serves  as  a  sub-model  to  the  larger  Montana 
model)  was  constructed  to  describe  the  process  whereby  the  conservation 
takes  place. 

Energy  Conservation  With  Insulation 

The  potential  energy  savings  from  insulating  were  determined  by  analy- 
zing a  "typical"  house  for  various  levels  of  insulation.  The  house  used 
for  the  analysis  is  described  in  Figure  1  and  Table  1.  The  insulation 
options  investigated  are  listed  in  Table  2.  These  options  represent  a 
very  wide  range  of  possibilities  and  include  most  of  the  cases  likely  to 
be  encountered  in  the  near  future. 


Cost  estimates  to  install  the  different  forms  of  insulation  for  both 
new  and  existing  houses  were  made.  These  estimates  were  based  on  published 
literature,  local  retail  prices,  and  discussions  with  local  dealers.  The 
resulting  estimates  are  shown  in  Table  3. 
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The  energy  savings  realized  by  installing  this  insulation  were  deter- 
mined using  standard  heat  transfer  calculations.  These  calculations  are 
somewhat  tedious  and  are   presented  separately  in  this  section.  The  resulting 
decreases  in  heating  requirements  are  shown  in  Table  4.  The  actual  energy 
saving  will  be  significantly  greater  than  the  Btu  savings  shown  as  they  do 
not  include  the  furnace  efficiency. 

Energy  Conservation  by  Life  Style  Change 

In  addition  to  the  savings  from  insulation,  significant  savings  in 
space  heating  requirements  can  be  achieved  by  actions  taken  by  the  occu- 
pants. Setting  the  thermostat  back  by  6°  reduces  energy  use  by  13%  in  a 
northern  climate."  A  10°  setback  from  10  p.m.  to  6  a.m.  can  result  in  a 
12%  savings.  Using  standard  heat  load  calculations  shown  later  in  this 
section,  it  was  determined  that  closing  off  an  unused  room  can  save  6%. 

These  are  typical  of  the  savings  which  can  be  achieved  with  relatively 
simple  measures.  More  drastic  life  style  adjustments  would  be  required 
for  larger  savings.  It  is  difficult  to  predict  to  what  extent  these 
energy  conserving  measures  will  be  employed.  However,  it  is  pretty  cer- 
tain that  the  less  drastic  measures  would  be  taken  if  heating  cost  became 
a  significant  fraction  of  total  income. 

Residential  Energy  Conservation 

The  effects  of  conservation  by  insulating  are  simulated  by  dividing 
all  housing  into  three  categories: 

Class  A  -  No  insulation 

Class  B  -  Marginal  insulation  consisting  of  6"  of  ceiling  insulation 

Class  C  -  Full  insulation  consisting  of  12"  of  ceiling  insulation, 

4"  of  wall  insulation  4"  of  floor  insulation,  storm  windows, 
and  storm  doors. 
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The  model  keeps  track  of  how  many  houses  are  in  each  class  and  based  the 
space  heating  energy  requirements  on  these  numbers.  The  number  of  houses 
in  each  class  may  change  by  new  houses  being  added,  existing  houses  being 
abandoned  and  existing  houses  being  upgraded.  The  structure  of  this  model 
is  shown  in  Figure  2. 

New  houses  are  divided  among  the  three  classes  according  to  the  eco- 
nomic gain  expected  from  adding  insulation  when  building.  The  return  is 
calculated  in  the  form  of  a  payback  period.  A  function,  shown  in  Figure 
3  is  used  to  determine  the  fraction  which  will  be  in  each  class.  This 
function  is  used  to  calculate  the  fraction  which  will  be  class  C  and  the 
fraction  which  will  be  class  A.  The  remainder  are  then  class  B. 

Existing  houses  are  upgraded  according  to  the  payback  period  for 
adding  insulation.  Again,  a  function  is  used  to  determine  the  rate  at 
which  insulation  is  upgraded.  The  function  is  shown  in  Figure  3.  Class 
A  houses  may  be  upgraded  to  either  class  B  or  class  C.  Class  B  houses 
may  be  upgraded  to  class  C.  The  rate  at  which  upgrading  takes  place  is 
determined  for  all  cases  by  the  same  function. 

The  rate  at  which  existing  houses  are  abandoned  is  determined  by  the 
average  life  span  of  a  house. 

The  actual  energy  for  space  heating  for  each  class  of  house  is  then 
calculated.  The  structure  of  this  part  of  the  model  is  shown  in  Figure  4. 
The  energy  per  household  is  calculated  by  modifying  the  "baseline"  energy 
use,  which  corresponds  to  no  insulation,  by  the  savings  which  results  from 
insulating.   In  addition,  there  is  a  savings  from  life  study  changes.  These 


- 
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are  determined  by  the  functions  shown  in  Figure  5.  Different  functions  are 
used  for  different  classes  of  housing,  since  there  is  less  to  gain  by  taking 
these  measures  in  houses  which  have  high  quality  insulation  and  thus  low 
energy  requirements  to  begin  with.  Also,  people  with  lower  incomes  will 
tend  to  be  more  affected  by  higher  prices  and  they  tend  to  have  the  houses 
with  poorer  insulation. 

In  order  to  use  the  conservation  model,  the  number  of  housing  units 

in  each  insulation  class  must  be  specified.  No  reliable  data  is  available 

to  determine  these  numbers.  Strictly  on  the  basis  of  judgment,  the  total 

number  of  households  in  Montana  were  divided  as  follows: 

30%  class  A,  40%  class  B,  30%  class  C. 

The  main  inputs  to  the  model  are:  the  price  of  energy,  the  number  of 
new  households,  and  the  "baseline"  energy  requirements  of  an  uninsulated 
house.  The  price  is  a  scenario  input  and  the  number  of  new  households 
comes  from  other  parts  of  the  Montana  model . 

A  calculation  of  heat  losses  from  our  "typical"  home  (if  uninsulated) 
shows  that  it  requires  about  250  x  10^  Btu/yr  for  space  heating.   If  we 
assume  a  furnace  efficiency  of  55%,  it  would  require  about  450  x  10° 
Btu/yr  of  gas  to  heat  this  home  in  Montana.  This  is  considerably  higher 
than  the  gas  demand  for  residential  space  heating  used  in  the  gas  demand 
section  of  this  report.  This  discrepancy  could  result  from  several 
causes.  The  size  and  structure  of  the  "typical  house"  may  not  be  typical. 
A  significant  number  of  households  live  in  multiple  family  units  which  use 
less  energy  for  space  heating.  Also,  many  people  may  supplement  their 
gas  heating  system  with  wood  or  electricity.  There  is  not  enough  data  to 
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determine  for  sure  how  important  these  factors  might  be.  To  account  for  the 
discrepancies,  the  model  is  scaled  to  give  space  heating  energy  requirements 
consistent  with  the  earlier  estimates  in  this  report  (see  section  on  resi- 
dential gas  demand). 


Heat  Transfer  Calculations 

The  characteristics  of  the  typical  house  are  described  in  Figure  1 
and  Table  1.  The  thermal  resistance  for  the  walls,  floors,  and  roof  were 
calculated  from  data  in  the  "ASHRAE  GUIDE  AND  DATA  BOOK."  The  component 
thermal  resistances  are  shown  below: 


Wall 

Roof 

Floor 

Component 

R-value 

Component 
Outside  Air 

R-value 
0.17 

Component 

R-value 

Outside  Air 

0.17 

Outside  Air 

0.17 

Siding 

0.81 

Shingles 

0.43 

Foundation* 

5.26 

Plywood 

0.63 

Plywood 

0.61 

Air  Space 

1.00 

Air  Space 

1.17 

Air  Space 

1.00 

Plywood 

0.78 

Plaster  Board 

0.45 

Plaster  Board 

0.45 

Flooring 

0.57 

Inside  Air 

0.68 

Inside  Air 

0.68 

Inside  Air 

0.61 

TOTAL 


3.91 


TOTAL 


3.33 


TOTAL 


8.37 


Windows  were  assumed  to  be  single  pane  and  undraped  which  gives  an 
R-value  of  0.88.  The  doors  consisted  of  1"  of  wood  giving  an  R-value  of 
1.25. 


The  total  heat  transfer  areas  of  each  type  of  surface  are  divided  by 
their  respective  R-values  to  determine  their  conductance. 
Surface  Area  R-value 


Wall 

Roof 

Floor 

Windows 

Doors 


680 
1200 
1200 

300 
63 


3.9 

3.3 
8.4 
0.88 
1.25 


TOTAL, 


Conductance 

174 
363 
143 
341 
50 

.   1072 


*  Foundation  assumed  to  have  little   infiltration. 
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Thus  the  total  conduction  heat  loss  is  approximately  1070  Btu/hr/degree 
of  temperature  difference. 

In  addition  to  the  conduction  losses,  there  are  losses  from  infiltra- 
tion. This  was  evaluated  on  the  standard  basis  of  1  1/2  room  air  changes 
per  hour.  This  corresponds  to  a  room  with  doors  or  windows  on  two  sides. 
The  total  heat  loss  is  determined  by: 


Q  =  CxpxVxA 


where: 


Q   is   the  heat  loss   per  unit  of  time  per  degree  of  temperature 
difference 

Cp   is   the  specific  heat  of  air 

p  is   the  density  of  air 

V   is   the  volume  of  the  house 

A   is   the  air  change   rate. 

For  our  house  we  have 

Q  =  0.24  Btu/lb/F  x  0.77   lb/ft3   x  7   ft  x  1200   ft3   x   1.5/hr 

=   233  Btu/hr/Degree  Temperature  Difference 

The  total    heat  loss   for  the  house  is   the  combination  of  the   infiltration 
and  conduction   losses  or  1300  Btu/hr/Deg  AT.     This  must  be  multiplied 
by  24  and  by  the  number  of  heating  degree  days   per  year  to  determine  the 
annual   heating  requirements. 

This  gives  us  our  baseline  estimate  for  a   house  with   no   insulation. 
To  calculate  the  effect  of  insulation,   its  R-value  is  added  to  the  existing 
R-value  of  the  surface.      Typical    insulation  has   an  R-value  of  3.7   per  inch 
of  thickness.     Also,    the  addition  of  storm  windows  and  doors  must  be  con- 
sidered.     The  extra  pane  of  glass   increases   the  window  R-value   to   1.82. 
The  storm  door  increases   the   total   door  R-value   to  2.25. 
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The  storm  windows  decrease  the  total  infiltration  by  1/3  and  the  storm 
doors  by  10/i.  The  following  table  summarizes  the  changes  in  the  overall 
heat  loss  for  each  insulation  option  in  Table  2. 


1 

Conduction 

Infilt 

ration 

Total 
Heat 

Wall 

Roof 

Floor 

Windows 

Doors 
50 

Windows 
233 

Loss 

None 

174 

364 

143 

341 

1300 

3"  ceiling 

174 

83 

143 

341 

50 

233 

1025 

6"  ceiling 

174 

47 

143 

341 

50 

233 

990 

6"  ceiling,  4"  walls 

36 

47 

143 

341 

50 

233 

850 

6"  ceil ing,  4"  wal  1  , 

4"  floor 

36 

47 

52 

341 

50 

233 

760 

6"  ceiling,  4"  wall  , 

4"  floor,  storm 

36 

47 

52 

165 

28 

133 

460 

windows  &  doors 

12"  ceiling,  4"  wall, 

36 

25 

52 

165 

28 

133 

440 

4"  floor,  storm 

windows  &  doors 

12"  ceiling,  6" 

26 

25 

39 

165 

28 

133 

415 

wall  ,  6"  floor, 

storm  windows  & 

doors 

It  was  desired  to  compare  masonary  walls  with  the  frame  wall  analyzed 
to  see  if  they  have  comparable  thermal  resistances.  The  typical  brick 
wall  below  is  evaluated  for  this  purpose. 

Common  Air  Space  Component  R- value 

ace  Brick   Brick   \    ,.  , , 
v        .    \    Wal  1 

>  1  /  Board 


Outside  Air 
Face  Brick  (4") 
Common  Brick  (4") 
Air  Space 
Wall  Board  (1/2") 
Inside  Air 


TOTAL , 


0.17 
0.44 
0.80 
0.97 
0.45 
0.68 

3.51 


It  was  also  desired  to  evaluate  the  effect  of  not  heating  an  unused 
room.  By  doing  this  the  interior  walls  of  the  room  add  a  thermal  resis- 
tance. The  room  analyzed  is  shown  on  the  next  page. 
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■ 


■ 


3'    x  5'   window 


/ 


10' 


12' 


( 


Door 


>y 


3'    x  5' 
window 


If  the  room  is  heated  the  conductance  for  that  section  of  the  house   is: 

A  AAA 

..     _     outside  wal  1  ,     window   ,     ceil  ing   ,     floor 

UH  "           R  R                  R                    R 

wall  window         ceiling         floor 

A  =    (12'    x  10')   *  7'    -   A  =   124   ft3 

ou  tside  wall  window 

A  =30  ft2 

window 

A       =  12'  x  10'  -  120  ft2 
ceil ing 

A     =  12'  x  10'  =  120  ft2 
floor 


A- 9  0 


Using  the  R-values  for  an  uninsulated  house  gives 

ii     _  124   ,      30         120   .    120       117  D.    ,,     /n  „    ,T 
UH   -   3^  +  08  +   373  +  O  "   117  Btu/h'VDeg  AT 


With   the  room  unheated,    the  effective   conductance  is   determined  by 
combining  the  conductance  of  the  outside  surface  with  that  of  the  interior 
wall    by: 


U, 


1 


1/U         +  i/u 
i .  w .  o .  s . 


The  interior  wall  is  a  2"  x  4"  frame  wall  on  16"  centers  with  1/2' 
wall  board  on  either  side.  The  overall  resistance  is  calculated  below 


Component 

R-value 

Still  Air 

0.68 

Plaster  Board 

0.45 

Air  Space 

1.17 

Plaster  Board 

0.45 

Still  Air 

0.68 

TOTAL 


3.43 


The  door  is   1"  of  wood  and  has  an   R-value  of  1.25.     The  overall    con- 
ductance of  the  room  is   then: 

A         A 

JL  +  -i 

R         R 


Ui.w. 


W  +  _d  _   154-21    +    21 


3.3         1.25 


55.6 


w         d 


With  the  door  closed  the  total  conductance  is 


37.7  Btu/hr/Deg  AT 


1   +  1 


55.6   117 

The  reduction  in  the  conductance  with  the  room  closed  is  then 
117  -  37.7  =  79  Btu/hr/Deg  AT 
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Compared  to  the  total  conductance  of  1300  for  the  house  initially  this 
represents  a  6%   decrease  in  heat  loss  due  to  lower  conduction  loss.  No 
decrease  in  infiltration  is  considered. 


m 


m 
- 
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1.      Coj  I     >n   Features  of  Typical  House 

Nominal  Dimensions   -    30'  x  40',     7  ft.    ceilings 
Windows   -    20,    wo  >d<  n    i  u     le  pane,    average   size  3'  x  5',    not 
weather   stripped 
,rs   -   3,    1"   '        '  d,    not  we         er  stri]  .    average  size 

3'  x  7' 

Table  2.      Insulation  Considered 

1 .  none 

2.  3"  ceiling 

3.  6"  ceiling 

4.  6"   ceiling,    4"  walls 

5.  6"  ceiling,    4"  walls,    4"  floor 

6.  6"  ceiling,    4"  walls,    4"  floor,    storm  windows,    storm 
doors 

7.  12"  ceiling,    4"  walls,    4"  floor,    storm  windows,    storm 
doors 

S.        12"  ceiling,    6"  walls,    6"  floor,    storm  windows,    storm 
doors 
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3.      Cost  Estimates  for  Insulati 


a)     CcjsI  for   units   of  insulation  insta!: 


I  ion 
3M  attic 

6"  attic 
12"  attic 
4"  wall 
6"  wall  * 
4"  fioor 
6:!  floor 
Storm  w  i  n  d  o  w  s 
Storrn  doors 


[ .  i  ..-  Construction 
20c/ft2 
25£/ft2 
35^/ft2 
25c/ft3 
25^/ft2 
259/ft2 
30^ /ft2 
$3  5  /window 
$120/door 


■  ;g  Structures 
25^/ft2 
30£/ft2 
40^/ft2 
50c/ft2 


35d/ft2 
40cl/ft2 
$40 /window 

$130/door 


*    6"    wall   insulation  requires  2  x  6   framing   and   is   thus  possible  o.nl\    on 
houses  which  already   have  this   or   new  construction.      Using;  2x6 
framing  on  24"   centers   actually    results    in  a.  cheaper,     stronger 
wall. 


1. 

0 

2. 

$240 

3. 

$300 

4. 

$490 

5. 

$790 

b)     Cost  for   insulation  options 
New  Old 

0 

$300 

$360 

$740 

$1  160 
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Talil  e   3.      Co;.:  i  nu  ed 


b)     Cost  for  insulation  options 


New 

6.  $1850 

7.  $1970 

8.  $2030 


Old 
$2350 
$24  7  0 


Table  4.      Energy  Savings  with  Insulation  for  Typical  House 
Insulation  Level  %  Savings  BTU/YR  Savings* 

1  0  0 


2 
3 
4 
5 
6 
7 
8 


21 

26 

36 

44 
66    1/2 
68   1/2 
70   1/2 


53  x   10'3  BTU/j  r 

66  x  10s  Bl  L'/yr 

91  x   10s  BTU/yr 

1  12  x  106  BTU/yr 

169  x  106  BTU/yr 

174  x  10*  BTU/yr 

179  x   10"  BTU/yr 


■■'■      Ba  ^c-d  on  8000  heating  degree  days  per  year.      Actual    heating 

L3  ement,    not  furnace  fuel   use.      Heating  requirements  with 
no  insulation  are  253  x  106    BTU/yr. 
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1/2"  Plywood 

1/2"  x  8" 
Drop  Siding  **  Ik 


2"  x  4"  Stud  Framing 
16"  Centers 


1/2"  Plaster  Board 


a)  wall 


Asphalt  Shingles 

1/2"   Plywood 

2"  x  6"  Rafters 
16"  Centers 


2"  x  6"  ceil ing 
Joist 


1/2"  Plaster  Board 


House  on  cement  block  foundation, 

no  basement.  Average  effective 

foundation  exposure  is  2  ft  c/om  u  a       a 
\i  I  i  i,  i  i  i  i  i  ii  :  i  i  i  i  i  i  i  J*  5/8  Hardwood 

y/Avwvy//WiAV^L  5/8«  Plywood 


2"   x  8"   Floor  Joist 
16"  centers 


c)    floor 


Figure  1 .--Construction  of  Typical  House 
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JUL 


Class    C 
Housing 


B    Housing 
Upgraded 


Life  Span 


New    Households 


Class    C 
Housing 


Life  Span 


C±    A    Housing  Upgraded 
1o     C 


B    Housing 
Upgraded    to  C 


Class    B 
Housing 


Ld? 


New   Households 


_<-^ 


K 


A     Housing 
Upgraded   to  B 


7T 


Class    A 
Housing 


STRUCTURE   OF 
HOME    INSULATION    MODEL 


Cost  for   B 

in    Existing  Savings  for  B 


Cost    for  C 

in  Existing         Savings  for  C 


Baseline 
E  nergy  Use 


Cost  for  C  Savings  for  C       Cost  for  B 

in  New  in   New 


4 


Energy    Price 


Baseline   Energy  Use 


Savings  for  B 


Figure  2 
A-97 


o 
1 


Fraction  of  existing  houses  upgraded  per  year 
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Reduction    in    Heating  Energy  Use  by  Lifestyle  Changes 


Fuel  Pries 


Baseline 
Energy  Use 


Savings   with 
Class    B 


Savings  with 
Class    C 


Figure  4 . --Calculation  of  Space  Heating  Fnergy  Requirements' 
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%   Reduction  in  Energy  Use  for  Space  Heating 


o 


o 


o 

CTv 


CO 
r+ 


Figure  5. --Energy  Savings   from  Lifestyle  Changes 
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APPENDIX  B 
Computer  Program  Listing 


■ 


■ 

■ 


a 


c 

c 

c 


COMMON  FOR  INPUT/OUTPUT  VARIABLES 


COM 

2  A3 

3  A I 
COM 
COM 
COx 
COM 
COM 
COM 
COM 
COM 
COM 
COM 
COM 
COM 
COM 
COM 


Mi  O.N 

0» 

1(10 

MON 

VI  O.N 

MON 

MON 

MON 

MON 

MON 

MON 

MON 

MON 

MON 

MON 

MON 

MON 


IRLI/A4>A5  »A6  »A9  »A12  *A14  »A15  »A17» A20>A26  *A27  >  A.-:  3  »A2 


» h 


! 

/AG  I 

/  F  0  R 

/GAS 

/COA 

/  M  I  N 

/SDE 

/OEM 

/OUT 

/L43 

/ca:j 

/ENR 

/LAN 
/HHT 


2 ( 25 ) ♦ A I  3 ( 25 ) » A 1 4 ( 3  5 ) ♦ A41 » A42 

/  A 1  »  A  2  »  A  3 

ST  I  /  A7*A3 

01/  A10»All-*A13 

LMI/  A15»A40 

EI/  A 1  8  ,  A  1  9 

MI/  A2l*A22»A23»A24»A25 

NDI/  A43 

PTI/  A39 

CRI/A33 

INI/  A  I  5  !  3  3  ) 

DMI/  AI5{ 5  ) 

031/  A31*A32 

01/    A33»A34>A36»A37 

CNI/  A44,A45,A4S 


c 

c 
c 


COMMON  FOR  PARAMETER  VARIABLES 


COM 
COM 

COM 
COM 
COM 

COM 
COM 
COM 

2  P 
COM 
COM 
CO 
COM 
COM 

2  P  A  8 
COM 

COM 


MON 
MON 
ION 

•10N 

-I  ON 

MON 

MON 

MON 

A63( 

MON 

MOM 

MON 

MON 

MON 

K25 

MON 

MON 


/GE.NRLP 
/AGP/  P 
/FORST? 
/GASO?/ 
/COALMrj 
/MINEP/ 
/SOEMP/ 
/DEMNDP 
25) |PA6 
/0U7PTP 
/LA30RP 
/CAPI\y 
/ENRDMP 
/DEMOGP 
) »PA82( 
/LAND?/ 
/HHTCNP 


/  IP 
1  >P2 
/  PI 
PI  8 
/  P2 
P30» 
P34» 
/  P3 
91  25 
/  IP 
/  P4 
/P51 
/  IP 
/P54 
40  )  » 
P62. 
/  PI 


43 
»P 

4* 

»P 

6  » 

IP 

P3 

7> 

) 

44 

5> 


»IP50»DT 

3»P4»P5»P6*P7»P3»P9»PlO>Pll»Pi; 

P15»P15»P17 

19 »P20 »P21 »P2  2  »D23  >P24>P25»P1 3' 
IP27*IP28*IP29*°93»P94, IP96»Plv 

31  »IP32  ♦  IP33»IP97 
5>P36 

P38  ♦  P39  » PI 02  » ? 40  »P41  ♦  I ?42  »PA66 


(23)  » 


,j:o 


c 
c 
c 


COMMON  FOR  STATE  VARIABLES 


b3»D64,P35»D55 
b»Pl07»P103»Pl09,DllO»Plll»Pll2 


PI 


COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
?.  S13(6 
COMMON 
COMMON 


/AGS/S1 
/FORSTS 
/GASOS/ 
/COALMS 
/MINES/ 
/S0EM5/ 
/OUTPTS 
/LABORS 
/CAP  INS 
/DEMOGS 
0  )  »S14( 
/LANDS/ 
/HHTCNS 


5, SI 

/S17 

S19* 

/S23 

S24, 

525 

/  S3 

/  S2 

/S3^ 

/  S7 

30) 

S29 
/  S4 


S20»S21*52  2»S1(50)»S2(50)»S3(20),S4(20) 

♦  S37»S5!20) 

S6(20) 


6»527»S23  >S39 

(35  )  »S33 (35) 

(6»25) »S3(60) »S9(60J  »S10(60) >S11(50) » 512(60) » 


»S30»S31»S32iS33»S36(3J 
0»S41»S42»S43 


B-l 


■ 

■ 


C  COMMON  FOR    AFNC 

COMMON  /AFNCP/  XFNC(75»30»2) til  .  ^  C  7 3  »  3  ) 

ICRX=5 

C  READ  SIMULATi  '.  CONTROL  VARIABLES 

READl ICRX»101)  TSTRT»TE.\D»DT»ISOiR»ISCON»I    :  , I  SPRM ,  I    r»i 
1C1  FORMATC3E12. 0,616) 
C  READ  AG  PARAMETERS 
REAO( ICRX ,100) 
LOO  FORMAT 16E12.0) 

READ  I  ICRX  ,100)  P7»P8»P9»P10»P11»P12 
READ  (ICRX  ,100)  P13 
C  READ  FOREST  PARAMETERS 

READ  (ICRX»100)  P14 »P15 »P16 »P1 7 
C  READ  GASOIL  PARAMETERS 

READ  (ICRX, 100)  Rl 3 »P19»P20 »P21 »P22 »P23 
READ  (ICRXtlOOJ  P24,P25,P104 
C  READ  COALMN  PARAMETERS 

READ  (ICRX,100)  P2 6 *P93 » P94 »P 1 03 
READ  (ICRX»2C0)  I P27 t IP28 ♦ I P29 , I P96 
200  F0RMAT(12I6) 
C  READ  MINING  PARAMETERS 
READ  (ICRXilOO)  P30 

READ  (ICRX, 200)  IP31 » IP32 » IP33 t IP97 
READ  SUPDEM  PARAMETERS 

READ  (ICRX»100)  P34»P3p,P36 
READ  DEMAND  PARAMETERS 

( ICRX  ,  ICO)     P3  7,P38,D39»P4  0,P4i ,P1C< 
(ICRX»200)     IP42»IP43 


C 
C 


■J  c   *  r» 

i ;.  E  A  D 

READ 

READ 

READ 

READ 

(PA66( J) »J=1 ,25 ) 
(PA67( J) >J=1»25) 
(PA68U) »J  =  1 ,25) 
(PA69( J) »J=1»25) 


( ICRX, 100) 
(ICRX  ,100) 
(  I  C  R  X  »  1 0  0  ) 
(  ICRX  ,100) 
READ  OUTPUT  PARAMETERS 

READ  !  ICRX  ,200)  IP44 
CO  10  1=1,25 
10  READ( ICRX ,100)  CPA70 ( I » J ) » J=l » 10 ) 

DO  20  1=1,25 
20  READl  ICRX  ,100)  ( PA71 ( I  * J > » J  =  l ,25 ) 
READ  LA33R  PARAMETERS 

READ  (IC?X»100)  P45  »P46»P47  >P43  »P49 
READ  (ICRX, 200)  IP50 
READ  (ICRX, 100)  ( p A72 ( J )  ,  J= 1  ,  3 5 ) 
READ  (ICRX,  200)  ( I PA73 ( J ) ♦ J=l > 35 ) 
READ  CAPINV  PARAMETERS 
READ  (ICRX  ,100)  P5  1 
READ  (ICRX, 100)  ( PA74 ( J )  , J=l ,35 ) 
READ  (ICRX,  100)  I PA75 ( J ) , J=l , 35 ) 
READ  tICRX»200)  (  I PA76 ( J )  » J  =  1  ,  3 5 ) 
READ  (ICRX,  100)  ( PA90 ( J ) » J=l ,25 1 
READ  EXRDMD  PARAMETERS 

READ  (ICRX  ,200)  IP52,IP53 

READ  (ICRX,  100)  ( { PA77 ( I » J ) , J  =  l , 10 > » i  =  1 , p ) 

READ  (ICRX, 100)  ( (PA73( I ,J) »J=1,10) ,1  =  1,5) 


(  (  PA79  (  I  ,  J)  *. 


REA 
READ  DEMOG  PARAMETERS 

READ  (ICRX, 100)  P54»P55 


3  5  ) 


=  i 


1  »  5  ) 


B-2 


■ 
■ 


(  i  c R  >; » .   I 

READ 

v  .  ■'.  X  1 1 0  0 ) 

R  E  A  D 

(ICRXj 

READ 

1 1 .  ;x»  ioo) 

•  EAO 

(ICRX* 100) 

READ 

( ICRX »100) 

■ 

( ICRX » 100) 

R  EAO 

!  ICRX  ,  ' 

READ 

! ICRX ♦200) 

READ< ICRX ,100) 


IP56»IP57*IP58»I: 

P  6  0 

I  :>  6  1 

(  PA801 J  )  »J=1  »4  ) 
CPA31 ( J) »J=1»25) 
(PA-;2  (J)  »J=1  »40) 
( PA83 ! J)  » J=l  .30) 
(  (  RASA!  I  » J)  »  J=l  >45  )  »  I 
I  ?  9 1  *  I P 
^95 


1»6) 


R  EAO  L  A  ND  PAP  AV  E  T  E  -:  5 


■ »  P  8  5 


'OP 


READ  (ICRX*  100) 
READ  HHTCO.N  PARAMETERS 

READ ( ICRX » 100)  P106»P107»P108»P109»P110»P111 

READ( ICRX > 100)  P112»P113 
READ  OTHER  PARAMETERS 

READ  (ICRX, 100)  P36 »P87 »P38 *P39 
READ  AG  STATE  VARIABLES 

READ  (ICRX»100)  S15»S16 
READ  FOREST  STATE  VARIABLES 

READ  (ICRX»100)  517,518 
READ  GASOIL  STATE  VARIABLES 

READ  (ICRX»100)  S19 *S20 »S21 *S22 

READ  (ICRX, 100)  ( S 1 ( J  )  » J  =  1  .  5C ) 

READ  (ICRX»100)  (S2( J)  ,J  =  1  ,50) 
READ  COAL.VN  STATE  VARIABLES 

READ  (ICRX, 100)  S23»537 

READ  (ICRX, 100)  ( S 5 ( J  )  * J= 1  ,  2 0 ) 
READ  MINING  STATE  VARIABLES 

READ  (ICRX  ,100)  S24 

READ  (ICRX,  100)  ( 5  6 ( J )  ,  J  =  i  ,  2  0 ) 
READ  SUPDEM  STATE  VARIABLES 

READ  (ICRX,  100)  S2  5 
READ  OUTPUT  STATE  VARIABLES* 

READ  (ICRX,  100)  S3B 
READ  LABOR  STATE  VARIABLES 

READ  (ICRX, 100)  S26»527»S23»S39 
READ  CAPINV  STATE  VARIABLES 

READ  (ICRX, 100)  (S3A( J)  »J  =  1  ,35 ) 

READ  (ICRX. 100)  ( S35 ( J )  , J= 1 , 35 ) 
READ  DEMOG  STATE  VARIABLES 

READ  (IC  X,100)  { <S7( I »J) #J=1»25 ) > I 

(S5( J)  »J  =  1 ,50) 
(S9( J)  ,J  =  1  ,50) 
CSIO(J) »J=1 >60) 
(Sll (J)  *J  =  1  ,60) 
( S12 (J) »J=1 ,60) 
(S13 (J) »J=1 ,60) 
<51M J)  »J  =  1 ,30) 


1,6) 


READ 
READ 
READ 
READ 
READ 
READ 
READ 


(  ICRX  ,100) 
(  ICRX, 100) 
(  ICRX, 100) 
(  ICRX,  ICO) 
(  ICRX, 100) 
( ICRX , 100) 


(  ICRX, ICO) 

REAr>  LAND  STATE  VARIABLES 

READ  (ICRX, 100)  S29  ,  530  ,S3 1 . 532 »S3 3 
REAO  (ICRX, 100)  ( S36( J)  ,  J  =  l ,3 ) 

READ  HHTCON  STATE  VARIABLES 

READt  ICRX  ,100)  SAO  ♦  S^  1  ,542 , S43 

READ  OTHER  SCENARIO  VARIABLES 

READt  ICRX, 100)  C  A  I  3  (  I  )  §1  =  1*25 ) 
READt  ICRX,  100)  ( A  I  6 (  I  )  ♦  I  =  1  ♦  35  ) 


B-3 


lil 


?E !■      t 1    <  1 1 00 )    o  »  A2  7  »  A 30  » A40 
READ!  IC  X  »100)  VI  »V2 
:  ?EAD  FUNC7I  )N  INFO 

READ  (ICRX»200)  NrNC 
0  0  2  5  0  1-1  ,\ri\C 

READ  (ICRXi200)  ( I FKCC I  * J ) » J=l *5 J 
L=IF\C( I  »4) 

READ  (ICRX»100)  ( ( XFNC ( I » J»K) » J=X »L ) tK=l *2 ) 
?53  CONTINUE 

33  4321  1*1 #20 
S3( I )=0i 
S4( I >=3. 
DO4321K=I#50 
J=K-I+1 
X  =  J 

CALL  AFNC(X»Y#3) 
S3( I J =S3( I  )+Sl ( J) *Y 
CALL  ArNC(X»Y»7J 
^321  SM  I)=SM  t  )+S2<  J)*Y 

CPRINT  ALL  INPUTS 


■ 


I =T  5T  ^T 

J=TEND 

PRINT  3Q00»I#J#DT 
3030  FORMAT ( ' 1STARTING  YEAR  ='»I5/'  ENDING  YEAR  =»»I5/ 
1  '  TIME  STEP  ='iF5.1»'  YEARS'///! 

PRINT  3001  »  I50ER* ISCON#ISDRC» ISPRM 
3331  FORMAT ('  SCENARIOS • »5X »4I3 ) 

IFCIPRNT.EQ.O)  GO  TO  7030 

PRINT  33  10 
3013  FORMAT ('  INPUT  VALUES'/) 

;<rs=i 

KTE  =  6 

PRINT  3020»KTS#KTE#P1»P2*P3»P4»P5#P6 

K  T  S = 7 

KTE=12 

PRINT  3020»KTS*KTE»P7»P8»P9»P10»P11.P12 

<TS=13 

PRINT  3030*KTS»P13 

,<TS=1A 

■<TE  =  17 

PRINT  3020»KTS#KTE#P14»P15»P16#P17 

KTS=1£ 

KTE=23 

PRINT  3320»:<T5»<TE»P18»D19,?2  0»P21»P2.?»P23 

K  T  S  =  2  A 

KTE=25 

PRINT  3020  » KTS  »KTE  »P24  »P25 

<TS=10^ 

PRINT  3331  »:<TS»P13A 

KTS=26 

PRINT  303C»KTS.P26 

K  T  S  =  1 0  3 

PRINT  33  31  »:<T5,P133 

KTS=27 

K  T  E  =  2  9 
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PRINT  3040»KTS,KTE, IP27. IP22 ,  IP29 

KTS=96 

PRINT  3060*.<TS»IP96 

KTS=30 

PRINT  3030»KTS»P30 

KTS=31 

KTE=33 

PRINT  3040fKTS»KTE» IP31, IP32 »IP33 

KTS=97 

PRINT  3060»KTS»IP97 

KTS=34 

KTE=36 

PRINT  3020»KTS»KTE»P34»P35 »P36 

KTS=37 

KTE=41 

PRINT  3020*KTS»KTE*P37»P38»P39»P40»P41 

KTS=102 

PRINT  303l»KTS»P102 

KTS=42 

KTE=43 

PRINT  3040»KTS»KTE*IP42»IP43 

KTS=66 

PRINT  3050»KTS»tPA66(J)»J=lf25) 

:<T5  =  67 

PRINT  3050 »KTS. CPA67I J) »J=1 ,25 ) 

KTS=68 

PRINT  3050»KTS>(PA68(J) »J=1»25) 

KTS=69 

PRINT  3050»KTSf (PA69( J) »J=1»25 ) 

KTS=44 

PRINT  3060»KTS»IP44 

KTS=70 

PRINT  3050».<TS.((  PA70  {  I  »  J  J  ,  J  =  l ,  10  )  ,  I  =  1 ,  25  ) 

KTS*71 

PPINT    3050»KTS»(  (PA7KI  »J)  »J=1,25J  »I=1,25) 

KTS=45 

KTE=49 

PRINT  3020»KTStKTE*P45  »D4S  >P4  7  ,P4S  »P49 

KTS=100 

KTE=101 

PRINT  3321 »KTS»KTE,P100»P10I 

KTS=5C 

PRINT  3060»KTS»IP50 

KTS=72 

PRINT    30  5CKTS,  (  PA  72  (  J  )  ,  J  =  l  ,  3  5  ) 

KTS=73 

PRINT  3C7C»<TS»(  IPA73( J)  ,J=1,3  5> 
<TS=51 

PRINT  3030,KTS»P51 
KTS=74 

PRINT  3050tKTS»(PA74(J)fJ=l»35J 

KTS=75 

PRINT  3050.KTS»(PA75<J)»J=lt35J 

KTS=76 

PRINT    3070»KTS»{  IPA76U)  »J=1»3  5) 

KT5=90 

PRINT  3C50»KTSt (PA90CJ) »J=1»25) 
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kts=52 

KTE=53 

PRINT  3040»KTS»KTE»IP52, IPS  3 

<TS=77 

INT  3050,  <TS 
■ TS=73 

PRINT  305C»:<TS 
KTS=79 

PRINT    305CKTS 
KTS=34 
K  T  E  =  5  5 

PRINT  3C20»KTS 
KTS=56 
K  T  E  =  5  9 

PRINT    3  0  43  » MS 
KTS=60 

PRINT  3030 »KTS 
<TS=61 

PRIM    3060*.<TS 
KTS=80 

PRINT    3050».<T5 
KTS=81 

PRIM    3050  »KTS 
K  T  S  =  3  2 

PRINT    3350  >MS 
KTS=33 

PRINT  3050 »KTS 
KTS=84 

PRINT  3050»KT5 
KTS=91 
KTE=92 

PRINT  3040*:<TS 
KTS=95 

PRINT  3030 »KTS 
PRINT  3080*P62 
PRINT  3085 »P86 
KTS=106 
KTE=111 

PRINT  3321  »MS 
KTS-112   . 
KTE=113 

121  »KTS 


3035 

3033 

3  02  0 
3021 
3030 
3031 
3040 
3050 
3060 
3070 


PRINT  3 
FORMAT ( 

FCRMATl 
FORMAT ( 
FORMAT! 

FORMAT ( 

FORMAT ( 
FO  >,  -1ATI 
FORMAT ( 
FORMAT ( 
KTS=3 
PRINT  3 


P  •  »  I  2 

P  •  ,13 
P  •  »I2 

P  •  »I3 


(  (PA77 (  I  » J)  »J=1 »10)  »I=1»5) 
((PA78(I*J»»J=l»10)fI=l»5) 
({PA79<I*J)»J=1,35)»I=1,  : 

KTE»P54»P55 

KTEtIP56»IP57»IP58»IP59 

P60 

IP61 

(PASO  !  J  )  >J=1  »A-J 

(PA31 (J ) »J  =  1 >25 ) 

(PA82J J) »J=1»40) 

CPA83(J)*J=lt30) 

( (PA84  1 I *J) *J=1 ,45) » 1=1,6) 

KTE»IP91,IP92 

P95 

P63»P64»P85»P65 

P37»P88,P89»P93,?94 

<T£»DI06»P10  7,PlC3>D10  9>rjiiG,?lll 


KT£fP112*PU3 
P86,    P87»    P83,    P89,       P93»    P94  •  »T4,0  »4G12.  5/T40  »2G12  .5  ) 

P62,    P63»    P6^»    P35,    P65 • »T40 , 5G12-5 ) 
'     -    P  •  »I2 >T40»6G12.5J 
'     -    p •  ,  13 »T40,6G12»5) 
T43  ,312.3) 
TA0,6G12.5 ) 


I3'  ,12, •    -    IP'  .I2.T40.6I12) 

PA*  »I2» (T40.6G12.5) ) 

IP'  ,12  »T40» I  12) 

I 3A' » 12. ( T40»6I 12 ) ) 


80»<TS, (AI3( J)  »J  =  1  ,25) 


KTS=6 

PRINT    3180»<TS»(AI6(J) *J=1»35) 

PRINT    3120»A26,A27*A30,A40»V1»V2 
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3120  FOR  \T('  A26»A?7»A30»A^0»V1»V2'  »T^-0,6G12*5) 
3180  FORv-'AT('  AI  •  ♦  I  2  ♦  (  T40  »5G12  -  5  )  ) 

PRINT  3090 
3090  FORMAT Ml') 

<TS=15 

K  T  r  =  1 6 

PRINT  3 100  »KT5  »KTE  »S15  »S1S 

KTS=17 

KT£=18 

PRINT  3100*KTS»KTEiS17»SlS 

KTS=19 
KTF=22 
PRINT  3100»KTS»KTE»S19»S20»S21»S22 

K  T  S  =  1 

PRINT  3110»KTS> (SI ( J)  »J  =  1  »50 > 

KTS  =  2 

PRINT  3110>:<TS»(S2(J)  >  J  =  1  »  5  0  ) 

KTS*3 

PRINT  3110*KT5» (S3 ( J) . J=l »20) 

KTS  =  ^ 

PRIMT  3110.KTS*  (SMJ)  »J=1»20J 

K  T  S  =  2  3 

PRINT  3I10»'<TS»S23 

KTS=37 

PRINT  3110»;<TS>S37 

KTS=5 

PRINT  3110*KTS> (S5 ( J ) > J=l >  2  3  ) 

KTS=24 

PRINT  3110  >'<TS»S24 

KTS=5 

PRINT  3110»:<TS»  (S6(  J)  »J=1  .20) 

KTS=23 

PRINT  3110»KTS»S25 

KTS=38 

PRINT  31I0»KTS»S3S 

KTS=26 

KTF.  =  23 

PRINT  3100»KTS»KTE»S26>S27»S28 

KTS=39 

PRINT  3110tKTS»S39 

KTS=3't 


PRINT  3110».<TS»  (S3M  J)  » 


j  =  i  »  j 


5  ) 


KTS=35 

PRINT  3110»KTS»(S35(J) »J=1»35) 

KTS=29 

KTF=33 

PRINT  3100tKTS»KTE»S29»S30»S31»S32»S33 

KT5=36 

PRINT  31l0»KTS»(S36( J) #J=1»3) 

KTS=37 

PRINT  3110,:<TS>S37 

KTS=40 

KTF=43 

PRINT  3100»KTS»KTE»S40»S4l  »SA-2  »S^3 
3100  FORMAT  ('  S'tI2»'  -  S  '  »  I  2  >T<VO  »6G12  .  5  ) 
3110  FCRVATt*  S'  ,I2i  (T40»6G12.5 )  ) 

PRINT  3090 
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PRINT  3: 

PRINT  3l40»NFNC 
3140  FORMAT  (  '  NF.NC  =  '  ,  I  •'.  ) 
1  3150  I=1,NFNC 

PRINT  3  1  7 ::  »  I  >  C  IFNCI  I  ,J)  »J  =  1»5) 

L=IFNC( I >4) 

PRINT  316C»<  (XFNC1  I  * J ♦< ) > J=l *L ) »K=1 , 2 ) 
50  CONTINUE 

60  ~ON"\T(»  XFNC  •  »  (  T4C6G12  .5  )  ) 
ro  FORMAT!'  I FNC«  > I  5 » < T40 »& 1 12 ) ) 
00  CONTINUE 


31 
3  1 
2  I 
7  0 

C 

C  E 

C 


no 


TIME=TSTRT 
ICOUNT=TSTRT+.0Q001 
3  00  CONTINUE 
C  DISCRETE  PARAMETER  CHANGES 

IF(  ICOUNT.EQ.1980)  P14=P14*1»1 

IF(  ICOUNT.EQ.1990  )  P14=P14.*1»05 

IF( IC0UNT.EQ.1930 )  P86=P8&*1*05 

IF( I COUNT. EQ. 198 5 )  P86=P36*1»05 

IF( ICOUNT.EQ.1930)  P15=P15*1*03 

IF( I COUNT. EQ. 1985)  P15=P15*1.03 

CALL  AFNC(TIME»A2»20J 

CALL  ACNC(T1-E,A:0,21  ) 

CALL  AFNC(TIME»A11»22) 

CALL  AFNC( TIME »AI5 ,23) 

CALL  AFNC<  TIME»A18  >24) 

CALL  -\~.\C(  TIME.A25  ,2  5) 

CALL  «\FNC(TIME»A31«26) 

CALL  AFNC ( TIME»A32  »27) 

CALL  AFNC(TIME»A36»28J 

CALL  AFNC(TIME,A37»29J 

CALL  AFNCCTIME»A4C»32) 

CALL  AFNC(TIMEtA»33) 

IF( ISOER.EQ.l)  GO  TO  350 

PA72(23)=A 

CALL  AFNC (TINE, A, 34) 

PA72(27)  =A. 

CALL  AFNC(TIME»A»35) 

PA72I26) =A 

CALL  AFNC(TIME>A»36) 

PA72(  1  )=A 

CALL  AFNC ( TINE, A, 37) 

PA72t2  )=A 

PA72 ( 3  )  =  A 

PA72(4)=A 

PA72( 5  )=A 

CALL  AFNC (TINE,  A,  33) 

PA72(25)-A 

CALL  AFNC(TI.vE,A,3?) 

pA72l 30) =A 

CALL  AFNC < TIME »A. 40) 

PA72(29) =A 

CALL  AFNC(TIME,A»41) 

PA72( 10)=A 
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3  50 


5 


t  A  >  4  3  > 


■'-  -  ! 


IE*  A, 45) 


[  6  )  =  A 
PA72(7)=A 
PA72( 3)=A 

•  72  < 9) -A 

Cai_l  AFNC(TIMEtA«42) 

PA72( 11 )=A 

CALL  AFNCIT 

PA72( 12 )=A 

CALL  AFMC(T 

PA72J  13)  =A 

CALL  AFNC(T 

PA72(  14  J =A 

CALL  AFNC(TIME»A»46) 

PA72< 15)=A 

PA72  1 16)=A 

CALL  AFNC(TIME»A»47) 

PA72( 17)=A 

CALL  AFNC(TIME»A,48) 

PA72 (1S)=A 

PA72( 19) =A 

PA72( 20! =A 

PA72(23)=A 

CALL  AFNC(TIME*A»49) 

PA72(21 ) =A 

PA72I 22)=A 

CALL  AFNC(TIME,A»50) 

PA72( 2A)=A 

CONTINUE 

IF( ISCON.EO.O)  GO  TO  3  60 

AA5=A3 

CALL  AFNC(TIViE»A44»62  ) 

CALL  AF\C{TI"E,Z>63) 

P109=P109*Z 

P110=?110*Z 

P111=P111*Z 

P112=P112*Z 

CALL  AFNC(TI:ME»P106»64) 

CALL  HHTCO.N 

PA77( 2»1 >=A45*.001 

CONTINUE 

CALL  AFNC(TIME»P100*511 

CALL  AFNC(TTME»P101»52 J 

CALL  AFNC(TIME.P45»53) 

A1=S30 

A3=523 

CALL  AG 

A7=S29 


CALL 
CALL 
CALL 
CALL 


FOREST 
GASOIL 
COALMN 
MINING 


A21=A8*P86 

A2?=A12*P36+A13*P87 

A23=A16*Pe3 

A24=A19*P89 

CALL  SUPDEM 

A26=S39 
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►  3  =  S  3  8 
CALL  DEMAND 
CALL  OUTPUT 
CALL  L 
CALL  CAP  IN V 


INDUSTRIAL  3WITG 


6000 


!rt  I 
IF(  I 
PA  79 
PA79 
CCS'T 
IF  (  I 
PA79 
PA79 
PA79 

6 ICO  CONT 
IF(  I 
PA79 
PA79 
PA79 
PA  79 
PA79 
PA79 
PA79 
CONT 
IF(  I 
PA79 
PA79 
CONT 
IF(  I 
PA79 

■400  comt 


6  2  0  0 


63  00 


5:j>C.-'':»0)  GO 
COUNT. NE. 1976) 
(2*25)=17.27 
{  2»9)=61.24 

:  iu 

COUNT. NE. 1977) 

(  2  » 2  5  )  =  1 7  •  1 

(2*30)  =  18  .4-9 

(2»29)=13.95 

INUE 

COUNT. NE.1973 ) 

(2*25) =17*04 

( 2  »9 >=18.63 

I  2  >  2  )  =  3  .  8 1 

<2>3)=3*31 

(  2  *  4  )  =  8  .  8 1 

(  2  >  5  )  =  8  .  3 1 

(2*30) =17.79 

INUE 

COUNT. NE. 1979) 

(2*25)=16.18 

(  2  » 3  0  )  =  1 7  . 1 

INUE 

COUNT. NE. 1930 ) 

(  2  » 3  0  )  =  1 6  .  4 

I  U 


"C  6 400 

CO  TO  6000 


iO  TO  6100 


GO  TO  6200 


GO  TO  6300 


:-0  TO  6  40  0 


INDUSTRIAL  SWITCHES 


C 

c 
c 


CALL  ENRDM 

SUBSTITUTION  IN  INDUSTRY 

IF(  ISPRM.EQ.O)  GO  TO  777 
IF( I COUNT. LT. 1980)  GO  TO 
IF(  ICOU.N:T.GT.1990)  GO  TO 
X29=PA79 (2  »29) 
X25  =  PA79( 2  .25) 
X2=PA79(2»2 ) 
X7=PA79( 2*7) 
GO  TO  777 
776  A=PA79(2»29) 


777 
776 


CALL  SU3ST(X29»A»1.».65E 

PA79(2»29)=A 

A  =  PA79(2,25  ) 

CALL  SU3ST(X25»Ai 1 . » .65E-06 ♦ .5  * 

PA79( 2  »25 ) =A 

/ s°A79(2»2 ) 

CALL  SUBST(X2»A» .3».65E-06».5»« 
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5».5».3»1.»TIME) 


*TIM£) 


» 1 . »  T I  ME ) 


p 
p 
p 
p 

p 

il 


- 


P  A  7  9  (  2  » 2  )  =  A 
PA79«  2  >3  >=A 

PA79  :;♦•>-•  a 

PA79(2»5 1 =A 
A=?A79(2 »7  ) 

C  ALL  SUB 5T{X7»Af.3»1.125! 
PA79( 2»7)=A 
7  77  CONTINUE 

CALL     DEf'.OG 
A33= (A30-V2 l/Dl 
V2=A3C 

A3^=(A33-V1 )/Di 
V1=A33 
CALL  LAND 

C 

C  TEMPORARY  OUTPUT  Or  RESULTS 

C 
3900  CONTINUE 

IF(IPRNT2.NE.O)  GO  TO  7300 

PRINT  4000»TIME 
4000  FORVAT( • 1RESULTS  THROUGH  •  ».-6«0///  ) 

IF(!PRNT.EQ»OJ  GO  TO  7100 

KTS=15 

KTE=16 

PRINT  4010»KTS»KTE»S15 »S16 

KTS=17 
KTE  =  l-3 
PRINT  4010»KTS»KTE»S17»S18 

KTS=19 
KTE=22 
PRINT  4010 ,KTS»KTE*S19»S20»S21 »S22 

K  T  S  =  1 

PRINT  40  2  0  ♦  K  T  S  ,  ( S 1 ( J )  » Jf 1 »50  J 

<TS=2 

PRINT  4020»KT5»(S2(J) »J=1»50) 

KTS=3 

PRINT  402C»KTS»{S3U)  *J  =  1»2C) 

KT5  =  4 

PRINT  4020»KTS»(S4{ J) »J=1»20) 

KTS^23 

PRINT  4C20»KTS»S23 

KT5=5 

PRINT  40  20»KTSt (S5( J) » J=l >20 ) 

KTS=24 

PRINT  4020»KTS»S24 

KTS=6 

PRINT  4020»KTS»(S6(J) »J=1»20) 

'<TS  =  25 

PRINT  ^020»KTS»S25 

KTS=33 

PRINT  ^020  »K.TS»S33 

KTS=26 

KTE=  2  3 

PRINT  4010»KTS»KTE»S26»S27fS28 

KTS=39 

PRINT  4020»KTS»S39 

KTS-3- 


1.  .TIN  ) 
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c  p  !  M  T    4020>:<TS»(S3MJ)»J  =  1»3  i  1 

If  T  S  = 

P  R  I  o     4  0  ?  0  ,  :<  T  S  >  (  S  3  5  (J  )  *  J  =  i  »  3  5  ) 

KTS=29 

KTl^33 

PR  PIT    ^01C*KTS*KTE»S29*S30»S31 »S32»S: 

K  T  S  =  3  5 

;   IT    4020tKTS» (S35( J} »J=1>3) 

:<TS=4 

K  T e  =  -V  3 

PRINT  4010»KTS»KTE*S40>S41 »S42»S43 
4010  FORMAT ('  5'. 12,'  -  S • *  12 » T4G »6G12. 5  J 
4020  FORMAT ('  S1 »I2» (T40t6G12*5) ) 
7  LOO  CONTINUE 

PRINT  5000 
5000  FOR.VAT(  '  1  '  ) 

KTS=1 

KTE  =  6 


* 


PRINT  5  010 

<TS=7 

KTE  =  9 

PRINT  5010 

KTS=10 

KTE=14 

PRINT  5010 

KTS---15 

KTE^--17 

PRIM  5010 

K7S=13 

KTE=20 

PRINT  50  10 

KTS=21 

KTE=25 

PRINT  5010 

KTS=1 

PRIM 7  50  2  0 

K75=26 

PRINT  5030 

KTS  =  3 

PRINT  5020 

KTS  =  2 

PRIM"  5020 

KT5  =  4 

PR  I  4T  5020 

KTS=27 

KTE=29 

PRINT  5010 

KTS=30 

P  R  I  N 7  5  0  3  0 

KTS  =  5 

PRINT  5020 

KTS--6 

PRINT  5020 

<TS=31 

<TE=32 

PRINT  5010 

KTS=33 


KTS»KTE»A1 »A2  »A3»A4»A5»A6 
KTS»:<TE»A7»A8»A9 
KTS»KTE»A10»A11*A12»A13»A14 
KTS»KTE»A15»A16  iA17 
KTS »KTE  »A 1 S  »A19  *A20 
KTS»KTE»A21 *A22  »A23 »A24 »A25 

kts, (Ai it j ) »j=i »io) 

KTS»A26 

KTS,  (AI3( J)  »J=1  ,25 ) 
KTS,  { AI2( J)  »J  =  1  ,25) 

KTS,  ( AI4( J)  »J  =  i  ,35  ) 

KTS,KTE,A27  iA23 »A29 
KTS,A30 

KTS,(AI5(J),J=1,5) 
KTS,:aI6(J)»J=1,3p) 

KTS,KTE»A31 ,A32 


B-12 


i 


5010 

5  0  2  0 
5030 
7  2  0  0 


7  3  GO 


740! 


<T  =  34 

[NT  501C»KTS»KTE»A33»A34 

KTS=3i 

p'->  I  \T  5310  »KTS,KTE  *A36  »A37  »A38  »A39 

KTS=40 

PRINT  5030»KT5»A40 

KTS=41 

K  T  E  =  4  3 

:  \!T  5010  »KTS»KT£»A41  *A42  >A^t3 
FORMAT ('  A«»I2»'  -  A1 » 12 »T40 ,6G12.5 ) 
FORMAT ( '  AI • »I2» (T40»6G12*5 ) ) 
FORMAT t '  A' * I2»T43,G12»3) 
CONTINUE 

IF(TIME.GE.TEND)  STOP 
TIME=TI-ME+DT 
ICO'JNT=TIME+.00001 
GO  TO  300 
CONTINUE 

PR  I. NIT    7400»TI.ME»AI5 
GO    TO    7200 

FORMAT ( lX»F6.0/«     AI5'*T20,6G14*5) 
END 
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C    S        IOUTINE    AG 

SU9ROUT  Ii\E     A3 
COMMON    /AG!/    AL»P»WR 

COMMON    /AGP/    PL,PO»OLP»FTR»RA,PLF,PTI »TF >SALV»ALEV I »AKDEP»Al  P >A 

C  J  '  V  0  N    /AGS/  A  <  0  >  A  !<  F 

COMM<    ,     /.       iRLI/     ft  IF  .:  >Y»TI  »XYZ123(  103) 

common   /genrlp/   xyz12m2)»dt 
c  agriculture  output 
tk=a:<o+akf 

Y=A*T;<**ALP*ALi 

ALR=Y/OLP 

ALE  =  ALR*V/R 

ET=FTR*AL#PL 

DOE=Y*PO 

3=RA*AKF 

C=PLF*AL*PL*RA 

E  I  =  n,  +  C 

D=Y*P 

ANFR=D-< DOE+EI- 
C  CAPITAL  FORMATION 

sx  =  o 


kLE+ET  ) 


I 


!F(ANFRiGT»0)    AX=ANFR*PTI 

IF(ANFR.GT.-TF*Y)     GO    TO    100 

SX  = f-ANFR~{ TF*Y) ) /SALV 

DX=»AX*ALEV! 

EX=A<F*AKDEP 

FX=AKO*AKDEP 

GX=FX+5X 

AKO=AKO+ ( AX-GX ) *DT 

AKF=AKF+(DX~EX !*DT 

TI=AX+DX 

A\FR!=ANFR-AX 

RETURN 

END 


« 
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C  SUBROUTINE  -ORES' 

SUBROUTINE  FOREST 

COMMON  /FORSTI/  TT»CR 

.  v.MCN  /FORSTP/  AC>FCR»DR>CT 

CO '•' XO.N  / PC R 5 7 S / T AC  * C ,\-> 

COMvo?'  /GENRLI/  XYZ123  {  3  )  »CFR  »XYZ124-  { 107  J 

COMVON  /GENRLP/  XYZ125(2)»DT 
C  FOREST  OUTPUT 

TC=AC*TT 

TCR=TAC*FCR 

ACN=TCR-CAP 

IF( ACN.LT.O  )  ACNaO 

DEP=CAP*OR 

CFR=(ACN/CT  )+DEP 

CRM=CAP/FCR 

DIF=TAC-CRM 

IF(D!F".GT.O)  CR=CRM 

IF(DIF.LE.O)  CR=TAC 

TAC=TAC+(TC-CR)*DT 

CAP=CAP+ (CFR-DEP ) *DT 

RETURN 

END 
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C  SUBROUTINE  GAS  IL 


C 

c 


I 


DO 


25 


5  0 


SUBROUTINE     GASDIL 

COVVON    /GASOI/    C^»GD,TGD 

COMMON    /GASO=V    GOR»OGR*CCP»DC*GOOE*ODGE»DPT»CCOISCR 

iOS/AFO»AFG»OOI  »GOI ♦ PRO (50) ,PRG(50)  .CAPO (23)  >CAPG(20) 
COMMON    /  LI/    XYZ123(4) »TOP,TOGI  »XYZ12M1C5J 

ERATIONS  INVOLVING  SUMMING  OF    PRODUCING  FIELDS 
-C'.'ANTITY  LEVELS 
ICAO=l. 
ICA3=l. 
TDD=AFO 
TOG=APG 
CO  3  0  1=1 » 50 
J  =  5  1  -  I 
A=PRO{ J)  /2» 
TD0=TD0+A 

IF{ ICAO.NE.l)  GO  TO  2  5 
CALL  AFNC(TDO»PCMOt54-) 
0°AQJ=GP/PCMO 
CALL  AFNC(0PADJ»X>4) 
ICO=X+l. 

IF( ICO.GT.J)  ICAO=ICC 
TDO=TDO-t-A 
A=PRG( J)  /2  . 
TDG=TDG+A 
IF<  ICAG-ME.l  )  GO  TO  50 

c  a  l  l  a  f  n  c  (  7  d  g  » p  ov  g  »  5  5  ) 

goadj  =  gp/pc-':g 

CALL  AFNC(GPADJ»Xt8J 
ICG=X+1. 

IF(  ICG.GE,J )  ICAG=ICG 


TDG=TDG+A 

AF 


7DO»PCVO»5^J 

:mo 

\'ccopadj>x,a) 

MG»55 ) 


CALL  AF 

ICO=X-4-l  . 

CALL  AFMC(TDG» 

GPADJ=GP/PCMG 

CALL  AFNC(GPADJ»X»8) 

ICG=X+1. 

TDO=AFQ 

TDG=A=G 

DOP=0. 

CGP=0. 

OIR=0. 

G  I  R  =  0  . 

OILCA°=0. 

GASCA?=0. 

C  IL3EK  =  0. 

GASPiX=0. 

OILINC=0. 

GASINC=0» 

DO  100  1=1.50 

J  =  5 1  - 1 

X  =  J 

D I SCNT  =  E XP ( -X*0 1  SCR ) 
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li 


4 


C 


« 


i: 


RO 


c 

r 


A =DR  0 ( J ) 
TDO=T0 3+A/2* 

IF(  J«"T.  !C  ■  JO    TO    60 

call   -\"-:c  (  .;»v  >1  ) 

DOP=00:>+Y*A 

CALL    AFNC  {  TOO  >  7.  ,  r>'-  ) 

CALL    AFNC(X»Y»3) 

CIR=OIR+Y*A*Z 

TOO  -A/2- 

I : ! J.GT. ICO)     GO    TO    70 

CALL    AFNC(X»W»1] 

CALL     AFNC(X*Y»2) 
L    AFNC(X>Z»3) 

0  I  L  I  NC=0  I  L  I  NC+W*OP*D  I  SCMT 

OILPEX=OILPEX+Y*DiSCNT*PCM0 

0  I L CAP=0 1 LCAP+Z *D I SCNT*PCMO 

A=PRGLJ) 

TDG=TDG+A/2« 

IF! J.3T. ICAG)  GO  TO  80 

CALL  AFNC(X,Y»5) 

DGD=DGP+Y*A 

CALL  AFNC  (TDG»Z»55  ) 

CALL  *FNCU»Y»7) 

C  I3»GIR+Y*A*Z 

T0G=TDG+A/2. 

IF( J*0T* ICG)  GO  TO  90 

CALL  AFNC  (X  »W»5) 

CALL  AFi\C(X»Y»6) 

CALL  AF.NC(X,Z*7) 

GASINC=GASINC+W*6P*DISCNT 

G A Sp EX = G ASP EX- Y*D I SC^T^PC-G 

GASCAP=GASCAP+Z*DISCNT*PCMG 

TDG=TDG+A/2» 

CONTINUE 
TOTAL  CAPITAL 

CAPOT=0. 

CAPGT=C* 

DO  110  I=lt?0 

CAPOT=CAPOT+CAPO( I ) 
11C  CAPGT=CAPGT+CAPG{ I ) 
CAPITAL  EXPENSE 

CEO=CAPOT*CCP 

CEG=CAPGT*CCP 
PRODUCTION 

GIP=GOR*DOp 

OIP=0GR*DGD 

TOP=DOP+OIP 

TGP=DGP+GIP 
INCOME  FOR  PRODUCTION 

OIN=DOP*OP+GIP*GP 

GIN=0GP*GP+OIP*0P 
PRODUCTION  PROFIT 

0PP=0ILI NC-OILPEX-OILCAP 

G?P  =  GASl>NC-GASr>EX-GASCAP 

EXPLORATION 


9  0 
100 
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c 

c 


.  .   .   :    .  "•    ACTIVITY 
CDR  =  0DI/OP  r*DCJ 
GO  i  =  GDI / CDPT*DC) 
DISCOVERY    RATES 

CALL    AFNCJ  TD0»DR0»9) 
CALL    AFNC(TDG»DRG*10) 

DD=ODR/DRO 
DGD=GOR/DRG 

*-  f>  r>  —  ^  r>  <">  t  *~  P\  O  C 

- 

TOD=DOD+ODG 

TGD^DGD+GDD 
1  INCOME •  EXPENSE  RETURN 

OE I  =DOD*OPP+GDO  *G  .D  P 

GEI=DGD*GPP+ODG*OPP 

CEPO=O0I  *CC 

CEPG=GDI*CC 

PREO=OEI-CEPO 

PREG=GEI-CEPG 

PROL=PREO/ODI 

PRGA=PREG/GDI 

CALL  AFNC(PROL»Y»ll J 

RIO=Y*C  : 

CALL  AFNC(PRGA»Y»11 ) 

R  IG  =  Y*GDI 
CAPITAL  ' DEPRECIATION • 

DPPO=ODI/0PT 

DPRG=GDI/OPT 
UPDATE  LEVELS 
PRODUCING  FIELDS 

I F { ICAO.GE*50) 

IF! ICAG.GE.50! 

CO  120  I ^2  »  50 

J=50-I+2 

IF { J.LE. I CAC )  PRO ( J ) =PRO ! J-l ) 

IF< J.LE. ICAG)  ?RGi J ) =PRG( J-l ! 

IF( J-l.EO. ICAO)  PRO( J) =  ?ROUl+PRO(  J-l ) 


AFO=AFO+PRO(50) 
AFG=AFG+PRG ! 50 ) 


ir 


(  J-l 


ICAG)  PRG( J>=PRG( J)+PRG( J-l ) 


13. 


120  CONTINUE 

PRO< l)=TOD 

PRG( 1) =TGD 
PRODUCTION  CAPATAL 

DO  130  1=2*20 

J=20~I+2 

CAPO< J)=CAPO(J-l  ) 

CAPG< J)=CA°G( J-l ) 

CAPO< l)=OIR 

CAPG< 1 )=GIR 
EXPLOR A  T ION  ! N VESTMEN T 

ODI=ODI+RIO-DQRO 

GDI=GDI+RIG-DPRG 

TOGI=OIR+GIR+RIO+RIG 

RETURN 

END 
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i  F  I  R  P 


>  s 


i.  .  ~ 


150 
100 


C  SUBROUTINE  COAL-' 

SUBRO'JT!     rOALM.N 

COMMON  /COALMI/  PRCS*PRST 

rr'-"QM    /COA'-iYP/CO.Kl  >K2»K3  »FDR? 

fOV^ONi    /COALVS/TM»PRC#AIN  (20  ) 

CCv-lON    /GENRL!  /    XYZ123  16)  »AO,AI  >XYZ12M103) 

CALL     V-  \C(  T'-'tCI  » '<  3  ) 

CALL    •'   '  !C(  TM»ACM*K1  ) 

As? ?C/ACv 

CALL    M  :NC(A»FO*K2 ) 

AACv=0. 

X  =  0. 

L=0 
100    CALL    AFNC(X»Y»K4) 

AACM=AACM+Y 

X  =  X  +  •  0 1 

IF(X.GT.FO)    GO    TO 

L=L  +  1 

IF(L.LE«99)  GO  TO 
150  AACM=AACM/L 

TI=0. 

DO  200  1=1 »20 

2  00  TI=TI+AIN<  I  ) 

CL0=TI/CI 

AO=CLO*FO 

GI=AO*?RC 

E.M=AO*AACM#ACM 

EC=TI*CC 

P=GI-E=V-EC 

Ir (TI .LE.0.0)  PRsO.O 

IF(TI.GT.O.O)  PR=P/TI 

CALL  AFMC(PR>RI »11) 

A!=RI^TI 

IF(AO/PRST.LT.0.5J  PSE=1.0 

IF(AO/PRST.GE.0.5)  PSE=1 .5-AO/PRST 

IF(D5E.LT.0o0)  PSE=0 

AI=AI*PSE 
C      MARKET  ENTRY  EFFECT 

USM=PRST-AO 

IF(USM.LT.0. J  USM=0* 

p  R  F = P  R 

IF(PRF.LT.0» )PRF=0 

IF(PRF.GT.10. J  PRF=10 

A  I=AI+USN!*CI*PRF*ENTRY 
C   UPDATE  LEVELS 

TMaTM+AO 

DO  330  1=2.20 

J=2^-I+2 

3  00  AIN( J)=AIN( J-l) 

AIN(1)=AI 
C       PRICE  ADJUSTMENT 

FCPD={AO-PRST)/PRST 

DPRC=0.0 

I F ( FORD.GT .0. )  DPRC=-FOPD*FDRP*PRC 

I F ( FOPD. LE .0. 0 )  DPRC= ( PRCS-P ?C ) *FI RP 

PRCsPRC+DPRC 


*Y 
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PI 
II 

VI 

il 


4 


sun 


IN 


SJ:< ^ojr i ,\f:  MitNii 

COMMON  /MINEi/  PIC  MO 

COMMON  /Ml  MEP/CCjKl  ,<2  »K3  »K4 

COMMOM  /MINES/TM»AIN (20) 

COMMON  /GENRLI/  XYZ123(3) »AI »XYZ12M 102) 


100 


150 


*i  '^  w 


:all  afnc< 


<3  ) 


CALL  AFNC(T 

•1»A 

1  ) 

a=prc/acm 

CALL  AFNC(A 

»FO 

»:<2 ) 

AACM=0« 

X=0* 

L  =  0 

CALL  AFNC(X 

»Y» 

<4) 

AACM=AACM+Y 

X=X+.01 

IF ! X*GT»FO ! 

GO 

TO 

150 

L=L  +  1 

IF(L*LE.99 ) 

GO 

TO 

100 

A ACM = A ACM /L 

TI=0. 

DO  200  1=1 » 

20 

TI=TI+AIN( I 

) 

CALL  AFNC(7 

I  »Z 

»56) 

CI=CI/2 

CLO=TI/CI 

AO=CLO#FO 

GI=AO*PRC 

EV=AO*AACM*AC 

FC=TI*CC 

p=Gl-EM-EC 

I  r ( 71  .LE.O. 

0) 

PR  =  C 

.0 

IF(TI.GT.O.O)  PR=P/TI 
CALL  AFNC(PR»R1 .11) 
AI=RI*7I 


C   UP 


>H   I 


E  LEVELS 


TM=TM+AO 
DO  300  1=2  i  20 
J=20~I+2 
3  00  AIN( J)=AIN( J-l 
AIN( 1)=AI 
RETURN 
END 
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C  SUBROUTINE  SUPDEV 

SUBROUTINE  SUPDEM 

CC'-'J'.  /SOEMI/  FROUT  »ONOUT»CLOU1  i  WOuTtAIMP 

CO-"-*M3N  /SDE.VP/P.MV.R  .FRMAjSCP 

COM-10.N  /SDEMS/PP.ROD 


'J-  <'■   '-1/  XYZ123.AGOUT,XYZ12M2)  »OILP *XYZ125 (9) .OSJ i  fM  , 

•  XYZ126<84)  "'"J,-J' 


100 


200 


3  0: 


OSU)  ^AoOUT 

CS(2) =FROUT 

OS(3)=Ci\OUT 

OSU  J=CLCUT 

C5! 5 )=OMOUT 

05(6)  =OiVOUT*PMMR 

IF(0ILP*FR.MA+AIM?.LT.PRR0D)  GO  TO 

IF(OILP*FRMA.GT.PRROD)  GO  TC  200 

OS(7) =PRROD 

GO  TO  300 

OSC  7 ) =01 LP*FRMA+A I  MP 

GO  TO  300 

0S(7)=0ILP*FR.MA 

D:^03^0ILD^F^;-, 

CONTINUE 
OS(7)=OS(7)*BOP 

RETURM 
END 


100 


' 
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c  sua  .   :  v  -  ■  .  :o 

SUBRC  n       !A.NO 

CO.yMON/DE.V.\D?/FS»SEXP,FEXP,EPP,RCC.REXC.IRAlLfPED(25),EX(25]  . 
?   SEO(25)»FED(25) 

COMMON/GENRLI/XYZ12316)  *  COAL  ,XYZ  1 24  !  2  )  ,AVHHI  »PT »XYZ128 < ? 1  , 

1  HT»XYZ125<  10), 

2  DM ( 25 ) »CAP I  2 5 ) , X YZ 1 26 ( 3 5 ) »  FG  »SG 
CO  wo?]    /DEMND!/       AVOUT 

COMMON    /GENRLP/    KD»XYZ127<2) 
EXPD=A>/HHI*HT*FS 
SG=AVOUT*SEXP 
FG=AVOUT*FEXP+PT*EPP 
CO    133     I=1,KD 
100    DM(I)«EXPD*PED< I)*(1.+EX(  I  )  )+SG*SED( I  >+FG*FEDC I >+CAP<  »  ! 
DM( IRAIL)=DM( IRAIL)+COAL*RCC+REXC 
PET' 
END 
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SU3R0UT  I  ■     ITPUT 

SUBROUTINE  OUTPUT 

:  /OUTPTP/  KS»ADS(25*10)»AD(25»25J 
-S^f'Cf-'-  W    XYZ^3<i*>'OSClO),DMC25),XYZl24<25>,ODC3P), 

co  •  ■  v.  i"  n  /  g  e  n  r lp /kd » k t » d t 
common  /outpti/  to 
co^n'on/outpts/avout 
dimension  a  (2-3) 

DO    110     1=1 tKD 

=  0. 
DO    103    J=1»KS 
6=B+AD5< I »J)#OS(J) 
M  !  )=3+DM<  I  ) 
DO    210     1=1 »KD 
3=0, 

DO    2  00    J=1»KD 
B=B+AD( I , J ) -At U) 
00  ( I )=B 
DO    300     1=1  ,<S 

J =:<  D+ 1 

od( j) =osi r ) 

TP  =  0, 

CO    400    1=1 »KT 

TO  =  TO+OD( I  J 

AvnUT=AVOUT+(TO-AVOUT)*DT/3 


1  C  0 

no 


2  00 

7  1  ': 


3  00 


4  00 


FND 
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C    SUBROUTINE    LABOR 

SUBROUTINE     LABOR 

COMM  )N    /LABOR  I  /    TE 

COMMON    /LABOR?/    W?R »USO »USP »USA »FLH »OER( 35)  »ILB(35) »FESR»SESR 

COM  ION/LA30RS/DV/»PW»HAL»AVHHI 

COMMON/GENRLI/ANFRI  »XY7.123  (9  J  »PT»UEM»WRM*HT  *XYZ124(  60  )  »OD(35)  , 

;  S  t  S  G  S 
C  V  p 


CO"  • 

/GENRLP/ 

LOYMENT 

AE=0 

PF  =  D 

DE=0 

DO  100 

1=1 »KT 

AE=AE+OD( I )*OER( I ) 
PE=PE+OD(  I  )*OER{  I  ) 
DE=DE-K)D( I )*OER( I ) 


m 


IF( ILB( I ).EQ»1) 
IF(  ILBU  ) .EQ.2 ) 

IF  (  ILB(  I  )  .E0.3) 
100  CONTINUE 
C  UNEMPLOYMENT  AMD  WAGE  CHANGES 

WP=PT*WPR 

PE=PE+FGS*FESR 

DE=DE+SGS*SESR 

TE=AE+PE+DE 

UEM=(WF-TE )/WF 

E0=0*05-UEM 

IF(ED,LT.0„0)  ED  =  0 

DOW=E0*USD*DW 

DPW=EO*USP*PW 

DHAL=ED*USA*HAL 
C  WAGES 

AWNHaANFRI /AE/t 1,-FLH) 

AW=AWNH* ( 1 .-FLH)+HAL#FLH 

T  I  =  AV*'*AE+PW*PE+DW  *DE 

i/.:RV=TI/TE 

HHI=TI/HT 
C  UPDATE  LEVELS 

DW=0W+9DW*DT 

PW=PW+DPW*DT 

HAL=HAL+OHAL*DT 

AVHH I =AVHH I + ( HH I -AVHH I ) *DT/3 

PETURN 

END 
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c  sua 

su: 
co- 
co- 

CO' 

2  X' 
CO' 

ov 

CP: 


I  ME 

ROUT 

■ 

'•'ON 
XON 
Z  1  2  6 

-3  0  1 


:  w 

:  N   CAPIN'V 

/CAPIiNI/  005(35) 

/CAPIKP/ATD»DPR(35  1 »CAR(35) > IM( 35 ) >CPS(25) 

/CAPIN5/RC< 3  5)  »A0( 35 ) 

/GE.NRLI/  XYZ123(2)  iTI  ,CFR  ,  XYZ 124  ,  TCG !  ,XYZ125,ACI  »AM1  t 

(40  )  .CAP (25 )  »OD<35 )  »XYZ127(2 ) 

/GENRLP/  <  •   .DT 

O.N  DAO(35)  »  01  -0(35)  »DIFR(35)  iDRC<  35)  >RI (35  )  .000 1  35  ) 

I  =  1  *  <  T 


c  =  o 


IF( IN( I ) *EO.O)  GO  TO  5  00 

DIFOt I >=0D( I )-A0( I ) 

IF(A3S(00S( I) ) *LT.1000)  005 < I )=1000 

0AO( I )=DIFO( I ) /ATO 

000 ( I )=(0D( I )-ODS( I ) )/(0D5( I )*DT) 
DIFR(  I  )=DOD( I >-RC<  I  ) 
DRC< I J=DIFR( I ) /ATO 
RI (I )=RC(I  >+QPR<  I  ) 

Ir(RKI)  .GT.OJ  C=AO<  I  )*RI  {  I  )*CAR(  I  ) 
UPDATE  LEVELS 

RC( 1 )=RC< I >+DRC( I >*0T 
A0(  I  ) =A0{  I  )+0A0(  !  )*DT 


500 

CP=CP+C 

CP=CP+TI+CFR+TOGI+Af 

llt-ACI 

00    6  0  0     1=1 »K0 

600 

CAP( ! )=CP*CP5i I ) 
00    700     1=1 »KT 

7  r%o 

005 ( I  )=O0(  1  ) 
RETURN 

END 
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c 

c 
c 
c 
c 


SUBROUTIN 

SUBROUTINE  ENRDM 
COMMON  /EN    :/  ENDEM(  5) 

COVMON  /ENRDMP/  MENR ,  IEND *AR ( 5 . 10 >  » SAT ( 5  ♦  10 ) »AI (5>35 ) 
COMMON  /GENRLI/  XYZ123 ( 13 ) »HT >XYZ124 (60  )  .00 ( 35) »XYZ127(  2 ) 
COMMON  /GEMRLP/  XYZ125 »KT »XYZ126 
•  1 —  ELECTR IC I  TV — KWH 


iAS-- 


• ?T#4^ 


c 


3 OIL 3TU 

COAL 4 BTU 

5 OTHER 3TU 

DO  ICO  1=1  »N< 
100  ENOEMI I ) =0 
RESIDENTIAL 

DO  2  00  I=1»NENR 

DO  200  J=1»NEND 
200  ENDEM<  I  )  =ENDE;M<  I  ) +AR  (  I  »J)*SAT<  I  »J)*HT 
INDUSTRIAL-COMMERCIAL 

DC  300  I=1»NENR 

DC  300  J=l iKT 
300  ENDEM( I )=ENDEM{ I )+AI { I »J)*OD(J) 

RETURN 

END 
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C  SUBROUTINE  DEMOG 

SUBROUTINE  DEMOG 

COMMON  /DEMOGI/  UENtta  ■: '. 

COMMON    /DEM0GP/UMC,WMC»IAM1  »  IAM2»IAF1  ♦  IAF2»S.MR,  >.AD»RM.(4]  »AMT  (25)  . 
wR(40) »DH(30) ,Bh( 5»45) »IS1»IB2*EhR 

CO  :MON    /DEMOGS/    PC  (6,2  5)  »PFF{60)  ,PYF(50)  ,PSF(50)  ,PSM(  50)  ,PY.M(60)  <. 
2       PF.M(50  )  »ELD(  30) 

COMMON    /GENRLI/    XYZ123  { 10  )  ♦?  »UEN'»V»  -  it  17  >XYZ124(  94  ) 

DIMI    \'SI    IN     PT(90)  »PTF(90)  > '-' T  M  <  90)  >BHT(6)  »:••"■'(  5  )  »HRM(4) 
C  -:f  »  ::•  *  *  »  ->  *  *  -  *  *  *  --  »  *  *  *  *  '■<■  *  #  *  **  *  »  *  »  *  *  *  »-  «  >  ■::■  *  »  :-  ~-  *  -*  *  -»  *■  *  ->  «■ ;;- ,- »  ■  fr  «-  *#; 

WMC1=WMC 

IF(UEM-UEN«G7.0.02)    WMC1=0«0 

As(1.01  +  (UEN-UEM)*UMC«M  t  WRM-WRN  )  /WRN  >*WMC1  ) 


HT=HT*A 
C  -x-  *  *  *  if  *  *  #*#  -:'-  *  K-*  #  #  *- £  *  *  -:<-  *  *  #  «■  «■  *  *  ;;- 

RETURN 

END 


■###*#«■##*•##**#**##*■*■#«■  #*#  ■<■<.;-;;•-.;- «-*-»*»#* 
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r 
C 

c 
c 
c 
c 


c 
c 
c 

c 

c 

c 
c 
c 
c 
c 
c 
c 


SUBROUTINE  LAND 

SUBROUTINE  LAND 

COMMON  /LANOI/  THHC »EMPCO*REFOR »PRESi 
COVMON  /LANDP/TONACR >REVEG » CD» EMPACR >RVYR3 
COMMON  /LANDS/FL»AGL  iU.NL  »URL  »PRL>  A.ML  (  3  ! 
COVMON  /GEKRLI /  XYZ123 ( 6 ) »COALPR »XYZ 124 { 1 
COVMON  /GENRLP/  XYZ125(2)»DT 
i.LPR-CDAL  PRODUCTION  TONS /YR) 
TSHC-",-:/,1  HOUSEHOLDS/YR 
:  >■:  CO-EMPLOYMENT  C:J-".N  3E  1JOB5/YR  ) 

F  C R-R E  FOR  ES  TAT  I  O.N  C  ACR E / YR ) 
PRESET -L  A '; 0    PRE SERVED  (  ACRE/  YR  ) 
PRL-PRESERVED  LAND(ACRE) 
[ST  LAND (ACRE) 


r  L- 


:OR^ 


ER  YEAR 


AGL-AGR I  CULTURE  LAND 

A.ML(3)/LAND  MINED  BEING  RECLA  I.MED  (  ACRE  )  3HRD  ORDER  DELAY 
U.NL-UNUSED  LARD  (  ACRE  ) 
URL-UR3AN  LARD (ACRE) 
TONACR-TONS  CAOL/ACRE 
REVEG-FRACTION  OF  LARD  MIRED  WHICH  CAR  3E  RECLAIMED 
CD-HOUSEHOLDS /ACRE 
EMPACR-EMPL0YEE5/ACRE 
RVYRS-YEARS  TO  REVEGETATE 
ACRLOS-ACRES  MIRED  PER    YEAR 
UR3LOS-ACRES  TO  URBAN  PER  YEAR 
ACR5T-ACRES  RECLAMATION  STARTED  0,v 
ACRGN-ACRES  RECLAIMED  PER    YEAR 
UMLOS-ACRES  TO  UNUSED  PER  YEAR 

ACRLOS=COALPR/TONACR 

UR3LOS=THHC/CD+EMPCO/EMPACR 

UNLOS=ACRLOS*( l.-REVEG) 

ACPST  =  ACRL05^-REVEG 

ACRGN=AML ( 3 ) / ( RVYRS*3. ) 
UPDATE  LEVFLS 

FL  =  FL-r(REFOR-DRESER)*DT 

AGL=AGL+(ACRGN-ACRLOS-URBLOS-REFOR ) 

UNL= UN L+UNLOS*DT 

URL=URL+UR3L05*DT 

PRL~DRL+PRESER*DT 

AMU 3)=AML ( 3)+(AML(2 )/ (RVYRS*3. J-ACRGN ) *DT 
AVL(2)=AML(2)+{AML< 1 )/(RVYRS*3. I-AMLC2 )/{RVYRS*3. ) )*DT 
AVL(  1  )=A.ML<  1  )+(ACRST-AML(l  >/(RVYRS*3.  )  )*DT 

RETURN 
END 


-DT 
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iU3R0iJTIN! 

SUBR0U7 : 


HHTCON 


COMMON  /HHTCMI/  PRCIN.DHH.AHE 

COV.VOM  /GENRLP/  XYZI23C2  )  ,D7 

COMMON    /H  ITO\P/    i-ASE,U3-3,USC,Ct3;),CCO,CB\»rr;,M:J< 

CO"i>.".ON    /HHTCNS/    H5  A  ♦ '! ^r!  .  HSC  »  PR  I  CL" 

TIAL     HEATING    ENERGY    USE 
—  D/(PRICE*3ASE*< 1.-US3) ) 

P3C3*(CC0-CBO)/(PI    ;  3ASE*(US3-USC)1 

P3CA=CCO/ (PRICE*3ASE*< 1 .-USC) ) 
C 

P3B=C3N/(PRICE#3ASE*(1«-US3J  J 

c8C  =  CC.\/(PRICZ-BASH*(l.-U3C)  ) 
C 

HRC=HSC/HRR 

HRB=HS3/HRR 

HRA=H5A/HRR 

RNH=DHH+HRC+HR3+HRA 


CALL     \FMC{P33A»A»57) 

RHA3C=HSA*A 

CALL  AFNC{P3CA.A,57) 

RHAC-HSA*A 

RHA3=RHABC-RHAC 


CALL    AFi\iC(PI 
RH3C=HSB*A 


v.  3  »  M  > 


57) 


CALL  AFNC{PBC»FC,53) 

CALL  AFNC(PB3»A»58) 

FA=1*-A 

F3---A-FC 

RMHA  -r<:;n*FA 

RNH3=RNH*F3 

RNHC=R,\H*FC 

CALL  AFNCfPRICE.A.59) 

EHHA=BASE*A 

CALL  AFNC(PRICE*A»60) 

EHHB=BA5E*A*USB 

CALL  AFNC( PRICE. A»61 ) 

EHHC=BASE*A*USC 

TH£=EKHA*HSA+EKHB*HS3+EHHC*H5C 

THH=HSA+HSB+HSC 

AHE=THE/THH 

HSA=HSA+(RNHA-HRA~RHABC) *OT 
HS3=HS3+ ( RNH3+RHA3-HR3~RHBC J  *DT 

HSC=HSC+(RNHC+RHBC+RHAC-HRC)*DT 
PR  ICE  =  PR  I  CE+(  PRCI  \-PR  I CE  )  *DT/2  .C 
RETURN 
END 
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C  SUBROUTIN      ;ST 

SUBROUTINE  SU3ST  (  XI  ,X  »F,CC  »EFFN  ,!  •  i  0»E 
C  SUBSTITUTION  i:.  THE  INDUSTRIAL  SECTOR 

CALL  A-.NC1  Tiv.E,PRE»65) 
•  CALL  AF.NC(TIME,PRS»66) 

S  -W  =PRE/EFFO- ( PRS  /EFFiN  +CC  )  *EFR 

SAV=SAV/(PRE/EFFO) 

CALL  '  ""  C  (SAV»RC»S7) 

A=X-X1*<  1  «,-,-) 

A=A*( l.-RC ) 

X=Xi*(i.-F)+A 

RETURN 

END 


FRiTIHE) 
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SUBROUTINE  AFNC 

SUBROUTINE  AFNC  (  XZ*Y»IZ) 

COMMON  /AF.NCP/  XFNC(75»30»2) »IFNC(75»5) 


w  a  -K-  "-  ii-  ''.    .  .  '.;-  >  v 


c   8- #  •*•**«■  •»•»■*•« 

C   REASSIGN  THE  INPUT  ARGUEMENTS 

x=xz 


■"- .  •  -> 


■■  ■"-  8  .  K- 


r 


rs    is   a  rfHicn   i  \terpol  \ 


7ES    :■     ' . 


.ED  fu  .z:  id: 


c 
c 
c 
c 


c 
c 
c 
c 
c 

r 

C 

c 


VALUES  BY  FITTING  A  P  OLYNCM  I AL  TO  A  NUMBER  OF  POINTS 

X  IS  THE  VALUE  OP  THE  INDEPENDENT  VARIABLE 

Y  IS  THE  VALUE  OF  THE  DEPENDENT  VARI ABLE  CALCULATED 

I  I  S  THE  FUNCT  ION  NU.M  }  ; 

TWO  ARRAYS  MUST  BE  SUPPLIED  TO  IT  THROUGH  COMMON 

THE  USER  MUST  SUPPLY  A  COMMON  STATEMENT  SUITABLE  FOR  THE  FUNCTIONS 

BEING  USED  " 

XFNC(I»J.K)  IS  THE  ARRAY  WHICH  CONTAINS  THE  TABLED  FUNCTION 

I  IS  THE  NUMBER  OF    THE    FUNCTION  BEING  USED 

J  IS  THE  INDEX  FOR  FUNCTION  POINTS 

K  IS  1  FOR  THE  INDEPENDENT  VARIABLE  AND  2  FOR  THE  DEPENDENT  VARIABLE 

IFNC(I«J)  DESCRIBES  THE  NATURE  OF  THE  FUNCTION 

I  IS  THE  NUMBER  OF  THE  FUNCTION 

THE  J  VALUES  ARE  USED  TO  CONVEY  INFORMATION  ABOUT  THE  FUNCTION 

1  STANDS  FOR  YES,  0  STANDS  FOR  NO 

J=l — EXTEND  FUNCTION  BELOW  MINUMUM  VALUE  OF  INDEPENDENT  VARIABLE 

J=2 — EXTEND   FU  CTION  BEYOND  MAXIMUM  VALUE  OF  INDEPENDENT  VARIABLE 

J=3 — EQUALLY  SPACED  VALUES  OF  INDEPENDENT  VA  IABLE 

J  =  4 — IS  USED  FO       THE  NUMBER  OF    DATA  POINTS  IN  THE  FUNCTION 

J  =  5 — IS  USED   0   THE  NUMBER  0   DA    POINT   TO  BE  USED  I  .  T Ht  IN  ERP 

IERRO   IS  A  SIGN  L  USED   0  INDICATE  A  MALFUNCTIO   IN  THE  IN  ERPOLATI 

3  INDICATES  NO  ERROR  o   1  IN  I  CAT    AN  ER    OR 


REASSIGN  OFTEN  U  ED  AR   Y  VALUES 
J4=IFNC(  I  >4) 
J5=IFNC ( I »5) 
X1=XFNC( I »1»1) 

xn-xfnc: i »J4»1 ) 

INITIA  IZE  ER  OR  PARAMETER 

I E  R  R  0  R = 0 
CHEC"'  TO  SEE  IF  X  IS  IN  THE  PROPER    RANGE 

IFU.GE.X1  )  GO  TO  103 

IF( IFNC( I »1 ) *EQ.l)  GO  TO  102 

V. R  I  T  E  «  6  ♦  1  0  1  )  I  »  X 
101  FORMATC24H  BELOW  RANGE  OF  FUNCT  I  ON  t  I  4*G15.4  J 

IERR0R=1 
10?  Y=XFNC( 1*1*2) 

RETURN 

103  IF(X.LE.XN)  GO  TO  106 

IF( IFNCl I .2) .EQ.l)  GO  TO  105 
W  R  I  T  E  (  6  *  1 0  M  I  » X 

104  FORMATC24H  ABOVE  RANGE  OF  FUNCTION » 1 4 »G15.4 J 
IERROR=l 

1C5  Y=XFNC( I »J4»2> 


RETURN 
FIND  STARTING  POINT  FOR  SEARCH  FOR  INDEPENDENT  VARIABLi 
105  IX=IFIX( ( (X-Xl ) / (XN-X1 ) ) *FLOAT ( U4-1 )  )  +  l 


:xe 
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I F ( I  X . EQ • J 4 )  I  X  =  J  4  - 1 
C   ARE  VALUES  OF  IN    .  1DENT  VARIABLE  EQUALLY  SPACED 
IF{ IFNCt I »3) .EO.O)  GO  TC  11C 

:  ■        '.  or  odd  number  of  data  points  to  to  3E  used 

1065  A=FLOAT ( J5 >*.5 

[  A  »  I  [  X  ( A  +  •  1 ) 

:  (A.GT.f="LOAT(  I  A)  )  GO  TC  107 

IX=IX-IA+1 

GO  TO  150 
C   WHICH  PART  OF  ::      \L  IG  POINT  IN 
107  J=0 

XI  1  =  XF.\C<  I  »  IX+1»1  J 

IF( (XI 1-X)/{XI1-XFNC<  I  »IX,1  )  )  .LT.0.5  )  J=l 

IX  = IX-IA+J 

GO  TO  150 
C   SEARCH  FOR  INDEXES 

110  J=0 
K=0 

111  IF(X-XFNC(  I  >IX»1) )  120»120»125 

120  IF(J.EQ.O)  GO  TO  121 
IX=IX-1 

GO  TO  10  55 

121  K  =  l 
IX=IX-1 

IF( IX.NE.O)  GO  TO  111 

IX  =  1 

GO  TO  1055 
125  IF(K.EO»l)  GO  T3  1055 

J  =  l 

I X  =  I X + 1 

GO  TO  111 
C   SET  INDEXES  FO'J.  INTERPOLATION 
150  IF(  IX.LT.l  )  IX  =  1 

J=J4-J5+I 

IF« IX.GT.J )  IX=J 

L=IX+J5-1 
C   MAKE  INTERPOLATION 

Y=0. 

DO  4  0  0  K=IX»L 

Y  L = 1  •  0 

DO  3  00  J=IX»L 

IF(J.EQ.K)  GO  TO  3  00 

XJ=XFNC( I »J»1) 

YL=YL*  !  X-X  J  )  /  (  XF^IC  (  I  ,'<  *  1  )  -XJ  ) 
3  00  CONTINUE 
400  Y=Y+YL*XFNC< I » < *2 > 

RETURN 

END 
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APPENDIX  C 
Baseline  Case:   Data  Input  and  Sample  Output 

A  list  of  all  baseline  variables  and  their  definitions 
is  presented  in  the  following  table.   The  variables  are 
listed  according  to  their  function  and  subroutine.   After  the 
definition  of  each  variable  the  proper  symbolic  name  used  in 
the  computer  program  is  shown  in  parenthesis. 

Since  the  model  presently  contains  a  simple  demographic 
sub-model  the  major  portion  of  the  demographic  variables  are 
not  being  utilized.   Also,  in  the  baseline  case,  employment/ 
output  ratios  were  held  constant  after  1972. 


SIMULATION  CONTROL  VARIABLES 

TSTRT  -  Starting  data  for  simulation  (yr) 
TEND   -  Ending  date  for  simulation  (yr) 
DT     -  Time  step  (yr) 


Scenario  Input 


ISOER  -  Employee/output  ratios 

If  ISOER  equals  1,  use  E/o  ratios  constant 
after  1972. 

ISCON  -  Residential  conservation 

If  ISCON  equals  1,  use  residential  con- 
servation subroutine. 

ISPRC  -  Firm  industrial  substitution 

If  ISPRC  equals  1,  use  firm  industrial  sub- 
stitutions . 

ISPRM  -  Estrapolated  industrial  substitution  of 
ISPRM  equals  1,  use  extrapolation. 


Print 

Control 

IPRNT 

IPRNT2 

1 
0 
1 
0 

0 
1 
1 
0 

Computer  Printout 

input  and  output 

output  of  demands 

input  and  output  of  demands 

input/output  variables 

C-l 


AGRICULTURAL  PARAMETERS 


Pl 
P2 


?7 


P10 
Pll 


r12 
P13 


Price  of  agricultural  land  (PL) 
Agricultural  operating  expense/out- 
put (OP) 
Agricultural  output/labor  (OLP) 
Property  tax  rate  (FTR) 
Agricultural  interest  rate  (RA) 
Fraction  of  agricultural  land 

mortgaged  (PLF) 
Propensity  to  invest  (PTI) 
Fractional  loss  before  disinvest- 
ment (TF) 
Salvage  value  fraction  (SALV) 
Leverage  on  cash  (ALEVI) 
Agricultural  depreciation  rate 

(AKDEP) 
Capital  weighting  constant  (ALP) 
Constant  relating  Agricultural  out- 
put to  capital  &  land  (A) 


$100.00/(acre) 

.675 

$58,500.00/(man/vr) 

.0096-1% 

.06-6% 

.06-6% 
.2 

.2 
.3 
1.0 

.077-7.7% 
.2 

10.126 


FORESTRY    PARAMETERS 


*'14 
Pl5 

p16 
p17 


Allowed  cutting  rate  (AC) 
Capital/cutting  rate  (FCR) 
Depreciation  rate  forestry  (DR) 
Investment  closure  time  (CT) 


9  0.0  (BF/acre/yr) 

.053 

.067-6.7% 

5.0  (yrs) 


GAS  AND  OIL  PARAMETERS 


F18 

P19 

?20 

p21 
P22 

?23 

P24 
P25 
P104 


Gas  produced  with  oil  for  oil  fields 

(GOR) 
Oil  produced  with  gas  for  gas  fields 

(OGR) 
Capital  cost  for  production  invest- 
ment (CCP) 
Drilling  cost  (DC) 
Gas  discovered  in  oil  expenditures 

(GDOE) 
Oil  discovered  in  gas  expenditures 

(ODGE) 
Exploration  investment  life  (DPT) 
Capital  cost  for  exploration  (CC) 
Discount  rate 


400.0  (ftVbbl) 

.000010  (bbl/ft3) 

.16  ($/$  capital/yr) 
5.0  ($/ft) 

0  (ft3/bbl) 

0  (bbl/ft3) 

5.0  (yr) 

.3  ($/$  capital/yr) 

.1 


COAL  MINING  PARAMETERS 


'26 
'93 


94 


103 


Coal  mining  capital  cost  (CC)  .16 
Price  sensitivity  to  saturation  of 

coal  market  .10 
Price  sensitivity  to  unsaturation 

of  coal  market  .30 

Rate  new  firms  enter  coal  market  .2 

C-2 


n 


■ 


■ 
■ 


COAL  MINING  PARAMETERS  (cont'd) 


IP 


27 


I?28 
IP29 

IP96 


Coal  mining  cost  function  (Kl)  #12 
Coal  mining  capital  function  (K2)  #13 
Coal  mining  cost  distribution  func- 
tion (K3)  #14 
Coal  cost  distribution  function  #30 


MINING  PARAMETERS 


P30 
IP31 

IP39 
IP33 


IP 


97 


Mining  capital  cost  (CC) 
Mining  cost  function  (Kl) 
Mining  capital  function  (K2) 
Mining  cost  distribution  function 

(K3) 
Mine  cost  distribution  function 


.16 

#15 
#16 

#17 
#31 


SUPPLY  AND  DEMAND  PARAMETERS 

P34     Primary  metals  output/mining  out- 
put (PMMR)  2.4  8 
p35     Fraction  of  oil  available  for  mining 

(FRMA)  .29 

P36     Base  oil  price  (BOP)  $4.00/bbl 


DEMAND  PARAMETERS 


.37 


39 


'40 
'l02 


IP 

IP 
PA 


42 
43 

66 


Fraction  of  earned  income  spent 
State  government  expenditures/house- 
hold expenditures  (RGPES) 
Federal  government  expenditures/ 
household  expenditures  (RGPEF) 
Rail  cost/ton  (RCC) 
Rail  export  constant  (REXC) 
Federal  government  expenditures/ 

person 
Rail  sector  index  (IRAIL) 
No.  of  supply  controlled  sectors  (KS) 
Household  expenditure  vector  (PED[25]) 


0.0 

.00477 

.01597 

.01041 

0.0 


03457 
00044 
02840 
06142 


01938 
00108 
04275 
00113 


00413 
04046 
16852 
00110 


.03211 
.00446 
.18643 
.00262 


.715-71.5% 

.1215 

.1013 

$.01/ (ton/mile) 

$100,000,000 


$786 

11 

23 


00 


.00276 
.01632 
.04549 
0.0 


PAg7    Exports/household  expenditure  vector  (EX [25]) 


0.0 

0.0 

.41430 

0.0 

0.0 


.57493 
0.0 

.73936 
0.0 


0.0 
0.0 
0.0 
0.0 


0.0 
0.0 

.62587 
0.0 


0.0 
0.0 
0.0 
0.0 


0.0 
0.0 
0.0 
0.0 


C-3 


n 


PA 


68 


State  government  expenditure  vector  (SED[25]) 


.22178 

.00217 
.00545 
.00834 
0.0 


00202 
00009 
01810 
02711 


00093 
00070 
00634 
00118 


00010 
00877 
00371 
00015 


00115 
00073 
01100 
01304 


.00321 
.00907 
.00123 
0.0 


PA 


69 


Federal  government  expenditure  vector  (FED  [25]) 


.04042 
.00659 
.00599 
.01306 
0.0 


00140 
00123 
00379 
02058 


.00387 
.00046 
.01228 
.00140 


.00042 

.10503 

-.00024 

.00184 


00066 
00420 
00363 
00141 


.00083 
.02819 
.01169 
0.0 


OUTPUT  PARAMETERS 


IP 

PA 


44 
70 


Number  of  demand  sectors  7 

Direct  coefficients  supply  output-demand  (ADS  [25,10]) 


.00652 

.01121 

0.0 

.00005 

0.0 

0.0 

.02365 

0.0 

.00314 

.00062 

.00007 

.00004 

.01444 

.00832 

.00023 

.00071 

.00051 

.00015 

.00287 

.00296 

.00543 

.00161 

.00935 

.04367 

.00237 

.00085 

.00645 

.01826 

.01977 

.00612 

.01485 

.00039 

.05611 


.00435 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

.00001 

0.0 

.00004 

0.0 

.00382 

0.0 

.00297 

0.0 

.00642 

0.0 

.00643 

0.0 

.01740 

0.0 

.00720 

0.0 

.00276 

0.0 

.01176 

0.0 

.02388 

0.0 

.00232 

0.0 

.01570 


.02582 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

.00001 

0.0 

.00282 

0.0 

.00142 

0.0 

.00568 

0.0 

.01093 

0.0 

.00169 

0.0 

.00170 

0.0 

.00090 

0.0 

.01043 

0.0 

.00442 

0.0 

.00395 

0.0 

.17407 


.00549 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

.00001 

0.0 

.00473 

0.0 

.00312 

0.0 

.00324 

0.0 

.01760 

0.0 

.00500 

0.0 

.00319 

0.0 

.00161 

0.0 

.02175 

0.0 

.01553 

0.0 

.00171 

0.0 

.05778 


.00749 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

.00004 

0.0 

.01078 

0.0 

.00068 

0.0 

.00250 

0.0 

.02012 

0.0 

.00648 

0.0 

.00315 

0.0 

.00769 

0.0 

.03140 

0.0 

.01066 

0.0 

.00308 

0.0 

.04581 


.00434 

0.0 

0.0 

0.0 

0.0 

0.0 

.00117 

0.0 

0.0 

0.0 

.00002 

0.0 

.00595 

0.0 

.00015 

0.0 

.00158 

0.0 

.00418 

0.0 

.01910 

0.0 

.00622 

0.0 

.00142 

0.0 

.05076 

0.0 

.01553 

0.0 

.00065- 

0.0 

.01010 


C-4 


n 

ii 
m 


m 
m 


PA70    Direct  coefficients  supply  output-demand 
(ADS [25])  cont'd 


.03410 

.00038 

.00060 

.01168 

.01677 

.00266 

.00099 

.00003 

.00029 

.00002 

.00001 

.00091 

.00179 

0.0 

0.0 

0.0 

0.0 


0.0 

.00031 

0.0 

.01069 

0.0 

.00226 

0.0 

.00030 

0.0 

.00013 

0.0 

.00309 

0.0 

0.0 

0.0 

0.0 


0.0 
.00030 

0.0 

.00549 

0.0 

.00759 

0.0 

.00021 

0.0 

0.0 

0.0 

.01303 

0.0 

0.0 

0.0 

0.0 
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LABOR  PARAMETERS 


^45 
P46 
P47 
p48 
"49 


IP 
PA 


50 
72 


Labor/population  (WPR)  .4 

Derivative  wages  sensitivity  (USD)  1.0 

Primary  wage  sensitivity  (USP)  1.0 

Hired  agriculutral  wages  (USA)  1.0 

Fraction  of  farm  labor  hired  (FLH)  .3 

Total  number  of  sectors  (KT)  30 
Employment/out 


000033 
000028 
000032 
000118 
000030 


000016 
000028 
000014 
000118 
000014 


.000016 

.000028 

.000078 

0.0 

.000011 


.000016 
.000028 

.000078 

0.0 

.000025 


000016 
000037 
000017 
000118 
000009 


.000028 

.000038 

.000118 

.0000407 

.000005 


I  PA 


73 


Sector  employment  types  (ILB[35]) 

1- agricultural 
2-primary  employment 
3-derivative  employment 


3 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

3 

3 

3 

3 

3 

3 

3 

3 

3 

2 

3 

3 

1 

2 

2 

2 

2 

2 

2 

CAPITAL  INVESTMENT  PARAMETERS 


Delay  times  (ATD) 

Sector  depreciation  rates 


.093 
.067 
.033 
.097 
.065 
0.0 


065 
082 
033 
097 
067 


5.0  (yrs) 


rates 

(DPR[35] ) 

065 

.065 

.065 

082 

.082 

.038 

052 

.052 

.047 

097 

.097 

.097 

067 

.067 

.082 

066 
038 
097 
065 
067 


PA 


75 


Sector  capital-output  ratios  (CAR[35]) 


.19 

.5 

3. 

.95 

.25 

0.0 


25 

59 
!. 
95 


.25 
.59 
.65 
.95 

.5 


.25 
.59 
.65 
.95 
.5 


.25 

.39 

1.76 

1.763 

1.05 

.95 

.95 

.25 

.59 

.5 

I  PA 


76 


Investment  calculation  index  (IN[35]) 


1 

1 
0 


1 
1 
0 


1 
1 

0 


1 

1 

0 


1 
1 
1 


1 
1 

1 


1 
1 


1 

1 


1 
1 


1 
1 


1 
1 


1 

0 


C-8 


■ 

n 


■ 


PA 
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Capital  expenditure  vector  (CPS[25]) 


.55070 

0.0 

.01408 

0.0 

0.0 


0.0 

.00082 
.03002 
0.0 


0.0 
0.0 

.05412 
0.0 


0.0 

.03365 
.05400 
0.0 


0.0 

.01502 
.00591 
0.0 


0.0 

.00810 
0.0 
0.0 


ENERGY  DEMAND  PARAMETERS 


IP52    No.  of  energy  types  (NENR)  (electri- 


IP 


53 


PA 


city,  natural  gas,  oil,  coal) 

No.  of  end  uses  (space  heating, 
kitchen  range,  refrigerator,  freezer, 
clothes  dryer,  air  conditioner,  other)  8 


77    Appliance  energy  use  (AR[r,10]) 


27400. 

460. 

6420. 

0.0 

0.0 

161900000 

0.0 

6420000. 

0.0 


4220. 

740. 

0.0 

0.0 

0.0 

88300000 

0.0 

0.0 


1180. 

0.0 

0.0 

141300000. 

0.0 

0.0 

0.0 

0.0 


1140. 

0.0 

3250. 

26500000 

0.0 

0.0 

0.0 

3250000. 


1200. 

93800. 

0.0 

0.0 

0.0 

0.0 

93800000 

0.0 


990. 

16300. 

0.0 

0.0 

0.0 

0.0 

16300000 

0.0 


PA 


78 


Appliance  saturations  (SAT[5,10]) 


- 


.023 
.134 
.193 

0.0 

0.0 

.019 

0.0 

.023 

0.0 


.179 

1. 

0.0 

0.0 

0.0 

0.0 

1. 

0.0 


.782 

0.0 

0.0 

.071 

0.0 

0.0 

0.0 

0.0 


1.051 

0.0 

.039 

0.0 

1. 

0.0 

0.0 

.003 


.653 
.833 

0.0 

0.0 

0.0 

0.0 

.049 

0.0 


.688 

.787 

1. 

0.0 

0.0 

0.0 

.029 

1. 


PA 


79 


Energy  coefficients  (AI[5,35]) 


^ 


.859 

19.911 

.859 

.859 

1.292 

0.0 

13.55 

31.4 

5.22 


.270 

.0417 

0.0 

.859 

.786 

0.0 

13.55 

103.8 

6.76 


.270 

.927 

.859 

.859 

.0768 

0.0 

13.55 

6.02 

6.76 


.270 

.0068 

.859 

.859 

.137 

0.0 

13.55 

0.0 

6.76 


.270 

.129 

.859 

.859 

8.708 

0.0 

15.7 

0.0 

6.76 


.199 

0.0 

.859 

.840 

1.335 

6.76 

44.7 

6.76 

6.76 


C-9 


■ 

■ 
■ 
■ 
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PA 


79    Energy  coefficients  (AT [5,35])  cont 


6.76 

6.76 

6.76 

6.76 

0.0 

18.73 

0.0 

6.43 

6.48 

29.8 

19.18 

0.0 

0.0 

0.0 

0.0 

0.0 

4872. 

7006. 

7006. 

7006. 

7006. 

7006. 

7006. 

7006. 

7006. 

7006. 

60131. 

6013.1. 

4872. 

1098. 

4872. 

4872. 

4872. 

4872. 

4872. 

4872. 

4872. 

4872. 

4872. 

31292. 

7006. 

7006. 

7006. 

7006. 

7006. 

7006. 

0.0 

0.0 

0.0 

0.0 

0.0 

1230. 

2340. 

2340. 

2340. 

2340. 

2340. 

2340. 

2340. 

2340. 

2340. 

2340. 

2340. 

1230. 

53440. 

1230. 

1230. 

1230. 

1230. 

1230. 

1230. 

1230. 

1230. 

1230. 

0.0 

2340. 

2340. 

2340. 

2340. 

2340. 

2340. 

0.0 

0.0 

0.0 

0.0 

0.0 

1466.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1466. 

0.0 

1466. 

1466. 

1466. 

1466. 

1466. 

1466. 

1466. 

1466. 

1466. 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

DEMOGRAPHY  PARAMETERS 


p54 
p55 
IP56 


IP 

ip; 
ip 


57 


59 


60 


IP 


PA 


61 


80 


Unemployment  migration  constant  (UMC)  .3 

Wage  migration  constant  (WMC)  .5 

Age  children  begin  leaving  home  (IAMl)  15 

Age  all  children  have  left  home  (IAM2)  24 

Age  family  formation  starts  (IAF1)  16 

Age  family  formation  ends  (IAF2)  36 
Single  family  females/married  females 

(SMR)  .086 
Marrying  age  difference  between  fe- 
males &  males  (MAD)  2  (yrs) 
Relative  age  mobilities  (RM[4]) 


1.0 


1.0 


PA 


81 


Vlatura 

tion  ra 

tes 

(AMT [25] ) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

.05 

.10 

.25 

.40 

.55 

.65 

.75 

.85 

.95 

1.00 

C-10 


- 


PA. 


82 


PA< 


PA 


84 


■ 

■ 


Marriage  rates  (AMR 

40]) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

.039 

.078 

.117 

.156 

.135 

.105 

.075 

.045 

.020 

.015 

.015 

.015 

.015 

.015 

.0035 

.0035 

.0035 

.0035 

.0035 

0.0 

0.0 

0.0 

0.0 

0.0 

O.C 

Death 

rates  (DH[30]) 

0.0 

.034 

.069 

.103 

.138 

.172 

.206 

.241 

.275 

.310 

.344 

.378 

.413 

.241 

.275 

.310 

.344 

.378 

.619 

.654 

.688 

.722 

.757 

.791 

.826 

.860 

.894 

.929 

.963 

1.0 

Birth 

rates  (BH[6,45]) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

.0233 

0.0 

0.0 

0.0 

.0466 

.0103 

0.0 

0.0 

.0699 

.0206 

0.0 

0.0 

.0932 

.0309 

.0048 

0.0 

.1164 

.0413 

.0097 

0.0 

.1050 

.0516 

.0145 

.0022 

.0815 

.0575 

.0242 

.0007 

.0018 

.0005 

.0699 

.0531 

.0291 

.0090 

.0026 

.0010 

.0582 

.0486 

.0339 

.0112 

.0035 

.0014 

.0466 

.0442 

.0317 

.0134 

.0044 

.0019 

.0349 

.0398 

.0274 

.0157 

.0053 

.0024 

.0233 

.0354 

.0271 

.0148 

.0061 

.0029 

.0116 

.0309 

.0249 

.0139 

.0070 

.0033 

0.0 

.0265 

.0226 

.0130 

.0079 

.0038 

0.0 

.0221 

.0204 

.0122 

.0074 

.0035 

0.0 

.0177 

.0181 

.0133 

.0068 

.0033 

0.0 

.0133 

.0158 

.0105 

.0063 

.0030 

0.0 

.0088 

.0136 

.0096 

.0058 

.0028 

0.0 

.0044 

.0113 

.0087 

.0053 

.0025 

C-ll 


PA 


84 


Birthrates  (BH[6,45J)  cont'd 


IP 
IP 


91 
92 


95 


0.0 
0.0 
0.0 
0.0 

0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 


Age  females  start 
Age  females  stop  g 
Number  of  households  per  elderly 
person 


.0090 

.0078 

.0047 

.0023 

.0068 

.0070 

.0042 

.0020 

.0045 

.0061 

.0037 

.0018 

.0023 

.0052 

.0032 

.0015 

0.0 

.0044 

.0026 

.0013 

0.0 

.0035 

.0021 

.0010 

0.0 

.0026 

.0016 

.0008 

0.0 

.0017 

.001.1 

.0005 

0.0 

.0009 

.0005 

.0003 

0.0 

0.0 

0.0 

0.0 

ving  birth 

(IB1) 

16 

ing  birth 

(IB2) 

45 

LAND  USE  PARAMETERS 


,62 
63 

85 


Tons    of    coal   mined/acre     (TONACR) 
Fraction   of    strip-mined   land    revege- 

tated    (REVEG) 
No.    of    households/acre    (CD) 
No.    of   years    to    revegetate 
No.  of  employees/acre    (EMPACR) 


28320     (tons/acre) 

.8-80% 

.354    (h/acre) 

10     (yrs) 

1.3    (emp/acre) 


HOUSEHOLD  CONSERVATION  PARAMETERS 


H 


106 


107 


108 


109 


110 


111 


112 


113 


Energy  use  in  an  average  home  heated 

with  no  insulation  (BASE) 
Fraction  of  energy  use  in  a  class  "B" 

insulated  home  (USB) 
Fraction  of  energy  use  in  a  class  "C" 

insulated  home  (USC) 
Cost  of  insulating  an  existing  home  to 

class  "B"  standards  (CBO) 
Cost  of  insulating  an  existing  home  to 

class  "C"  standards  (CCO) 
Cost  of  insulating  a  new  house  to 

class  "B"  standards  (CBN) 
Cost  of  insulating  a  new  house  to 

class  "C"  standards  (CCN) 
Life  span  of  an  average  house  (HRR) 


135840.  (ft3) 

.74 

.315 

$360 

$2,470. 

$300. 

$1,970. 
75  (yrs) 


OTHER  PARAMETERS 


;s6 

'87 


Forestry  output/timber  cut 
Base-gas  price 


003  ($/ftJ) 


C-12 


n 


OTHER  PARAMETERS  (cont'd) 


'88 
'89 


Base  coal  price 

Base  raining  product  price 


3.4  ($/ton) 
1.0  ($/ton) 


AGRICULTURAL  STATE  VARIABLES 


•15 
;16 


Owned  agricultural  capital  (AKO) 
Financed  agricultural  capital  (AKF) 


$210,000,000 
$210,000,000 


FORESTRY  STATE  VARIABLES 


17 


Timber  allowed  to  cut  (TAC) 
Forestry  capital  investments  (CAP) 


1,500,000,000 
$80,000,000 


(BF/yr) 


GAS  AND  OIL  STATE  VARIABLES 


S,q     Abandoned 
S20     Abondoned 

oil 

fields  (AFO) 

0 

gas 

fields  (AFG) 

0 

Sot     Exploration  investment  for  oil 

(ODI) 

$10,000,000 

S       Exploration  investment  for  gas 

(GDI) 

$5,000,000 

sf             Producing 

oil 

Fields  (PRO[50]) 

20,000,000   0.0 

30 

,000,000   0.0 

10,000 

,000   0.0 

50,000,000   0.0 

90 

,000,000   0.0   100,000 

,000   0.0 

100,000,000   0.0 

125 

,000,000   0.0   125,000 

,000   0.0 

125,000,000   0.0 

125 

,000,000   0.0   125,000 

,000   0.0 

100,000,000   0.0 

95 

,000,000   0.0    i 

90,000 

,000   0.0 

85,000,000   0.0 

80 

,000,000   0.0 

75,000 

,000   0.0 

70,000,000   0.0 

65 

,000,000   0.0    1 

50,000 

000   0.0 

55,000,000   0.0 

So      Producing 

gas 

fields  (PRG[50]) 

80,000,000,000 

0.0 

100,000,000,000 

0.0 

130,000,000,000 

0 

0 

130,000,000,000 

0.0 

150,000,000,000 

0.0 

150,000,000,000 

0 

0 

150,000,000,000 

0.0 

150,000,000,000 

0.0 

130,000,000,000 

0 

0 

120,000,000,000 

0.0 

90,000,000,000 

0.0 

80,000,000,000 

0 

0 

70,000,000,000 

0.0 

60,000,000,000 

0.0 

50,000,000,000 

0 

0 

40,000,000,000 

0.0 

30,000,000,000 

0.0 

30,000,000,000 

0 

.0 

25,000,000,000 

0.0 

25,000,000,000 

0.0 

20,000,000,000 

0 

0 

20,000,000,000 

0.0 

20,000,000,000 

0.0 

15,000,000,000 

0 

0 

15,000,000,000 

0.0 

COAL  MINING  STATE  VARIABLES 


!23 
'37 


Total  coal  mined  (TM) 
Actual  coal  price  (PRC) 


0 
3 


(tons) 
($/ton) 


C-13 


n 


■ 


Coal  investment  (A1N[20]) 


3,150,000 
2,520,000 
1,890,000 
1,260,000 


3,050,000 
2,420,000 
1,780,000 
1,150,000 


2,940,000 
2,310,000 
1,680,000 


2,840,000 
2,210,000 
1,580,000 


2,730,000 
2,100,000 
1,470,000 


2,630,000 
1,990,000 
1,360,000 


MINING  STATE  VARIABLES 


^24 
S6 

30,000 
24,000 
18,000 


Total  mined  (TM) 
Mining  investment 


000 
000 
000 


12,000,000 


29,000,000 
23,000,000 
17,000,000 
11,000,000 


:  (A1N[20]) 

28,000,000 
22,000,000 
16,000,000 


27,000,000 
21,000,000 
15,000,000 


0.0 (tons) 


26,000,000 
20,000,000 
14,000,000 


0 
0 
0 


SUPPLY  AND  DEMAND  VARIABLES 

S25     Oil  refinery  capacity  (PRROD) 


59,000,000  (bbl/yr) 


OUTPUT  STATE  VARIABLES 

Soo     Average  output  for  Montana  (AVOUT) 


$5,390,000,000 


LABOR  STATE  VARIABLES 


^26 
b27 

^28 
S39 


Derivative  wage  rate  (DW) 
Primary  wage  rate  (PW) 
Hired  agricultural  wage  rate  (HAL) 
Average  household  income  (AVHHI) 


5,072.00  ($/yr) 
9,480.00  ($/yr) 
3,340.00  ($/yr) 
10,500.00  ($/yr) 


CAPITAL  INVESTMENT  STATE  VARIABLES 

S-34     Averaged  rate  of  change  (RC[35]) 


'35 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

Averaged  outputs   jr  each  sector  (AO[35]) 


1. 

Construction 

398,000,000 

6 

2. 

Meat  Process 

106,000,000 

7 

3. 

Dairy  Process 

44,000,000 

8 

4. 

Grain  Process 

40,000,000 

9 

5. 

Other  Food 

67,000,000 

10 

C-14 


Print  &  Publish 
Chemical 

Fabricated  Metals 
Stone  &  Clay 
Other  Manu- 
facturing 


32,000,000 
25,000,000 
14,000,000 
42,000,000 

225,000,000 


Averaged  outputs  for  each  sector  (AO[35j)  cont'd 


35 

11.  Railroad  Trans- 
portation 

12.  Motor  Freight 

13.  Communications 

14.  Utility 

15.  Wholesale 

16.  Retail  Trade 

17.  Finance,  Insurance 
&  Real  Estate 

18.  Personal  Service 

19.  Business  Repair 

20.  Professional  Ser- 
vices 


179,000,000   21 


135 

108, 

213 

173 

513 

533 

119 

86 


,000,000 
,000,000 
,000,000 
,000,000 
,000,000 

000,000 
000,000 
000,000 


185,000,000 


22 

23 
24 
25 
26 

27 
28 

29 
30 


Federal  Enter- 
prizes 

State  &  Local 
Enterprizes 
Dummy 

Agricultural 
Forestry 
Petroleum  & 
Natural  Gas 
Coal  Mining 
Non-Energy 
Mining 

Primary  Metals 
Petroleum  Re- 
fining 


18,000,000 

29,000,000 

39,000,000 

931,000,000 

299,000,000 

108,000,000 
17,000,000 

177,000,000 
440,000,000 

233,000,000 


DEMOGRAPHY  STATE  VARIABLES 


s7     I 

un  child  at 
3278. 

age  J  (PC[6,25] ) 

1756.         1054. 

463. 

4683. 

234. 

4488. 

314  2. 

19  0  8. 

1010. 

449. 

224. 

4336. 

3035. 

1843. 

976. 

434, 

217. 

4577. 

3204. 

19  4  5. 

1030. 

458. 

229. 

4746. 

3322. 

2  017. 

10  6  3. 

475. 

237. 

5196. 

3637. 

2203. 

n.69_. 

1258. 

.  §.2Q_._ 
559. 

260. 
280. 

5592. 

3914. 

2377. 

5926. 

414  9. 

2519. 

1333. 

5  93. 

296. 

6007. 

^235. 

2553. 

1352. 

1319. 

601. 

300. 

618?. 

4327. 

2627. 

618. 

309. 

6486. 

4  2  60. 

2  5  8  7. 

1369. 

609. 

304. 

6182. 

4327. 

2627. 

1319. 

618. 

309. 

6381  . 

4467. 

2712. 

1436. 

533. 

319. 

6221  . 

4355. 

2644. 

1400. 

622. 

311  . 

6177. 

4  3  24. 

2625. 

1390. 

618. 

309. 
302. 

6172. 

<^320. 

2623. 

1389. 

617. 

5  59S. 

3988. 

2421. 

1282. 

570. 

285. 

5082. 

3557. 

2160. 

1143. 

_  50^. 
345. 

254. 
173. 

3449. 

24  14. 

1466- 

776. 

18  2  2. 

1275. 

774. 

410. 

132. 

91- 

76?  . 

538. 

32  6. 

173. 

77. 

39. 

262. 

184. 

112. 

59. 

26. 

13. 

B7. 

61. 

37. 

20. 

9. 

5. 

8. 

5. 

3. 

2. 

1. 

.5 

0 

0 

0 

0 

0 

0 

Number  of  married  females  age  J  (PFF[60]) 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

0 

0 

0 

33. 

962. 

1633. 

2473. 

3251  . 

3740. 

4205. 

4159. 

3417. 

3705. 

3759. 

4023. 

3603. 

3635. 

C-15 


8 

3686 
3210 
3602 
3786 
3420 


Number  of  married  females  age  J  (PFF@60])  cont'd 


3467 
3315 
3502 
3505 
3448, 


3418 
3400 
3550 
3838 
3311, 


3525. 
3466. 
3528. 
3737. 
3224. 


3418 

3775 

3612, 

3633 

3146, 


3486 
3583 
3547 
3591 
3019 


Number  of  young  single  females  age  J  (PYF[60]) 


0.0 

0.0 

0.0 

0 

.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 

.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 

.0 

0.0 

72. 

742. 

1089. 

1437. 

1399. 

1140. 

937. 

669. 

614. 

567. 

494. 

395. 

303. 

30  7. 

288. 

241. 

249. 

241. 

247. 

0.0 

0.0 

0.0 

0 

.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 

.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 

0 

0.0 

0.0 

0.0 

S10 

Number  of 

mature  single  femal 

es  age  J 

(PSF[60] ) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Q 

0._. 

0  .  . 



0    .  - 

c 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

0 

c 

0 

0 

0__ 

0  _ 

0  - 

o 

o 

242. 

250. 

216. 

2  61. 

224. 

27C. 

271. 

264. 

267. 

165. 

272. 

267. 

285. 

2  64. 

269. 

381  . 

273  . 

270 . 

257. 

250. 

249. 

2^2. 

237. 

227. 

Sll 

Number  of 

mature  single 

s  males 

age  J  (PSM[6  0]  ) 

0 
0 

0 
0 

0 
0 

0 

0 

0 
n 

0 
0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4R3  . 

49?. 

512. 

522. 

563. 

539. 

542. 

527. 

53^. 

531. 

544. 

-  534. 

570. 

532. 

57S. 

563. 

54  7. 

541. 

515. 

519. 

493. 

455. 

474. 

454. 

S12     Number  of  young  single  males  age  J  (PYM[60]) 


0 
0 
875. 


0 
0 
1285. 


0 
0 
0 
16  94. 


0 
0 
0 
1965 


0 
0 

36. 
1473. 


143. 
1265. 


C-16 


i 


"12 


Number  of  young  single  males  age  J  (PYM[60])  cont'd 


i 


940. 
599. 

0.0 
0.0 
0.0 
0.0 

S13 

907.        865. 
563.        512. 
0.0         0.0 
0.0         0.0 
0.0         0.0 

Number  of  married  males 

813.        681. 
528.        512. 
0.0         0.0 
0.0         0.0 
0.0         0.0 

age  J  (PFM[60] ) 

591. 

523. 

0.0 

0.0 

0.0 

0 
0 
0 

0 
0 
0 

0 
0 
0 

0 
0 

0 

0 
0 
2. 

0 
0 
3  86. 

15  5  0. 
3146. 

3  3  9'-. 

227. 

3412. 

3192. 

2994. 

3461. 
3147. 

34444. 

3704. 

3246. 

3205. 

3263. 

345  5. 

3  8  7  2. 
3322. 
314  7. 

3  3  3  0. 
3347. 
3  210. 

2969. 
3  3  31. 

3501. 

3  0  66. 
3239. 
3237. 

3144. 
3233. 
3549. 
3062. 

derly  age  J 

5772. 
38^6. 
3135. 
2198. 
885. 

3491. 
3340. 
3359. 

3314. 
3230. 
3320. 

3162. 

S 
14 

6538. 
o,  v  5  2  . 

3189. 

Number  of  el 

5810. 
4374* 

2  9  5  2. 

(ELD[30] ) 

5265. 
4076. 

2  909. 

5050. 
3350. 

2792. 

4940. 
3399. 

32C7. 
265^. 
885. 

3207. 
2323. 
885. 

3039. 
1911  . 
325. 

2963. 
1592. 
885. 

2551. 
1524. 
253. 

LAND  USE  STATE  VARIABLES 


|29 

;30 
!31 
|32 
|33 
'36 


Forest  land  (FL) 

Agricultural  land  (AGL) 

Waste  land  (UNL) 

Urban  land  (URL) 

Wilderness  land  (PRL) 

Mined  land  (third  order  delay) 


(AML[3J) 


16,660,000 

62,900,000 

10,562,000 

818,000 

2,145,000 


(acres) 
(acres) 
(acres) 
(acres) 
(acres) 


1,333.3 


1,333.3 


1,333.3 


HOUSEHOLD  CONSERVATION  STATE  VARIABLES 


140 


'41 

!42 

'43 


65100. 
86800. 
65100. 
.0000014 


C-17 


MISCELLANEOUS  VARIABLES 


AI. 


Demands  from  capital  investment  (CAP [25]) 


.000156 

.000004 

0.0 

0.0 

AI, 


0.0 
.00000023 


0.0 
.0000085 


0.0 

.0000095 

.0000153 


0.0 

.0000043 

.0000017 


Lagged  outputs  for  each  sector  (ODS[35]) 


0.0 
.0000023 


398000000 
25000000 

108000000 
86000000 

299000000 


106000000 
14000000 
213000000 
185000000 
108000000 


44000000 
42000000 
173000000 
18000000 
17000000 


40000000 
255000000 
513000000 

29000000 
177000000 


67000000 
179000000 
533000000 

39000000 
440000000 


32000000 
135000000 
119000000 
931000000 
233000000 


A 
A 
A 
V 

vJ 


26 
27 

30 

40 


Earned  income  (AIN)  or  (TI) 

Total  Montana  population  (PT) 

Total  Montana  households 

Limit  on  tons  of  coal  mined  per  year 

Previous  total  employment 

Previous  total  households 


$2,400,000,000.00 
694,000 


217,000 
300,000,000 
265,700 
217,000 


(tons) 


FUNCTIONS 

The  follov;ing  functions  are  used  as  a  means  to  inter- 
polate between  given  points.   AFNC  (X,Y,Z)  is  the  notation 
used  in  the  program  where  X  is  the  value  of  the  independent 
variable,  Y  is  the  value  of  the  dependent  variable  calculated 
and  Z  is  the  function  number. 

AFNC  (TIME,  A2,20)  -  Price  of  agricultural  products  (P) 


1970. 


1972, 

1. 


1974 
1.5 


1976, 
1.3 


1978, 
1.3 


1980, 
1.3 


AFNC  (TIME,  A   ,21)  -  Oil  price  (OP) 


1970, 
2000, 
15.0 


1975. 

1980. 

1985. 

1990. 

1995 

4.0 

7.0 

10.0 

11.0 

12.0 

20.0 

AFNC  (TIME,  A1]L,22)  -  Gas  price  (GP) 


1970. 

1975. 

1980. 

1985. 

1990. 

1995 

2000. 

.0002 

.0004 

.001 

.002 

.003 

.004 

.035 

C-li 


II 

n 


■ 


AFNC 


(TIME,    A15,23)    -    Coal    price    (PRC) 


1970. 

1975. 

1980. 

19S5. 

1990. 

1995 

2000. 

3.0 

4.0 

5.0 

7.0 

9.0 

11.0 

15.0 

« 


7AFNC 


(TIME,  A   ,24)  -  Price  of  mined  products 


12 


1970. 

1972. 

1974. 

1976. 

1978. 

1980. 

1982. 

1984. 

1986. 

1988. 

1990. 

1992. 

1.0 

.9 

1.0 

1.1 

1.2 

1.25 

1.3 

1.35 

1.4 

1.45 

1.5 

1.5 

AFNC     (TIME,  A25, 25)- Oil    imports    available    (AIMP) 


1970. 
25000000, 


1980. 
20000000, 


1990. 


2000, 


50000000. 


50000000. 


AFNC     (TIME,    A       ,  26)-National    unemployment    (UEN) 


1970. 
.07 


1975. 
.07 


1980. 


1985. 


.05 


.08 


AFNC    (TIME,    A       , 2  7) -National   wage   rate     (WRN] 


1970. 
6800. 


1975. 
7000. 


1980. 


1990. 


6500. 


6700. 


AFNC 


(TIME,    A35,28)     -    Reforestation    rate    (REFOR) 


1970. 


1980. 


1990. 


0.0 


0.0 


0.0 


AFNC 


(TIME,    A3?,29)     -    Preservation   rate     (PRESER) 


1970. 


1980. 


1990, 


0.0 


0.0 


AFNC 


(TIME,    A.     ,32)    -    Limit   on    tons    of    coal   mined 


1970. 
25000000. 


1975. 
35000000. 


16  2 

1980.  1985.  1990. 

100000000.  150000000.  200000000. 


0.0 


1995. 
200000000. 


AFNC    (TIME,    A, 33)    -    Employee/output   ratios    for   non-energy  mining 

Sector    (28) 


0 


1972, 


1980. 


1985, 


.000025 


.000018 


.000016 


- 


C-19 


i 
■ 


* 


• 


AFNC     (TIME,    A, 34) 


-    Employee/output   ratios    for   coal   mining 
Sector    (27) 


0 


T  1972. 

1980. 

1 

AFNC 

(TIME, 

A, 35) 

1  1972. 

i 

1 
1980. 

J 

AFNC 

(TIME, 

A, 36) 

1  1972. 

i 

1 
19S0. 

AFNC 

(TIME, 

A, 37) 

^  1972. 

i 

1 
1980. 

AFNC 

(TIME, 

A, 38) 

■""' 

i 

1 
1980. 

AFNC 

(TIME, 

A, 39) 

l 

11972. 

i 

1 
1980. 

■ 

AFNC 

(TIME, 

A, 40) 

1 1 972. 

i 

1 
1980. 

% 

AFNC 

(TIME, 

A, 41) 

i1972- 

i 

1 
1980. 

AFNC 

(TIME, 

A, 42) 

1985, 


.000011 


.000009 


.000007 


Employee/output   ratios    for  petroleum  and 
natural    gas   Sector    (26) 


0 


1985, 


.0000] 3 


.000012 


Employee/output   ratios    for  construction 
Sector    (1) 


0 


1985. 


.000033 


.000024 


.000020 


-   Employee/output   ratios    for   food   processing 
Sectors     (2,3,4,5) 


0 


1985. 


.000016 


.000013 


.000011 


-    Employee/output   ratios    for   forestry 
Sector    (25) 


0 


1985. 


.000030 


.000021 


.000017 


-    Employee/output   ratios    for   petroleum  refining 
Sector    (30) 


0 


1985. 


.000005 


.00000- 


.000004 


Employee/output    rabios    for    primary   metals 
Sector     (29) 


0 


1985, 


.000009 


.000007 


.000006 


Employee/output   ratios    for   other   manufacturing 
Sectors     (6,7,8,9,10) 


0 


1985, 


.000028 


.000023 


.000021 


-    Employee/output   ratios    for   railroad    trans- 
portation  Sector    (11) 


0 


11972 


1980. 


1985 


.000037 


.000022 


.000017 


C-20 


■ 
■ 


■ 
■ 


AFNC 


(TIME,    A, 43)    -    Employee/output    ratios    for   motor    freight 
Sector    (12) 


0 


1972, 


1980. 


1985. 


.000038 


.000033 


.000014 


AFNC 


(TIME,    A,  44)     -    Employee/output   ratios    for   conimunications 
Sector    (13) 


0 


1972. 


1980. 


1985. 


.000032 


.000021 


.000017 


AFNC 


(TIME,    A, 45)    -    Employee/output    ratios    for   utilities 
Sector     (14) 


0 


1972, 


1980. 


1985. 


.000014 


.000010 


.000008 


AFNC 


(TIME,    A, 46)    -    Employee/output   ratios    for  wholesale   and 
retail    trade   Sectors    (15,16) 


0 


1972. 


1980. 


1985. 


.000078 


.000058 


AFNC 


(TIME,    A, 4  7)     -    Employee/output    ratios    for    finance,    insurance, 
and   real    estate   Sector    (17) 


1 


0 


1972, 


1980. 


1985. 


.000017 


.000014 


.000013 


AFNC 


(TIME,    A, 48)    -    Employee/output   ratios    for    services 
Sectors     (18,19,20,22) 


0 


1972. 


1980. 


1985, 


.000118 


.000100 


.000092 


AFNC 


(TIME,    A, 49)    -    Employee/output   ratios    for    enterprizes 
Sectors    (21,22) 


0 


1972. 


1980. 


1985. 


0.0 


0.0 


0.0 


AFNC 


1972. 


AFNC 


(TIME,    A, 50)    -    Employee/output    ratios    for    agriculture 
Sector    (24) 

10  3  2 

1980.  1985.  .0000407  .0000305  .0000265 

(TIME,    A.., 62)    -    Price    for   residential    gas 


1978.  1980. 

.0000014  .0000016 

.000004  .0000045 


1 

1 

0 

9 

1970. 

1975. 

1976. 

1977. 

1985. 

1990. 

1995. 

.000001 

.000002 

.0000025 

.000003 

.0000035 

C-21 


1 

I 
■ 
■ 


AFNC  (TIME,  Z,63)  -  Cost  escalation  of  insulation 


1973. 

1974. 

1980. 

1985. 

1990. 

1995 

.677 

1.0 

1.0 

1.0 

1.0 

1.0 

AFNC  (TIME,  P;lo6'  64)  -  Baseline  energy  consumption  for  space 

heating 


1970. 

1975. 

1980. 

1985. 

1990. 

1995. 

135840000 

135840000 

130000000 

125000000 

123000000 

120000000 

1972, 


1972 


AFNC  (TIME,  P,00,51)  -  Federal  government  spending 


1980. 


1985. 


0 
.0000114 


.000091 


AFNC 


AFNC 


(TIME,    P-^q-^,52)    -    State   government   spending 

10  3  2 

1980.  1985. 

(TIME,    P45,5  3)    -    Labor/population   ratios     (WPR) 


,000084 


1972. 
.48 


0.0 
1.25 


1 

,3 

4.0 

i 

0.0 

1 

0.0 

1.25 

1975. 
.51 


1980. 
.515 


1985. 
.52 


1990. 


AFNC 


AFNC 


(TDO,    PCMO,54)    -    Increase    in   production   costs    as   oil 

supply    is    depleted 

115  2 

10000000000.  20000000000.  30000000000.  40000000000.  1.0 

1.5  2.0  3.0  5.0 

(OPADJ,    X,4)     -    Incremental    production   costs    of   oil 
supply    is   depleted 


0 


12 


2.51 

5.0 

5.0 


2.52 

6.5 

10.0 


2.53 

8.0 
15.0 


2.54 
10.0 
20.0 


3.0 

15.0 

25.0 


AFNC 


(TDG,  PCMG,55)  -  Increase  in  production  costs  as  gas 

supply  is  depleted 

115  2 

1000000000000.   2000000000000.   3000000000000.   4000000000000.   1.0 
1.5  2.5  4.0 


C-22 


s 

■ 
■ 


II 
II 


AFNC  (GPADJ,  X,8)  -  Incremental  cost  of  producing  gas 


0 

1 

1 

12 

2 

.000050 

.000051 

.000052 

.000053 

.000054 

.000080 

.000140 

.006200 

.000300 

.000450 

.000600 

.000900 

0.0 

5.0 

10.0 

15.0 

20.0 

25.0 

30.0 

35.0 

40.0 

45.0 

50.0 

55.0 

I 


AFNC  (X,  Y,3)  -  Oil  investment  expense 


0 


26 


1.0 

3.0 

5.0 

7.0 

9.0 

11.0 

13.0 

15.0 

17.0 

19.0 

21.0 

23.0 

25.0 

27.0 

29.0 

31.0 

33.0 

35.0 

37.0 

39.0 

41.0 

43.0 

45.0 

47.0 

49.0 

51.0 

.0080 

.0324 

.0570 

.0862 

.0803 

.0404 

.0269 

.0117 

.0117 

.0117 

.0110 

.OLIO 

.0110 

.0110 

.0110 

.0110 

.0110 

.0110 

.0110 

.0110 

.0110 

.0110 

.0110 

.0110 

.0110 

.0110 

.0110 

.0110 

I 


AFNC  (X,  W,l)  -  Oil  production 


26 


1.0 

3.0 

5.0 

7.0 

9.0 

11.0 

13.0 

15.0 

17.0 

19.0 

21.0 

23.0 

25.0 

27.0 

29.0 

31.0 

33.0 

35.0 

37.0 

39.0 

41.0 

43.0 

45.0 

47.0 

49.0 

51.0 

.00067 

.00337 

.00853 

.01572 

.02246 

.02583 

.02807 

.02876 

.02919 

.02876 

.02807 

.02583 

.02246 

.01797 

.01382 

.0114 

.00898 

.00786 

.00651 

.00517 

.00449 

.00337 

.00247 

.00112 

.00022 

.0 

AFNC  (X,  Y,2)  -  Oil  operating  expense 


I 
I 

1 


0 

0 

1 

26 

3 

1.0 

3.0 

5.0 

7.0 

9.0 

11.0 

13.0 

15.0 

17.0 

19.0 

21.0 

23.0 

25.0 

27.0 

29.0 

31.0 

33.0 

35.0 

37.0 

39.0 

41.0 

43.0 

45.0 

47.0 

49.0 

51.0 

.00168 

.00843 

.02133 

.03930 

.05615 

.06458 

.07018 

.07190 

.07298 

.07298 

.07018 

.07103 

.06738 

.05990 

.05115 

.04671 

.04191 

.03930 

.03581 

.03102 

.03143 

.02528 

.01976 

.01008 

.00220 

.0 

AFNC  (X,  Y,5)  -  Gas  production 


26 


1.0 

3.0 

5.0 

7.0 

9.0 

11.0 

13.0 

15.0 

17.0 

19.0 

21.0 

23.0 

25.0 

27.0 

29.0 

31.0 

33.0 

35.0 

37.0 

39.0 

41.0 

43.0 

45.0 

47.0 

C-23 


I 

I! 
I 
I 

I 
II 
II 
I 
II 
II 
II 
II 
II 
II 
II 
II 
II 
II 
II 


AFNC 

(X,  Y,5) 

-  Gas  production 

(cont'd) 

$9.0 

51.0 

.00067 

.00337 

.00853 

.01572 

.02246 

.02583 

.02807 

.02876 

.02919 

.02876 

.02807 

.02583 

.02246 

.01797 

.01382 

.0114 

.00898 

.00786 

.00651 

.00517 

.00449 

.00337 

.00247 

.00112 

.00022 

.0 

7AFNC     (X,    Y,7)     -    Gas    investment 


26 


1.0  . 

3.0 

5.0 

7.0 

9.0 

1.1.0 

13.0 

15.0 

17.0 

19.0 

21.0 

23.0 

25.0 

27.0 

29.0 

31.0 

33.0 

35.0 

37.0 

39.0 

41.0 

43.0 

45.0 

47.0 

49.0 

51.0 

.OOOOOOSO 

.00000324 

.0000057 

.00000862 

.00000808 

.00000404 

.00000269 

.00000117 

.00000117 

.00000117 

.00000110 

.00000110 

.00000110 

.00000110 

.00000110 

.00000110 

.00000110 

.00000110 

.00000110 

.00000110 

.00000110 

.000001 10 

.00000110 

.00000110 

.00000110 

.00000110 

7AFNC     (X,    Y,6)    -    Gas   operating    expense 


• 


0 

0 

1 

26 

3 

1.0 

3.0 

5.0 

7.0 

9.0 

11.0 

13.0 

15.0 

17.0 

19.0 

21.0 

23.0 

25.0 

27.0 

29.0 

31.0 

33.0 

35.0 

37.0 

39.0 

41.0 

43.0 

45.0 

47.0 

49.0 

51.0 

.0000000335 

.0000001685 

.000004265 

.000000786 

.000001123 

.000001369 

.000001488 

.000001524 

.000001547 

.000001524 

.000001684 

.000001808 

.000001797 

.000001797 

.000001797 

.000001710 

.000001616 

.000001572 

.000001497 

.000001448 

.000001572 

.000001348 

.000001112 

.000000560 

.000000121 

7AFNC  (TDO,  DRO,9)  -  Drilling  required/unit  of  oil  discovered 

as  supplies  are  depleted 


I 


0 

0 

1 

11 

2 

2000000000. 

2200000000. 

2400000000. 

2600000000. 

2800000000. 

3000000000 

3200000000. 

3400000000. 

3600000000. 

3S00000000. 

4000000000. 

.04 

.04 

.06 

.08 

.10 

.12 

.14 

.16 

.18 

100. 

AFNC  (TDG,  DRG,10)  -  Drilling  required/unit  of  gas  discovered 

as  supplies  are  depleted 


i 
* 


0 
14000000000. 
26,000000000. 
.0000015 
.00003 


0 
1600000000000. 
2800000000000. 
.000002 
.00010 


1 
1800000000000. 
3000000000000. 
.000004 
.100 


11 

2000000000000. 
3200000000000. 
.000008 


C-2  4 


2200000000000. 
3400000000000. 
.000015 


2400000000000. 

.0000010 

.00002 


1 
I 

I 
I 

II. 
I 

I 

I 

( 

■ 
■ 

II 

I 
I 

■I 

II 

I 
I 
I 
J 


AFNC    (TM,    CI,    K3)    -    Coal    mining   cost,  distribution    function    #14 
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AFNC  (PBBA,  A, 57)  -  Rate  insulation  upgraded  versus  payback 

period 
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AFNC  (TIME,  PRS,66)  -  Price  of  coal  for  industry 
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APPENDIX  D 
Price  Scenario  Input 

After  the  baseline  case  was  established,  variations  were 
made  in  some  assumptions.   The  model  then  computes  the  effects 
of  these  new  assumptions  on  natural/synthetic  gas  demand  in 
Montana.   Changes  in  assumptions  were  implemented  through 
changes  in  input  to  the  model.   This  appendix  lists  the  input 
changes  used  for  the  scenarios  in  this  report.   All  input 
values  not  listed  here  are  unchanged  from  those  used  to  es- 
tablish the  baseline  case. 

Case  1.   Firm  industrial  inter-fuel  substitution  (before  1981) 
Input : 

ISPRC  =  1 

Case  2 .   Firm  industrial  inter-fuel  substitution  (before  1981) 
and  extrapolated  industrial  inter-fuel  substitution 
(after  1980) .   Input: 
ISPRC  =  1 
ISPRiM  =  1 

Case  3.   Residential  conservation.   Input: 
ISCON  =  1 

Case  4.   Firm  industrial  inter-fuel  substitution  and  extrapo- 
lated industrial  inter-fuel  substitution  and  resi- 
dential conservation;  with  the  residential  gas  price 
increasing  to  $2.00/MCF  in  the  year  2000.   Industrial 
gas  increases  in  price  to  $1.50/MCF  in  the  year  2000. 
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Input : 

ISCON  =  1 

ISPRC  -  1 

ISPRM  =  1 

Function  62  (Price  of  residential  gas)  v/as 

modified  as  follows: 

Function  control  -  1,  1,  0,  6,  2 

Independent  variable  -  1975,  1976,  1977,  1980,  1985, 
(year)  2000 

Dependent  variable  -  1.4E-6,  1.6E-6,  1.6E-6,  1.8E-6, 
($/million  Btu)     1.8E-6,  2.0E-6 

Function  65  (Price  of  industrial  gas)  was 

modified  as  follows: 

P'unction  control  -  1,  1,  1,  6,  2 

Independent  variable  -  1975,  1976,  1977,  1980,  1985, 
(year)  2000 

Dependent  variable  -  .9E-6,  1.1E-6,  1.1E-6,  1.3E-6, 
($/million  Btu)     1.3E-6,  1.5E-6 

Case  5.   Firm  industrial  inter-fuel  substitution  and  extrapo- 
lated industrial  inter-fuel  substitution  and  resi- 
dential conservation;  with  the  residential  gas  price 
increasing  to  $3/MCF  in  the  year  2000.   Industrial 
gas  increases  in  price  to  $2.50/MCF  in  the  year  2000. 
Input: 

ISCON  =  1 

ISPRC  =  1 

ISPRM  =  1 

Function  62  (Price  of  residential  gas)  was 
modified  as  follows: 
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Function  control  -  1,  1,  0,  6,  2 

Independent  variable  -  1975,  1976,  1977,  1980,  1985, 
(year)  2000 

Dependent  variable  -  1.4E-6,  1.6E-6,  2.0E-6,  2.25E-6, 
($/million  Btu)     2.35E-6,  3.0E-6 

Function  65  (Price  of  industrial  gas)  was  modi- 
fied as  follows: 
Function  control  -  1,  1,  1,  6,  2 

Independent  variable  -  1975,  1976,  1977,  1980,  1985, 
(year)  2000 

Dependent  variable  -  .9E-6,  1.1E-6,  1.5E-6,  1.75E-6, 
($/million  Btu)     1.85  E-6,  2.5E-6 

Case  6.   Firm  industrial  inter-fuel  substitution  and  extrapo- 
lated industrial  inter-fuel  substitution  and  resi- 
dential conservation;  with  the  residential  gas  price 
increasing  to  $6.50/MCF  in  the  year  2000.   Industrial 
gas  increases  in  price  to  $6.00/MCF  in  the  year  2000. 
Input : 

I SCON  =  1 

ISPRC  =  1 

ISPRM  =  1 

Function  62  (Price  of  residential  gas)  was 

modified  as  follows: 

Function  control-  1,  1,  0,  6,  2 

Independent  variable  -  1975,  1976,  1977,  1980,  1985, 
(year)  2000 

Dependent  variable  -  1.4E-6,  1.6E-6,  2.0E-6,  3.5E-6, 
($/million  Btu)     4.5E-6,  6.5E-6 
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Function  65  (Price  of  industrial  gas)  was  modi- 
fied as  follows: 
Function  control  -  1,  1,  1,  6,  2 

Independent  variable  -  1975,  1976,  1977,  1980,  1985, 
(year)  2000 

Dependent  variable  -  .  9E-6,  1.1E-6,  1.5E-6,  3.0E-6, 
($/million  Btu)     4.0E-6,  6.0E-6 
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